Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 985,855 A→T 7.6% intergenic (‑146/+39) aspC ← / ← ompF aspartate aminotransferase, PLP‑dependent/outer membrane porin 1a (Ia;b;F)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913985,8550AT 7.6% 254.7 / 7.5 90intergenic (‑146/+39)aspC/ompFaspartate aminotransferase, PLP‑dependent/outer membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base A (44/38);  new base T (3/4);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTT  >  NC_000913/985707‑986003
                                                                                                                                                     |                                                                                                                                                    
cATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCaaa                                                                                                                                                     >  2:198076/1‑150 (MQ=255)
 aTGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAg                                                                                                                                                    <  1:691876/150‑1 (MQ=255)
  tGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGt                                                                                                                                                   >  2:1528948/1‑150 (MQ=255)
  tGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGt                                                                                                                                                   >  1:1685205/1‑150 (MQ=255)
   gACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAAACAGGACCAAAGt                                                                                                                                                   >  1:1565947/1‑150 (MQ=255)
    aCGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTcc                                                                                                                                                 <  1:962661/150‑1 (MQ=255)
      gatgTTCCATTATGGTTACAGAAGCGAAGTCCGCTATCAGGGTAACTGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTTTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTg                                                                                                                                               <  1:121577/147‑1 (MQ=255)
      gAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTg                                                                                                                                               <  1:185329/150‑1 (MQ=255)
         gTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACTTTGGTCCTGttt                                                                                                                                            <  1:1307257/150‑1 (MQ=255)
            ccATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACTTTGGTCCTGtttttt                                                                                                                                         <  1:81579/150‑1 (MQ=255)
             cATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACTTTGGTCCTGttttttt                                                                                                                                        <  2:375244/150‑2 (MQ=255)
              aTTATGGTTACAGAAGGGAAGTCCGCTATCATGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCg                                                                                                                                       <  1:232859/150‑1 (MQ=255)
              aTTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCg                                                                                                                                       <  1:658632/150‑1 (MQ=255)
               ttATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCgg                                                                                                                                      >  2:1417004/1‑150 (MQ=255)
                  tGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGTGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGTCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCAt                                                                                                                                   <  1:1452016/150‑1 (MQ=255)
                   ggTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCAtt                                                                                                                                  >  1:1214616/1‑150 (MQ=255)
                    gTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGt                                                                                                                                              >  2:887905/1‑137 (MQ=255)
                    gTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGt                                                                                                                                              <  1:887905/137‑1 (MQ=255)
                        cAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACa                                                                                                                             >  1:363045/1‑150 (MQ=255)
                          gAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACaaa                                                                                                                           >  1:117909/1‑150 (MQ=255)
                          gAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACaaa                                                                                                                           >  2:946133/1‑150 (MQ=255)
                             gggAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAagag                                                                                                                        <  1:520463/150‑1 (MQ=255)
                               gAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgt                                                                                                                      >  1:1682708/1‑150 (MQ=255)
                                aaGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCAAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtg                                                                                                                     <  1:1522425/150‑1 (MQ=255)
                                  gTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTTCTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTATCAAAGAGGtgtg                                                                                                                   <  1:490448/150‑1 (MQ=255)
                                    ccGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCt                                                                                                                 >  1:1665153/1‑150 (MQ=255)
                                    ccGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCt                                                                                                                 >  2:1186958/1‑150 (MQ=255)
                                    ccGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCt                                                                                                                 >  1:955221/1‑150 (MQ=255)
                                      gCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGAAGGATGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCGTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTAt                                                                                                               >  1:851321/1‑150 (MQ=255)
                                        tATCAGTGTAACGGGAGATTTACAAAAATCCAACTATAACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTa                                                                                                             <  2:209792/150‑1 (MQ=255)
                                        tATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCGAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTa                                                                                                             <  2:542100/150‑1 (MQ=255)
                                        tATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTa                                                                                                             >  2:1572239/1‑150 (MQ=255)
                                         aTCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGTAGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAg                                                                                                            <  1:151521/150‑1 (MQ=255)
                                            aGGGTTACTGTAGATTTACAAAATTCCAACTATTACTGATGAAAACGCTGGCTTTTTTTGCAAGACGTGAGTTTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAAc                                                                                                         <  2:340753/150‑1 (MQ=255)
                                            aGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAAc                                                                                                         <  1:198076/150‑1 (MQ=255)
                                              ggTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTg                                                                                                       >  2:1129238/1‑150 (MQ=255)
                                               gTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTATCAAAGAGGTGTGCTATTAGAACTgg                                                                                                      >  1:1551034/1‑150 (MQ=255)
                                               gTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTgg                                                                                                      <  2:637433/150‑1 (MQ=255)
                                               gTAAAGGGAGATTTTCAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTTGAAGGGAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTgg                                                                                                      <  2:302338/150‑1 (MQ=255)
                                                ttaCGGGAGTTTTTCAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAAGAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGt                                                                                                     <  1:1239840/148‑1 (MQ=255)
                                                     ggAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtg                                                                                                                     <  2:1465033/129‑1 (MQ=255)
                                                     ggAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtg                                                                                                                     >  1:1465033/1‑129 (MQ=255)
                                                     ggAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCGGCTTTTGCAAGCCGGGAGATTGCTCTGGAAGGTAT‑AAAAAAAACACGCCCAAAGTCCTGTTTTTTCGGCATTTAACAAAGCGGCGTGCTATTAGAACTGGTAAACg                                                                                                <  2:1169644/150‑1 (MQ=255)
                                                         aTTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATAc                                                                                            <  1:1417004/150‑1 (MQ=255)
                                                            tACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCca                                                                                         <  2:1230494/150‑1 (MQ=255)
                                                              cAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCcaca                                                                                       >  2:1697129/1‑150 (MQ=255)
                                                               aaaaTTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAg                                                                                      >  2:1534912/1‑150 (MQ=255)
                                                                aaaTTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGc                                                                                     >  1:635002/1‑150 (MQ=255)
                                                                  aTTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCaa                                                                                   <  1:946133/150‑1 (MQ=255)
                                                                        aCTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGgtgt                                                                             <  1:1528948/150‑1 (MQ=255)
                                                                        aCTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGgtgt                                                                             >  1:935888/1‑150 (MQ=255)
                                                                         cTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTgtc                                                                            >  2:490716/1‑150 (MQ=255)
                                                                          tATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTgtcg                                                                           <  2:955221/150‑1 (MQ=255)
                                                                                 gatgaAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAAc                                                                    >  1:1066458/1‑150 (MQ=255)
                                                                                       aaCGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGaa                                                                     <  1:1213512/143‑1 (MQ=255)
                                                                                       aaCGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTAGTCTGaa                                                                     >  2:1213512/1‑143 (MQ=255)
                                                                                         cGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGccc                                                            >  2:773086/1‑150 (MQ=255)
                                                                                           cAGGCTGTTTTTGCAAGCCGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATAAGAACTGCTAAACGATACCCACAGAAAAGGTGTCGTCTCAACCTACGGCCAg                                                          >  1:1611113/1‑150 (MQ=255)
                                                                                               ctctttttGCATGACGTGAGATTGCTCTGGAAGGAAT‑AAAAAAAACAGGACCACAGTCCGGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTttg                                                      <  1:1091109/147‑1 (MQ=255)
                                                                                               cTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTttg                                                      <  1:566665/150‑1 (MQ=255)
                                                                                                    tttGCATGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTc                                                 <  2:1123930/150‑1 (MQ=255)
                                                                                                     ttGCAAGTCGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCa                                                <  1:1350407/150‑1 (MQ=255)
                                                                                                      tgtaaGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAg                                               <  2:1570889/147‑1 (MQ=255)
                                                                                                       gCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTAGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGACCCTACGCCCAGTTTGTTGTCAGa                                              >  2:1001354/1‑150 (MQ=255)
                                                                                                         aaGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTAAGCCCAGTTTGTTGTCAGAAt                                            >  1:1002161/1‑150 (MQ=255)
                                                                                                          aGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATc                                           <  2:802193/150‑1 (MQ=255)
                                                                                                                  gaTTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGt                                   >  1:590974/1‑150 (MQ=255)
                                                                                                                       ctctGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtgatg                              >  2:746964/1‑150 (MQ=255)
                                                                                                                       ctctGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtgatg                              >  2:1353705/1‑150 (MQ=255)
                                                                                                                       ctctGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtgatg                              <  2:635002/150‑1 (MQ=255)
                                                                                                                       ctctGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtgatg                              >  2:265323/1‑150 (MQ=255)
                                                                                                                       ctctGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGATATTAGAACTGGTAAAAGATACCAACAGCAACGGTGTCGTCTGAACCTACGCCAAGTTTGTTGTCAGAATCGATCTGGTtgatt                              >  2:1598679/1‑149 (MQ=255)
                                                                                                                        tctGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtgatga                             <  2:117909/150‑1 (MQ=255)
                                                                                                                         ctGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTgatgat                            >  2:242133/1‑150 (MQ=255)
                                                                                                                            gAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTa                         >  1:1712227/1‑150 (MQ=255)
                                                                                                                                gTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGTATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCa                     >  1:887291/1‑150 (MQ=255)
                                                                                                                                gTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCa                     >  1:1543227/1‑150 (MQ=255)
                                                                                                                                 tat‑aAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCaa                    >  1:1247470/1‑150 (MQ=255)
                                                                                                                                  at‑aAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAAc                   <  1:488018/150‑1 (MQ=255)
                                                                                                                                  at‑aAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAAc                   <  2:404942/150‑1 (MQ=255)
                                                                                                                                   gccgaaaaaaCAGGACTTTGGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGtt                                  >  2:149550/5‑135 (MQ=255)
                                                                                                                                   gccgaaaaaaCAGGACTTTGGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGtt                                  <  1:149550/131‑1 (MQ=255)
                                                                                                                                    ccgaaaaaaCAGGACTGTGGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACa                  >  1:1305977/4‑150 (MQ=255)
                                                                                                                                     aaaaaaaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATAAACACAGCAACAGTGACGTCTGAACCTACGCCAACTATGTAGTCAGAATCGATCTGCTTAATTATGAAGTCAACAt                 >  2:1257451/1‑150 (MQ=255)
                                                                                                                                      aaaaaaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATa                >  2:1710536/1‑150 (MQ=255)
                                                                                                                                         aaaaCAGGACTTTGGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATTATGTAGTCAACATAGGt             >  1:1703146/1‑150 (MQ=255)
                                                                                                                                           aaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTgg           <  1:490716/150‑1 (MQ=255)
                                                                                                                                               ggACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACAt       >  1:396296/1‑150 (MQ=255)
                                                                                                                                                 aCCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGt     >  2:753067/1‑150 (MQ=255)
                                                                                                                                                  ccAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGatgatg                           <  2:337759/127‑1 (MQ=255)
                                                                                                                                                  ccAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGatgatg                           >  1:337759/1‑127 (MQ=255)
                                                                                                                                                    aaaGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGtttt  <  1:1572239/150‑1 (MQ=255)
                                                                                                                                                    aaaGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGtttt  <  2:1066458/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTAT‑AAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTT  >  NC_000913/985707‑986003

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑