Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,434,310 C→T 5.1% intergenic (‑86/+648) ynaE ← / ← ttcC cold shock protein, Rac prophage/pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,434,3100CT 5.1% 222.2 / 2.3 78intergenic (‑86/+648)ynaE/ttcCcold shock protein, Rac prophage/pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (35/39);  new base T (2/2);  total (37/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 = 9.43e-01

CCACTTTAGGGTGTCTTTTGATTTCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCC  >  NC_000913/1434198‑1434459
                                                                                                                |                                                                                                                                                     
ccACTTTAGGGTGTCTTTTGATTTCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTaga                                                                                                                  >  1:1566061/1‑150 (MQ=11)
ccACTTTAGGGTCTCTTATGATTTCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATAGTATGTCTGATTCAGGATATATTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCATAATTTTATTaga                                                                                                                  >  2:1605813/1‑150 (MQ=11)
         ggTGTCTTTTGATTTCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGAATGACTa                                                                                                         >  2:367580/1‑150 (MQ=17)
           tgtCTTTTGATTTCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTaaa                                                                                                       >  1:1558249/1‑150 (MQ=32)
            gtCTTTTGATTTCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGACAATCCCTAATTTTATTAGAGAATGACTaaaa                                                                                                      >  1:804545/1‑150 (MQ=17)
              cTTTTGATTTCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAAt                                                                                                    <  1:1849698/150‑1 (MQ=32)
                tttGATTTCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGc                                                                                                  <  2:835560/150‑1 (MQ=32)
                tttGATTTCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGc                                                                                                  <  1:1595682/150‑1 (MQ=32)
                  tGATTTCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTg                                                                                                <  2:1805713/150‑1 (MQ=32)
                    ttttCATGAATATATTACTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCATAATTTTATTAGAGCATGACTAAAAATGCTGaa                                                                                              <  1:1286977/149‑1 (MQ=12)
                    ttttCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATTTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGaa                                                                                              <  1:1293788/149‑1 (MQ=17)
                    aTTTCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGaa                                                                                              <  2:145859/150‑1 (MQ=32)
                    aTTTCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGaa                                                                                              <  2:1348357/150‑1 (MQ=32)
                     tttCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAat                                                                                             <  1:1577870/150‑1 (MQ=32)
                     tttCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCACTAATTTTATTAGAGCATGACTAAAAATGCTGAat                                                                                             >  1:802048/1‑150 (MQ=17)
                      ttCATGAATATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAata                                                                                            <  1:1784305/150‑1 (MQ=17)
                      ttCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAata                                                                                            >  1:955308/1‑150 (MQ=32)
                      ttCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAata                                                                                            >  1:1553107/1‑150 (MQ=32)
                         aTGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGa                                                                                         >  2:1628775/1‑150 (MQ=32)
                           gAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGAAATGTTGACTTATACATCGGTTTTGTCTGGTTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATa                                                                                       >  2:986902/1‑150 (MQ=17)
                                tatTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCATGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCTTGACTAAAAATGCTGAATATGATAAGGAg                                                                                  <  1:191378/150‑1 (MQ=17)
                                 atTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATTCATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGc                                                                                 <  1:731405/150‑1 (MQ=17)
                                 atTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGc                                                                                 <  1:457936/150‑1 (MQ=32)
                                   tCCTTAAAATATTATTGCTTTTCATTGCGATATTGTATGTCTGCTTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGa                                                                               <  2:175921/150‑1 (MQ=21)
                                    ccTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGaa                                                                              >  2:6398/1‑150 (MQ=35)
                                    ccTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGaa                                                                              >  1:558300/1‑150 (MQ=35)
                                    ccTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGaa                                                                              >  2:374461/1‑150 (MQ=35)
                                       tAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTg                                                                           <  2:1098477/150‑1 (MQ=37)
                                        aaaaTATTATTGATTTTCATTGCGATATTGTCTGTCTGATTCATGTTATGTAGTCTTTTACATCTGTTTTGTCTGGGTTATTGGATAAGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGa                                                                          <  1:1171009/150‑1 (MQ=11)
                                          aaTATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGGTATGGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGAtt                                                                        <  1:1422220/150‑1 (MQ=255)
                                                 ttGATTTTCATTGCGATCTTGTATGGCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTa                                                                 <  1:456314/150‑1 (MQ=255)
                                                   gATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGAGAAGTGATTATCAGTATg                                                               >  1:757487/1‑150 (MQ=17)
                                                    aTTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGc                                                              <  1:1379419/150‑1 (MQ=37)
                                                          aTTGCGATATTGTATGTCTGTTTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCa                                                        <  2:1535491/150‑1 (MQ=37)
                                                          aTTGCGATATTGTATGTCTGATTCAGGATATGTTGAGTTATACATCGGTTTTGTCTTGGTTATTGGATATGCCAATCCCTAATTTTCTTAGAGCATGACTAAAAATGCTGAATATGATAAGGATCGAAGTGATTATCAGTATGCTGTTCa                                                        <  2:616573/150‑1 (MQ=255)
                                                          aTTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCa                                                        >  2:1776918/1‑150 (MQ=37)
                                                            tGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTTGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCata                                                      >  2:1493880/1‑150 (MQ=37)
                                                            tGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCATAATTTTATTAGAGCATGACTAAAAATTCTGAATATGATAAGGAGCGAATTGATTATCAGTATGCTCTTCata                                                      >  1:33574/1‑150 (MQ=255)
                                                                atatTGTATTTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGTGTTATTCGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGc                                                  <  1:312034/150‑1 (MQ=255)
                                                                   ttGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTc                                               >  1:1454342/1‑150 (MQ=37)
                                                                    tGTATGTCTGATTCAGGCTATGATGACTTATACATCGGTTTTCTCTGGCTTCTTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGAAAAGGAGCGAAGTGATTATCAGTATGCTGTTCATAAAGCCTCg                                              <  2:1153404/150‑1 (MQ=11)
                                                                    tGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCg                                              >  2:570714/1‑150 (MQ=37)
                                                                    tGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCg                                              <  2:1583524/150‑1 (MQ=37)
                                                                    tGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATCAAAAGGAGCGAAGTGATTATCAGAATGCTGTTCATATAGCCTCg                                              <  2:1851831/150‑1 (MQ=255)
                                                                     gTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGa                                             <  1:1068082/150‑1 (MQ=37)
                                                                       aTGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAAt                                           <  1:1730191/150‑1 (MQ=37)
                                                                         gTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTTTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTa                                         >  1:264115/1‑150 (MQ=37)
                                                                         gTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTTTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTa                                         >  1:1099233/1‑150 (MQ=37)
                                                                         gTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCACTCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGATGTTCATATAGCCTCGAAtaa                                         >  2:158746/1‑148 (MQ=255)
                                                                         gTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTa                                         >  1:418245/1‑150 (MQ=37)
                                                                         gTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTa                                         <  1:1359408/150‑1 (MQ=37)
                                                                          tCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAg                                        >  1:1080134/1‑150 (MQ=37)
                                                                           ttGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGt                                       <  1:1079832/149‑1 (MQ=37)
                                                                              aTTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAAt                                    <  2:1600498/150‑1 (MQ=37)
                                                                                 taGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAAtgtg                                 >  2:651954/2‑150 (MQ=37)
                                                                                  cggATATGTTGACTTATACATCGGTTTTGTCTGGCTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAATATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATgtgt                                <  1:884325/149‑1 (MQ=255)
                                                                                  aGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATgtgt                                <  2:378506/150‑1 (MQ=37)
                                                                                   ggATATGTTGACTTATACATCGGTTTTGTTTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGtt                               <  2:1099233/150‑1 (MQ=37)
                                                                                   ggATATGTTGACTTATACATCGGTTTTGTTTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGtt                               <  2:264115/150‑1 (MQ=37)
                                                                                    gATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTTTta                              >  2:232603/1‑150 (MQ=37)
                                                                                    gATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTta                              >  2:1037223/1‑150 (MQ=37)
                                                                                       atGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTtatat                           >  2:730525/1‑150 (MQ=37)
                                                                                        tGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTtatata                          >  2:27373/1‑150 (MQ=37)
                                                                                          ttGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATg                        <  1:676300/150‑1 (MQ=37)
                                                                                           tGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGa                       <  2:1212783/150‑1 (MQ=37)
                                                                                           tGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGa                       >  1:808085/1‑150 (MQ=37)
                                                                                            taCTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGat                      <  1:778527/149‑1 (MQ=37)
                                                                                            gACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGat                      <  2:272777/150‑1 (MQ=37)
                                                                                            gACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGat                      >  2:536853/1‑150 (MQ=37)
                                                                                             aCTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGata                     <  2:655881/150‑1 (MQ=37)
                                                                                              cTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGatat                    >  1:194260/1‑150 (MQ=37)
                                                                                                tataCATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAg                  <  2:235504/150‑1 (MQ=37)
                                                                                                 ataCATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATTACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGt                 <  1:981344/150‑1 (MQ=37)
                                                                                                   aCATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTg               >  2:1058895/1‑150 (MQ=37)
                                                                                                    cATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCTTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGa              >  1:901696/1‑150 (MQ=37)
                                                                                                    cATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGAGTTATATATGATAAAGTTGa              >  1:235962/1‑150 (MQ=255)
                                                                                                     aTCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGAc             >  1:1483851/1‑150 (MQ=37)
                                                                                                       cGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACaa           >  2:513392/1‑150 (MQ=37)
                                                                                                                cTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTTTTATATATTATATAGTTGACAATTTTTATca  >  2:1382240/1‑149 (MQ=255)
                                                                                                                |                                                                                                                                                     
CCACTTTAGGGTGTCTTTTGATTTCATGAAGATATTCCTTAAAATATTATTGATTTTCATTGCGATATTGTATGTCTGATTCAGGATATGTTGACTTATACATCGGTTTTGTCTGGGTTATTGGATATGCCAATCCCTAATTTTATTAGAGCATGACTAAAAATGCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCC  >  NC_000913/1434198‑1434459

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