Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,152,553 G→C 5.6% M106I (ATG→ATC) ‡ metC → cystathionine beta‑lyase/L‑cysteine desulfhydrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,152,5530GC 5.6% 425.1 / 11.5 124M106I (ATG→ATC) ‡metCcystathionine beta‑lyase/L‑cysteine desulfhydrase
Reads supporting (aligned to +/- strand):  ref base G (58/59);  new base C (4/3);  total (62/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 138 total aligned reads displayed.

GGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCC  >  NC_000913/3152404‑3152702
                                                                                                                                                     |                                                                                                                                                     
ggACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCTATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTtgatg                                                                                                                                                       <  2:2103483/150‑1 (MQ=255)
ggACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTtgatg                                                                                                                                                       <  2:2395101/150‑1 (MQ=255)
ggACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTtgatg                                                                                                                                                       <  1:2322865/150‑1 (MQ=255)
 gACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTtgatga                                                                                                                                                      >  1:959681/1‑150 (MQ=255)
 gACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGACGTGTGAACTGGTAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTtgatga                                                                                                                                                      <  2:1736968/150‑1 (MQ=255)
  aCGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGAc                                                                                                                                                     <  2:1662793/150‑1 (MQ=255)
  aCGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGGGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGAc                                                                                                                                                     <  1:2333508/150‑1 (MQ=255)
    ggCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCa                                                                                                                                                   <  2:2257734/150‑1 (MQ=255)
    ggCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCa                                                                                                                                                   <  2:1441044/150‑1 (MQ=255)
    ggCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCa                                                                                                                                                   <  2:1834051/150‑1 (MQ=255)
           aaCGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGGTCATCAacacatg                                                                                                                                            <  1:812551/150‑4 (MQ=255)
              gTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGGTCATCAacacatgatc                                                                                                                                         <  1:1906616/150‑7 (MQ=255)
               ttAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGGTCATCAacacatgatcg                                                                                                                                        <  1:701476/150‑8 (MQ=255)
                 aaCCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGa                                                                                                                                      >  2:1239671/1‑150 (MQ=255)
                  aCCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGaa                                                                                                                                     >  2:1297903/1‑150 (MQ=255)
                   cccATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAAc                                                                                                                                    <  1:277103/150‑1 (MQ=255)
                   cccATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAAc                                                                                                                                    >  1:687828/1‑150 (MQ=255)
                    ccATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAAcc                                                                                                                                   >  1:177906/1‑150 (MQ=255)
                    ccATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAAcc                                                                                                                                   <  2:884102/150‑1 (MQ=255)
                     cATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCg                                                                                                                                  >  1:1229621/1‑150 (MQ=255)
                      aTTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGa                                                                                                                                 >  2:2043155/1‑150 (MQ=255)
                          ctcCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAAc                                                                                                                                    <  1:1054717/143‑1 (MQ=255)
                          ctcCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAAc                                                                                                                                    >  2:1054717/1‑143 (MQ=255)
                           tcCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAg                                                                                                                            <  2:1409545/150‑1 (MQ=255)
                           tcCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAg                                                                                                                            >  1:214916/1‑150 (MQ=255)
                            ccTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAgg                                                                                                                           <  2:1590531/150‑1 (MQ=255)
                            ccTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAgg                                                                                                                           <  1:1414416/150‑1 (MQ=255)
                            ccTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAgg                                                                                                                           >  2:804998/1‑150 (MQ=255)
                            ccTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAgg                                                                                                                           >  1:460781/1‑150 (MQ=255)
                                  aacaaGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCt                                                                                                                     >  2:507822/1‑150 (MQ=255)
                                   acaaGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTg                                                                                                                    <  2:2196490/150‑1 (MQ=255)
                                     aaGCGATGTGTGACCTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTa                                                                                                                  <  2:1192766/150‑1 (MQ=255)
                                      aGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAg                                                                                                                 <  2:2187988/150‑1 (MQ=255)
                                      aGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAg                                                                                                                 <  1:2464787/150‑1 (MQ=255)
                                       gCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGc                                                                                                                <  1:116526/150‑1 (MQ=255)
                                        cGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCa                                                                                                               >  1:1503628/1‑150 (MQ=255)
                                          aTGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCACCACCGCCTATGAACCGAGTCAGGATTTCTGTAGCaaa                                                                                                             >  1:1672670/1‑150 (MQ=255)
                                              gtgAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATcc                                                                                                         <  1:1297903/150‑1 (MQ=255)
                                              gtgAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATcc                                                                                                         >  2:686965/1‑150 (MQ=255)
                                                gAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTc                                                                                                       <  2:266877/150‑1 (MQ=255)
                                                  aCTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAg                                                                                                     <  2:280963/150‑1 (MQ=255)
                                                   cTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGc                                                                                                    >  1:587758/1‑150 (MQ=255)
                                                     ggAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCaa                                                                                                  <  1:1899422/150‑1 (MQ=255)
                                                        aGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACt                                                                                               >  2:1549307/1‑150 (MQ=255)
                                                             gcgcAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCg                                                                                          <  2:273460/150‑1 (MQ=255)
                                                             gcgcAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCg                                                                                          >  1:454375/1‑150 (MQ=255)
                                                              cgcAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGTCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGt                                                                                         <  1:1284935/150‑1 (MQ=255)
                                                              cgcAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGt                                                                                         <  2:1647855/150‑1 (MQ=255)
                                                               gcAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCa                                                                                                      <  1:1967659/136‑1 (MQ=255)
                                                               gcAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCa                                                                                                      >  2:1967659/1‑136 (MQ=255)
                                                                  ggCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACg                                                                                     <  2:1286871/150‑1 (MQ=255)
                                                                           cTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATgg                                                                            <  2:1622416/150‑1 (MQ=255)
                                                                            tATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGt                                                                           <  1:804998/150‑1 (MQ=255)
                                                                              tttCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGttt                                                                         >  1:1439828/1‑150 (MQ=255)
                                                                              tttCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGttt                                                                         <  2:1872843/150‑1 (MQ=255)
                                                                              tttCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGttt                                                                         >  2:2320439/1‑150 (MQ=255)
                                                                                  ccTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATc                                                                     <  1:470894/150‑1 (MQ=255)
                                                                                    tGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCg                                                                   >  1:633548/1‑150 (MQ=255)
                                                                                      cGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCAAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCt                                                                 <  2:775569/150‑1 (MQ=255)
                                                                                                gCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGcc                                                       <  2:2430326/150‑1 (MQ=255)
                                                                                                gCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGcc                                                       <  1:1384174/150‑1 (MQ=255)
                                                                                                gCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGcc                                                       >  1:2140928/1‑150 (MQ=255)
                                                                                                  ggTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACTAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGa                                                     >  2:1421705/1‑150 (MQ=255)
                                                                                                  ggTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGa                                                     >  1:1351862/1‑150 (MQ=255)
                                                                                                     tGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGatat                                                  >  1:143745/1‑150 (MQ=255)
                                                                                                     tGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGatat                                                  >  2:533493/1‑150 (MQ=255)
                                                                                                        tAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGt                                               <  2:643971/150‑1 (MQ=255)
                                                                                                         aaTTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGtt                                              <  1:2043155/150‑1 (MQ=255)
                                                                                                           ttCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTaa                                            >  2:820186/1‑150 (MQ=255)
                                                                                                           ttCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTaa                                            >  1:1505510/1‑150 (MQ=255)
                                                                                                             ccATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGc                                          >  2:330671/1‑150 (MQ=255)
                                                                                                               aTTCTTGCTTTTATCGAACCGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCAt                                        <  2:1687666/150‑1 (MQ=255)
                                                                                                                   ttGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGc                                    <  1:1661516/150‑1 (MQ=255)
                                                                                                                      cTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGc                                 >  2:1980552/1‑150 (MQ=255)
                                                                                                                      cTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGc                                 >  1:1682541/1‑150 (MQ=255)
                                                                                                                           aTCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAac                            >  1:1003638/1‑150 (MQ=255)
                                                                                                                           aTCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAac                            <  1:2320439/150‑1 (MQ=255)
                                                                                                                             cGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTTGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAacac                          <  1:2267672/150‑1 (MQ=255)
                                                                                                                             cGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAacac                          >  1:1066311/1‑150 (MQ=255)
                                                                                                                              gAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACt                         >  1:1905814/1‑150 (MQ=255)
                                                                                                                                aCAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTaa                       <  1:2017458/150‑1 (MQ=255)
                                                                                                                                     gCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCg                  <  1:799889/150‑1 (MQ=255)
                                                                                                                                     gCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCg                  <  2:1177888/150‑1 (MQ=255)
                                                                                                                                      ataggcggtgtTGGTCATCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAATGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCgt                 >  1:2062102/8‑150 (MQ=255)
                                                                                                                                       gATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCgtg                <  2:1883579/150‑1 (MQ=255)
                                                                                                                                       gATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCgtg                <  2:1492318/150‑1 (MQ=255)
                                                                                                                                        aTCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCgtgt               <  2:285817/150‑1 (MQ=255)
                                                                                                                                         ggcggtgtTGGTCATCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGtt              >  1:2369630/5‑150 (MQ=255)
                                                                                                                                            tgtgTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCt           >  2:609067/1‑150 (MQ=255)
                                                                                                                                             gtgtTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTg          >  2:1015041/1‑150 (MQ=255)
                                                                                                                                             gtgtTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTg          >  1:889924/1‑150 (MQ=255)
                                                                                                                                             gtgtTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGCGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTg          <  1:1312456/150‑1 (MQ=255)
                                                                                                                                              tgtTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTgg         <  2:1229621/150‑1 (MQ=255)
                                                                                                                                              tgtTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTgg         <  2:975958/150‑1 (MQ=255)
                                                                                                                                                 tgatgaCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAAt      >  1:2326703/1‑150 (MQ=255)
                                                                                                                                                 tgatgaCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCACTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAAt      >  2:540746/1‑150 (MQ=255)
                                                                                                                                                 tgatgaCCAACACAGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAAt      >  1:1458746/1‑150 (MQ=255)
                                                                                                                                                  gatgaCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATc     <  1:533493/150‑1 (MQ=255)
                                                                                                                                                  gatgaCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATc     >  2:799518/1‑150 (MQ=255)
                                                                                                                                                    tgaCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGc   >  2:906616/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCC  >  NC_000913/3152404‑3152702

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