Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,020,956 A→G intergenic (‑37/‑182) sulA ← / → sxy SOS cell division inhibitor/CRP‑S‑dependent promoter expression factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,020,9560AG98.8% 247.3 / ‑4.6 82intergenic (‑37/‑182)sulA/sxySOS cell division inhibitor/CRP‑S‑dependent promoter expression factor
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (43/37);  total (44/38)
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.12e-01

CACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATC  >  NC_000913/1020817‑1021090
                                                                                                                                            |                                                                                                                                      
cACTGATATGCCCCGCTGGAGTGTTGTCCGTAGAGACACGCGCAATTTTACGTGCTGCGGATTAGAACGTCGCAGATCGATGTGCATAGTCTCAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATgtgcag                                                                                                                                     <  1:111274/143‑1 (MQ=255)
cACTGATAAGCTGTGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATgtgcag                                                                                                                                     <  1:487249/143‑1 (MQ=255)
cACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATgtgcag                                                                                                                                     <  2:1118087/143‑1 (MQ=255)
 aCTGTTAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATgtgcagt                                                                                                                                    <  1:1049852/143‑1 (MQ=255)
 aCTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGCGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATgtgcagt                                                                                                                                    <  2:728254/143‑1 (MQ=255)
 aCTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATgtgcagt                                                                                                                                    >  2:309590/1‑143 (MQ=255)
    gacaaGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTGTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATgtgtagtaca                                                                                                                                 <  2:334390/140‑1 (MQ=255)
      tAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATc                                                                                                                               <  2:426374/143‑1 (MQ=255)
        aGCCCGGCTGTAGTGTTTTCCGTAGCGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCa                                                                                                                             >  1:1423922/1‑143 (MQ=255)
         gCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGGGCAGTACATCCAg                                                                                                                            >  1:576095/1‑143 (MQ=255)
          cccGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGCTACGGGATGGCTGTGCAGGACATACAGt                                                                                                                           >  1:409839/1‑143 (MQ=255)
           gctgCTGTAGTGTTTTCAGGAGAGACACGCGCAATTTTCCCTGCTGCGGATGAGAACGACGAAGCACGATGTGCATAGCCTGAAGTGTACATAATCAATACAGAACCTGTC‑CCTGACGGTATCGATGTCCAGTACATCCAGTg                                                                                                                          <  1:1037423/140‑1 (MQ=255)
                  tAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGCGCAGTACATCCAGTGCcaacaa                                                                                                                   >  1:981122/1‑143 (MQ=255)
                   aGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGCACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACaaa                                                                                                                  >  1:862019/1‑143 (MQ=255)
                    gtgtTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCCATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAg                                                                                                                 >  2:126110/1‑143 (MQ=255)
                    gtgtTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAg                                                                                                                 >  2:871231/1‑143 (MQ=255)
                         ttCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGCAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTGCTGTATGGGTGTGCAGTACATCCAGTGACCACAAAGATCaa                                                                                                            >  2:451427/1‑143 (MQ=255)
                         ttCAGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCaa                                                                                                            <  1:1008371/143‑1 (MQ=255)
                           ccGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAAcc                                                                                                          <  2:416038/143‑1 (MQ=255)
                             gTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCACTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAGCCAAggtacacccg                                                                                                        >  1:991134/1‑134 (MQ=255)
                                   acacGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTc                                                                                                  >  2:1441419/1‑143 (MQ=255)
                                      cgcgCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGa                                                                                               <  2:977641/143‑1 (MQ=255)
                                      cgcgCAATTTGACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGa                                                                                               <  2:44440/143‑1 (MQ=255)
                                       gcgcAATTTTTCTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAATGATCAACCCTATTTTCGGaa                                                                                              <  1:39124/143‑1 (MQ=255)
                                       gcgcAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGaa                                                                                              >  2:1417436/1‑143 (MQ=255)
                                       gcgcAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGaa                                                                                              >  1:1303042/1‑143 (MQ=255)
                                        cgcAATTTTACTTGCTGCGTATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGaaa                                                                                             >  2:866559/1‑143 (MQ=255)
                                        cgcAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGaaa                                                                                             <  1:546979/143‑1 (MQ=255)
                                         gcAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAAACAATCCAGCACCTGAG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGAACAACCCTATCTTCGGAAag                                                                                            >  1:145571/1‑143 (MQ=255)
                                          cAATTTTACTTGCTGGGGATGAGAACGACGAAGAACGATGTGCAGAGCCTGAAGTGTACATAATCAAGCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAaga                                                                                           <  2:858836/143‑1 (MQ=255)
                                            actttACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGc                                                                                         <  1:927842/141‑1 (MQ=255)
                                              tttACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCt                                                                                       <  1:309590/143‑1 (MQ=255)
                                               tttcTTGGTGCGTATGAGATCGACGAAGAACGATGTGCATACCCTTAACTGTATTTCATCAATCCATCTCCTGTG‑AGTTACTTTATGTATTTGCAGTACATCCAGTGTCAATAAATATCATCCCTATTTTCGGAAAGAGCCTc                                                                                      <  1:1281898/140‑1 (MQ=255)
                                                tACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCg                                                                                     <  2:530640/143‑1 (MQ=255)
                                                   ttGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCaa                                                                                  >  2:857250/1‑143 (MQ=255)
                                                       tgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATGTTCGGAAAGAGCCTCGCAAAttt                                                                              >  1:1044301/1‑143 (MQ=255)
                                                        gcGGAGGAGAACGACGAAGCACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAAtttt                                                                             <  2:185089/143‑1 (MQ=255)
                                                               agaACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTg                                                                      <  2:1419356/143‑1 (MQ=255)
                                                               agaACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTg                                                                      <  2:1224217/143‑1 (MQ=255)
                                                                 aaCGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTggg                                                                    >  2:220597/1‑142 (MQ=255)
                                                                 aaCGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGt                                                                    >  2:547830/1‑143 (MQ=255)
                                                                   cgacgaAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGa                                                                  >  2:368545/1‑143 (MQ=255)
                                                                   cgacgaAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGa                                                                  >  2:697159/1‑143 (MQ=255)
                                                                    tacaaagacggATGTGCGTAGCCGGGATTGTGCATTATCATTGCA‑CCCCCGTGTAGTTACTGTATGGATTTGCAGTAGATCCAGTGACAACAAAGATCAACCCTATTTTCGTAAATATCTTCGCAAATTTTGTCGTTGGTGAc                                                                 <  1:1146394/133‑1 (MQ=255)
                                                                      cgcagacCGATGTGCGTAGCCGGAAGTGTACATAATCAATCCAGCCCCTGGG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACgg                                                               <  1:623623/140‑1 (MQ=255)
                                                                      cgaAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACgg                                                               <  2:1423922/143‑1 (MQ=255)
                                                                              taTGTGCATATCCTGAAGTGTACATAATCAATCCATCCCTGGTG‑AGTTCCTGTATGGATGTGTAGCACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAACCTCACAAATTTTGTCGTTGGTGACGGGAAAACAt                                                       <  1:286250/142‑1 (MQ=255)
                                                                              gctgtgCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGCATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACAt                                                       <  1:1002095/141‑1 (MQ=255)
                                                                              gATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGCGCCTCGCACTTTTTGTCGTTGTTGAGAGGAAAACAt                                                       >  2:207147/1‑143 (MQ=255)
                                                                               aTGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATa                                                      >  2:1011417/1‑143 (MQ=255)
                                                                                   gCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAAtt                                                  >  2:1334392/1‑143 (MQ=255)
                                                                                    cATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGCATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAACTTTTGTCGTTGGTAACGGGAAAACATAACttc                                                 >  2:531974/1‑142 (MQ=255)
                                                                                         ccTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTGACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCtt                                            >  2:935425/1‑143 (MQ=255)
                                                                                          cTGAAGTGTACATGAGCAAGCCAGCCCCTGTG‑AGGTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAGCCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTg                                           <  1:44117/143‑1 (MQ=255)
                                                                                          cTGAAGTGTACATAATCAATCCAGCACATGAG‑TATTCCGCAAGGGAGGTGCAGTACATCCAGTGACAACTAAGATCAACCCTATCCCCGGAAAGAGCCTCGCAAAGGTTGTCGCAGGTTACGGGAAAACATAAATTAATCTTg                                           >  1:876432/1‑143 (MQ=255)
                                                                                           tGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAGACATAAATTAATCTTGc                                          >  2:1317013/1‑143 (MQ=255)
                                                                                           tGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCATATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCtttc                                          >  2:801888/1‑141 (MQ=255)
                                                                                                tgtACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGATAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCtt                                     >  2:1440162/1‑143 (MQ=255)
                                                                                                tgtACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATACATTAATCTAGcccact                                     >  1:14942/1‑140 (MQ=255)
                                                                                                 gtACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAGAACATAAATTAATCTTACCCCTTa                                    >  2:1331776/1‑143 (MQ=255)
                                                                                                  tACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTaa                                   <  1:483161/143‑1 (MQ=255)
                                                                                                  tACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTaa                                   >  2:583027/1‑143 (MQ=255)
                                                                                                  tACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTGTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACACAAATTAAGCTTGCCCCTTaa                                   >  1:34799/1‑143 (MQ=255)
                                                                                                   aCATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAACTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCGTTacg                                  >  1:1227315/1‑141 (MQ=255)
                                                                                                    cATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGGCGGGAAAACATAAATTAATCTTGCCCCTTAGGa                                 >  2:88870/1‑143 (MQ=255)
                                                                                                      tAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTCTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAAt                               >  2:574086/1‑143 (MQ=255)
                                                                                                        atcagtcCAGCCCCTGTG‑AGTTACTGCATCGATGTGCAGTACATCCAGGGACAACAAAGATCAACCCTATTTCCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATaa                             <  2:409839/143‑1 (MQ=255)
                                                                                                           actcCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGtt                          <  2:377275/141‑1 (MQ=255)
                                                                                                           actcCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCGCGCAAATGTTGTCGTTGGTGCCGGGAACACATCAATTAATCTTGCCCCTTAAGAATAAGtt                          <  1:880636/141‑1 (MQ=255)
                                                                                                           aatcCAGCCCCTGTG‑AGTTACTGTATGGATGCGCAGTACATCCAGTGACAACAACGTTCCACCATATTATTTGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGtt                          <  1:1175686/143‑1 (MQ=255)
                                                                                                            atcCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTCAGAATAAGTTg                         >  2:763183/1‑143 (MQ=255)
                                                                                                             tcCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTCATCTTGCCCCTTAAGAATAAGTTGc                        >  1:519670/1‑143 (MQ=255)
                                                                                                             tcCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGc                        >  1:1159316/1‑143 (MQ=255)
                                                                                                             tcCAGCCCCTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGc                        <  1:282100/143‑1 (MQ=255)
                                                                                                             cccAGCCCTTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGc                        <  1:287794/142‑1 (MQ=255)
                                                                                                                  ccccTGTG‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTAtt                   >  2:210018/1‑143 (MQ=255)
                                                                                                                       gtg‑AGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGt              <  2:437430/143‑1 (MQ=255)
                                                                                                                        tg‑AGTTACTGGATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCTTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTa             <  2:352129/143‑1 (MQ=255)
                                                                                                                         tgcgTTACTGGATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTa             <  1:1394339/140‑1 (MQ=255)
                                                                                                                           aggtACTGTATGGATGTGCAGTACATCAAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGt           <  2:19164/140‑1 (MQ=255)
                                                                                                                             ttACTGTATGGATGTGCAGTACATCCAGTGGCAACAAAGATCAACCATATTTCCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCGTATTTTCGTAGTTa         >  2:548697/1‑143 (MQ=255)
                                                                                                                             ttACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTa         >  2:708965/1‑143 (MQ=255)
                                                                                                                              ttgtcactGGATGTACTGCACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTaa        >  2:1405935/8‑143 (MQ=255)
                                                                                                                              tACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTaa        <  1:1417436/143‑1 (MQ=255)
                                                                                                                                cTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTCTATCCCGAAAGATCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACg      <  2:1435197/143‑1 (MQ=255)
                                                                                                                                 tGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCATATTTTCGGAAAGAGCCTCGCAAATTATGTCGTTGGTGACGGGAAAACATAAATTAATCTTACCCCTTAAGAATAAGTTGCCTCTTTTCAAAGATAAcag     >  1:492900/1‑141 (MQ=255)
                                                                                                                                  gTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGa    <  2:1360056/143‑1 (MQ=255)
                                                                                                                                   tATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGAt   <  2:290452/143‑1 (MQ=255)
                                                                                                                                    aTGGATGTGCAGTACATCCATTGGCGACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATGTTGGCGTTGGTGACGGGAAAACATAAATTAATCTTGCCGCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATc  <  1:631778/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                      
CACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTG‑AGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATC  >  NC_000913/1020817‑1021090

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