Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,378,977 T→A 6.2% intergenic (‑33/+30) ampC ← / ← frdD penicillin‑binding protein; beta‑lactamase, intrinsically weak/fumarate reductase (anaerobic), membrane anchor subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,378,9770TA 6.2% 394.5 / 10.9 138intergenic (‑33/+30)ampC/frdDpenicillin‑binding protein; beta‑lactamase, intrinsically weak/fumarate reductase (anaerobic), membrane anchor subunit
Reads supporting (aligned to +/- strand):  ref base T (66/63);  new base A (7/2);  total (73/65)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

Only 100 of 138 total aligned reads displayed.

GTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGGT  >  NC_000913/4378830‑4379125
                                                                                                                                                   |                                                                                                                                                    
gttgCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAg                                                                                                                                                    <  2:712239/150‑1 (MQ=255)
gttgCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAg                                                                                                                                                    <  2:577116/150‑1 (MQ=255)
     tctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGGCGTTTTGAACATAGGGGCTGGTTTCCATACAAAACGGCCCGCCATAGGCggg                                                                                                                                               <  2:1178712/150‑1 (MQ=255)
      ctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGc                                                                                                                                              >  2:391983/1‑150 (MQ=255)
      ctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGc                                                                                                                                              >  2:392785/1‑150 (MQ=255)
       tataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGCGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGcc                                                                                                                                             <  2:1583807/150‑1 (MQ=255)
       tataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGcc                                                                                                                                             <  2:1561242/150‑1 (MQ=255)
       aattaGCGGGGTAATTGTGCGTTGCACAATATGGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGcc                                                                                                                                             <  2:1568181/146‑1 (MQ=255)
        ataAGCTGGGTTATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCg                                                                                                                                            <  2:223614/150‑1 (MQ=255)
        ataAGCGTGGTATTTGTGCGATGCACAATATCGTTGATTTGTTGTGCGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCCCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCg                                                                                                                                            <  1:81763/150‑1 (MQ=255)
        ataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTATGGCGGGCCg                                                                                                                                            <  2:507018/150‑1 (MQ=255)
          aaGCGGGGTAATTGCGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTATGGCGGGCCgtt                                                                                                                                          <  1:474325/150‑3 (MQ=255)
           cgCGGGGTAATTGTGCGATGCACAATATCGGTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGGCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAt                                                                                                                                         <  2:1695462/149‑1 (MQ=255)
            gctgggTATTTGTGCGATGCACATTATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCATGAGGCGGGAATTAATAAGGCGCAGAGCGTCGTGTTGAACATAGGGTCTGGTATCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                                                                                                                                        <  1:343292/147‑1 (MQ=255)
            gCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                                                                                                                                        <  2:1654857/150‑1 (MQ=255)
             cgctgTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAttt                                                                                                                                       <  2:199623/146‑1 (MQ=255)
             caacgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGAGGGCCGGAttt                                                                                                                                       >  2:136821/8‑150 (MQ=255)
             cGGGGTAATTGTGCGATGCACAATATGGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAttt                                                                                                                                       <  2:1632147/150‑1 (MQ=255)
              aacgataTTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTa                                                                                                                                      >  2:585892/7‑150 (MQ=255)
                cgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                                                                                                                                        <  2:1226679/142‑1 (MQ=255)
                cgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                                                                                                                                        <  2:261989/142‑1 (MQ=255)
                cgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                                                                                                                                        >  1:261989/5‑146 (MQ=255)
                cgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                                                                                                                                        >  1:1226679/5‑146 (MQ=255)
                cgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACa                                                                                                                                    >  2:28118/5‑150 (MQ=255)
                        tgCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAAATAATAAGGCGCAGAGCGTCGTTTTGAACATAGGTTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACAATGGCCAt                                                                                                                            >  2:1513552/1‑150 (MQ=255)
                          cGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGCTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGc                                                                                                                          <  1:286666/150‑1 (MQ=255)
                           gATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCg                                                                                                                         <  1:512196/150‑1 (MQ=255)
                             tGCACAATATCGTTGATTTGTTGCGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGtt                                                                                                                       <  2:991218/150‑1 (MQ=255)
                             tGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGtt                                                                                                                       <  2:1182589/150‑1 (MQ=255)
                             tGCACAATAACGATGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACAGCCCGCCATAGGCGGGCCGGATTTACATTCgcgttgcgtt                                                                                                                       >  2:624957/1‑150 (MQ=255)
                              gCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGAGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACAATGGCGATGCGtaa                                                                                                                      >  1:141954/1‑148 (MQ=255)
                                acaATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGa                                                                                                                    >  2:141338/1‑150 (MQ=255)
                                   atatCGTTTATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTg                                                                                                                 <  2:193398/150‑1 (MQ=255)
                                   atatCGTTGATTTGTTGCGGGGCAGCAACTGTGGAGCAAGAGGCGGTAATTAAAAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTg                                                                                                                 <  1:362161/150‑1 (MQ=255)
                                       cGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAAc                                                                                                             >  2:1197610/1‑150 (MQ=255)
                                          tGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAACCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGac                                                                                                          <  1:1458734/150‑1 (MQ=255)
                                           taTTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGaca                                                                                                         <  1:809801/149‑1 (MQ=255)
                                            aTTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGTTTGTAACGacac                                                                                                        >  2:697149/1‑150 (MQ=255)
                                             tttGTTGCGGGGCCGCAAATGGGGAGCAAGAGGCGGTAATTAATAAGGCGCTGAGCGGCGTTTTGATCATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACAcc                                                                                                       <  2:1053541/150‑1 (MQ=255)
                                                gttgAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCAGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTTTAACGACAccttt                                                                                                    >  2:327539/1‑147 (MQ=255)
                                                  tgAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCa                                                                                                  >  1:566975/1‑150 (MQ=255)
                                                     ggggCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCg                                                                                               <  2:720706/150‑1 (MQ=255)
                                                       ggCATCAAATGTGGATCAAGAGGCTGTAATTAATAAGGCGCAGAGCGTCGTTTTTAACATAGGGTCTGGTTTCCATCCAAAACGGCCCGCCATTGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTg                                                                                             <  2:1444901/150‑1 (MQ=255)
                                                         cagcaAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATATGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGAAACCAATCAGCGTGaa                                                                                           >  1:535798/1‑149 (MQ=255)
                                                           gcaAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACaa                                                                                         >  2:1731295/1‑150 (MQ=255)
                                                            caAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCTGATTTACATTGTCGATGCGTTATATTGTAACGACACCAATCAGCGTGACAAc                                                                                        >  1:1127837/1‑150 (MQ=255)
                                                            caAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAAc                                                                                        >  1:681092/1‑150 (MQ=255)
                                                                tgtgGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTc                                                                                    >  1:62564/1‑150 (MQ=255)
                                                                 gtgGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCa                                                                                   >  2:1197343/1‑150 (MQ=255)
                                                                  tgGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAg                                                                                  >  1:735632/1‑150 (MQ=255)
                                                                   ggAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAgg                                                                                 <  1:585892/150‑1 (MQ=255)
                                                                     aGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGAt                                                                               >  1:1120794/1‑150 (MQ=255)
                                                                        aaGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATagc                                                                            <  1:1365477/150‑1 (MQ=255)
                                                                         agagGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATagca                                                                           >  1:155161/1‑150 (MQ=255)
                                                                            ggCGGTAATTAATAATGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGcc                                                                        <  2:567051/150‑1 (MQ=255)
                                                                            ggCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGcc                                                                        <  2:978262/150‑1 (MQ=255)
                                                                                gTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGTTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGAGGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCCGGATAGCAGCCAGAc                                                                    >  2:53184/1‑150 (MQ=255)
                                                                                 taattaatAAGGCGCAGAGCGTCGTTTTGAACATAGGTGCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGAcc                                                                   <  1:1297772/150‑1 (MQ=255)
                                                                                 taattaatAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGAcc                                                                   <  1:1197610/150‑1 (MQ=255)
                                                                                  aattaatAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCg                                                                  >  2:500309/1‑150 (MQ=255)
                                                                                    ttaatAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTa                                                                >  1:1340438/1‑150 (MQ=255)
                                                                                     taataaGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGTATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTat                                                               >  2:14696/1‑149 (MQ=255)
                                                                                     taataaGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACCTTGGCCATGCGTTATATTGTAACGACAACAATCAACGTGAAAACAGTCAGGATAGCAGCCAGACCGTat                                                               >  1:1271905/1‑149 (MQ=255)
                                                                                     taataaGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGCCACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAg                                                               >  2:1258452/1‑150 (MQ=255)
                                                                                       ataaGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGaa                                                             <  1:91106/150‑1 (MQ=255)
                                                                                         aaGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGaaaa                                                           <  2:1676980/150‑1 (MQ=255)
                                                                                          aGGCGCAGAGCGTCGTTTTGAACATAGGGCCTGGTGTCCATACAAAACGGCCCGCCATAGGCGGGACGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAAc                                                          <  1:872995/150‑1 (MQ=255)
                                                                                           ggCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAAcc                                                         >  1:1034256/1‑150 (MQ=255)
                                                                                                agagCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCAttt                                                    <  1:392785/150‑1 (MQ=255)
                                                                                                 gagCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAACCCCAGTTg                                                   >  1:764970/1‑150 (MQ=255)
                                                                                                 gagCGTCGTTTTGAACATAGGGTCTGGGTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTg                                                   <  2:836895/150‑1 (MQ=255)
                                                                                                       cgtTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGAAGAAAACCCATTTGCCCGCa                                             >  2:1341581/1‑150 (MQ=255)
                                                                                                         ttttGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAACCCCATTTGCCCGCAgg                                           >  1:828339/1‑150 (MQ=255)
                                                                                                            tGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTAc                                        >  2:1276572/1‑150 (MQ=255)
                                                                                                               aCATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGAGAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTg                                     <  1:130091/150‑1 (MQ=255)
                                                                                                                cATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTgg                                    >  2:804525/1‑150 (MQ=255)
                                                                                                                    gggTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGAttt                                <  2:929786/150‑1 (MQ=255)
                                                                                                                     ggTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGAtttt                               >  1:1492604/1‑150 (MQ=255)
                                                                                                                       tCTGGTTTCCTTACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATCCGTTAGATTGTAACGACACCAATCAGCGGGACAAGTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCa                             <  1:323106/150‑1 (MQ=255)
                                                                                                                        cTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAg                            <  2:728371/150‑1 (MQ=255)
                                                                                                                          ggTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGAt                          >  2:189838/1‑150 (MQ=255)
                                                                                                                          ggTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTAAATTGGCGATGCGTTATATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAAAAAACCCATTTTCCCGCAGGTAAGTGGATTTTAAGAt                          >  2:1524218/1‑150 (MQ=255)
                                                                                                                           gTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGaac                         >  2:639962/1‑148 (MQ=255)
                                                                                                                             ttCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGt                       <  1:1025392/150‑1 (MQ=255)
                                                                                                                                  tACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTACCGCCACCAATCAGCGTGAGACCTGTCAGGATTGGAGCCAGACCGTAGCAAACCCATGTGCCCGCAGGGAGGTGGATTTACAGATCGTGCTTc                  >  2:102562/1‑150 (MQ=255)
                                                                                                                                    cAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATcgc                >  1:788908/1‑150 (MQ=255)
                                                                                                                                    cAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATcgc                <  1:811360/150‑1 (MQ=255)
                                                                                                                                    aaatccGGCCCGCCTATGGCGGGCCGTATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATcgc                >  1:764323/6‑150 (MQ=255)
                                                                                                                                      atccGGCCCGCCTATGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATGAGCGTGACAACTGTCAGGTTAGCAGGCAGACCGTAGCAAACCCATTTGCCAGCAGGTACGTGGATTTTCAGATCGTGCATGGCgt              >  1:1047102/4‑150 (MQ=255)
                                                                                                                                      atccGGCCCGCCTATGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTAAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATAGTGCATCGCgt              >  2:936578/4‑150 (MQ=255)
                                                                                                                                        aCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtgg            <  2:750292/150‑1 (MQ=255)
                                                                                                                                         cGGCCCGCCTATGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGGACGTGGATTTTCAGATCGTGCATCGCgtggt           >  1:1192245/1‑150 (MQ=255)
                                                                                                                                         cGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggt           >  2:1073678/1‑150 (MQ=255)
                                                                                                                                          ggCCCGCCTATGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggtg          >  1:1554888/1‑150 (MQ=255)
                                                                                                                                          ggCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggtg          <  2:798342/150‑1 (MQ=255)
                                                                                                                                           gCCCGCCTATGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGc         >  1:1683107/1‑150 (MQ=255)
                                                                                                                                           gCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGc         >  1:598395/1‑150 (MQ=255)
                                                                                                                                             ccGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCAt       >  2:1054564/1‑150 (MQ=255)
                                                                                                                                              cGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATa      <  2:789085/150‑1 (MQ=255)
                                                                                                                                                  aTAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGgt  <  1:5395/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGGT  >  NC_000913/4378830‑4379125

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: