Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,439,399 A→C S37S (TCT→TCG)  yhdN ← DUF1992 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,439,3990AC98.6% 246.6 / ‑3.5 69S37S (TCT→TCG) yhdNDUF1992 family protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (28/40);  total (28/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

TGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCA  >  NC_000913/3439170‑3439625
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                  
tGTCAGCATAATCGCCGCATAAAAATTCGGTACTCAATCCAGCTTTTCGCAGATTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAGTTCCGGTGGCACGTGCGAatca                                                                                                                                                                                                                              <  2:385519/236‑1 (MQ=255)
    aGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAatcat                                                                                                                                                                                                                             <  1:699763/233‑1 (MQ=255)
    aGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAatcat                                                                                                                                                                                                                             >  2:699763/1‑233 (MQ=255)
          aTCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAatcatcatc                                                                                                                                                                                                                         <  1:424780/231‑1 (MQ=255)
          aTCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATcc                                                                                                                                                                                                                        <  2:362397/232‑1 (MQ=255)
          aTCGACGCCTAAAAAGTCGTTAACAAAGCCAGCTTGACGCAGCTTCAATTCCAGTAATGACAATCGACGGCGAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATcc                                                                                                                                                                                                                        <  2:140941/232‑1 (MQ=255)
            cgccgcGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCCTTCTTCAGCAAGCGATACACCGCACGGAATTCCGGTGGCACGTGCGAATCATCCTCCAgaaac                                                                                                                                                                                                                  >  2:320826/1‑232 (MQ=255)
                          tCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCt                                                                                                                                                                                                    <  1:720531/236‑1 (MQ=255)
                           cGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTa                                                                                                                                                                                                   <  2:39828/236‑1 (MQ=255)
                                  cAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCa                                                                                                                                                                                            <  1:518568/236‑1 (MQ=255)
                                  cAAGCCAGTTTGTAGCAGCTTCAAATCCAATAATGACAATCTACTGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCa                                                                                                                                                                                            <  2:385921/236‑1 (MQ=255)
                                     tCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGCATCAATGGTTCGCCGCTACGTGCCAgg                                                                                                                                                                                          >  2:528311/1‑235 (MQ=255)
                                             gTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCaaa                                                                                                                                                                                  >  2:162363/1‑235 (MQ=255)
                                             gTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAAc                                                                                                                                                                                 <  1:162363/236‑1 (MQ=255)
                                              tCGCAGCTTCAATTCCAGTAAGTACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATTCCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACt                                                                                                                                                                                <  1:336639/236‑1 (MQ=255)
                                                    cTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAAc                                                                                                                                                                                 >  1:388444/1‑229 (MQ=255)
                                                           tCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTg                                                                                                                                                                   <  1:623851/236‑1 (MQ=255)
                                                                    gACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGc                                                                                                                                                          >  2:668895/1‑236 (MQ=255)
                                                                    gACAAGCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCtt                                                                                                                                                              <  1:591089/232‑1 (MQ=255)
                                                                          cGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATg                                                                                                                                                    <  2:221319/236‑1 (MQ=255)
                                                                            aCGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATgg                                                                                                                                                   >  1:702736/1‑235 (MQ=255)
                                                                                cTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATgg                                                                                                                                                   >  2:51031/1‑231 (MQ=255)
                                                                                    ccTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGcc                                                                                                                                          <  1:51031/236‑1 (MQ=255)
                                                                                    ccTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGcc                                                                                                                                          <  1:649617/236‑1 (MQ=255)
                                                                                      ttttGATGTTGCTGTTCATCGTTACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTACCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCa                                                                                                                                        <  2:702736/234‑1 (MQ=255)
                                                                                      tctTGATAGTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCGCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCa                                                                                                                                        <  1:358742/236‑1 (MQ=255)
                                                                                               tGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGct                                                                                                                                               <  1:631286/220‑1 (MQ=255)
                                                                                               tGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGct                                                                                                                                               >  2:631286/1‑220 (MQ=255)
                                                                                                   gATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAAc                                                                                                                           <  1:668895/236‑1 (MQ=255)
                                                                                                          gTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATaa                                                                                                                    <  1:487580/236‑1 (MQ=255)
                                                                                                           tGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATaaa                                                                                                                   <  1:371777/236‑1 (MQ=255)
                                                                                                              cGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAatcatcatc                                                                                                                                                                                                                         <  1:763842/131‑1 (MQ=255)
                                                                                                              cGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAatcatcatc                                                                                                                                                                                                                         >  2:763842/1‑131 (MQ=255)
                                                                                                                  tACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCa                                                                                                                                                                                            <  2:311494/156‑1 (MQ=255)
                                                                                                                  tACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCa                                                                                                                                                                                            >  1:311494/1‑156 (MQ=255)
                                                                                                                       ttGAGAATATACAAAAGCTGAATTGCTTCGCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGa                                                                                                       <  2:86305/236‑1 (MQ=255)
                                                                                                                        tGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATcc                                                                                                                                                                                                                        <  2:93120/122‑1 (MQ=255)
                                                                                                                        tGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATcc                                                                                                                                                                                                                        >  1:93120/1‑122 (MQ=255)
                                                                                                                            aaTATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTgcgc                                                                                                                                                                <  2:145153/174‑1 (MQ=255)
                                                                                                                            aaTATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTgcgc                                                                                                                                                                >  1:145153/1‑174 (MQ=255)
                                                                                                                               atCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTc                                                                                                                                  >  2:253738/1‑201 (MQ=255)
                                                                                                                               atCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTc                                                                                                                                  <  1:253738/201‑1 (MQ=255)
                                                                                                                                 ccAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGGGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATGCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACaaaaa                                                                                             <  1:528311/236‑1 (MQ=255)
                                                                                                                                    gAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAagag                                                                                          <  1:371790/236‑1 (MQ=255)
                                                                                                                                         tGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACaaa                                                                                               <  2:605497/226‑1 (MQ=255)
                                                                                                                                         tGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAg                                                                                     <  1:600086/236‑1 (MQ=255)
                                                                                                                                          gAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGa                                                                                                       <  2:643684/217‑1 (MQ=255)
                                                                                                                                          gAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGa                                                                                                       >  1:643684/1‑217 (MQ=255)
                                                                                                                                            aTTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCACGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTa                                                                                  <  2:454551/236‑1 (MQ=255)
                                                                                                                                                cTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGt                                                                              >  2:225839/1‑236 (MQ=255)
                                                                                                                                                     ctcCGTTGCTCAAGTTCTGGCGGTAAGCACCCGGCATTCTTCAGCAAGCGATTCCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATa                                                                         <  1:342705/236‑1 (MQ=255)
                                                                                                                                                        cGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGt                                                                      <  1:48552/236‑1 (MQ=255)
                                                                                                                                                             cTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGa                                                                 >  1:687528/1‑236 (MQ=255)
                                                                                                                                                                  gTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGAtacta                                                            >  1:410169/1‑236 (MQ=255)
                                                                                                                                                                       ggCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACaaa                                                                                               >  1:437007/1‑196 (MQ=255)
                                                                                                                                                                       ggCGGTAAGCAACCGGCATTCTTCAGCAAGAGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACaaa                                                                                               <  2:437007/196‑1 (MQ=255)
                                                                                                                                                                              aGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATcc                                                                                                                                                                                                                        <  2:663837/68‑1 (MQ=255)
                                                                                                                                                                              aGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATcc                                                                                                                                                                                                                        >  1:663837/1‑68 (MQ=255)
                                                                                                                                                                                          tcttcAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAagag                                                                                          <  1:592233/182‑1 (MQ=255)
                                                                                                                                                                                          tcttcAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAagag                                                                                          >  2:592233/1‑182 (MQ=255)
                                                                                                                                                                                                          ccccGAACTTAATGCCGGTGGAACGTGCAAATCATCATCCAGTATCAATGGTTCGCCGTTACCTGCCAGGTTATCAAACTCCCCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGt                     <  2:492730/235‑1 (MQ=255)
                                                                                                                                                                                                               cACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTc               >  1:458405/1‑236 (MQ=255)
                                                                                                                                                                                                               cACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTc               >  1:17785/1‑236 (MQ=255)
                                                                                                                                                                                                               cACGTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGGTTTTTTACGGTGCGTc               >  1:492730/1‑236 (MQ=255)
                                                                                                                                                                                                                  gTAATTCCGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTCTTTTTACGTTGCTTCAg             >  2:570045/1‑235 (MQ=255)
                                                                                                                                                                                                                         cGGTGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTActct     >  2:250200/1‑236 (MQ=255)
                                                                                                                                                                                                                            tGGCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgca  >  1:132798/1‑236 (MQ=255)
                                                                                                                                                                                                                              gCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGa                                                                                                       >  2:493088/1‑133 (MQ=255)
                                                                                                                                                                                                                              gCACGTGCGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGa                                                                                                       <  1:493088/133‑1 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                  
TGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCGTGACGGATACCTTTGAGAATATCCAGAAGCTGAATTGCTTCTCTCCGTTGCTCAAGTTCTGGCGGTAAGCAACCGGCATTCTTCAGCAAGCGATACCCCGCACGTAATTCCGGTGGCACGTGAGAATCATCATCCAGTATCAATGGTTCGCCGCTACCTGCCAGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCA  >  NC_000913/3439170‑3439625

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑