Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,362 Δ1 bp 5.4% intergenic (+156/+16) ydbL → / ← feaR DUF1318 family protein/transcriptional activator for tynA and feaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,446,3620A. 5.4% 309.4 / 9.5 92intergenic (+156/+16)ydbL/feaRDUF1318 family protein/transcriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base A (47/40);  new base . (2/3);  total (49/43)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

ACCGCAATATCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGG  >  NC_000913/1446214‑1446510
                                                                                                                                                     |                                                                                                                                                    
aCCGCAATATCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACa                                                                                                                                                      <  1:385230/149‑1 (MQ=255)
   gCAATATCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGc                                                                                                                                                   <  1:400753/149‑1 (MQ=255)
      atatCGTGATAACGCTACAGTCGAAGGCGATTGATCTATTTTCCTGAAACAATGTTAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAAc                                                                                                                                                <  7:327945/149‑1 (MQ=255)
      atatCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAAc                                                                                                                                                <  8:63880/149‑1 (MQ=255)
      atatCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTCAATATTCAAAAACTCCTGTCAAATTGACTTTTGCCATGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAAc                                                                                                                                                <  7:327949/149‑1 (MQ=255)
       tatCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACg                                                                                                                                               >  4:391266/1‑149 (MQ=255)
        atCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGc                                                                                                                                              >  1:205470/1‑149 (MQ=255)
        atCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGc                                                                                                                                              >  3:388367/1‑149 (MQ=255)
         tCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCTTTGCATGTAC‑TGCAACGc                                                                                                                                              >  3:29074/1‑148 (MQ=255)
         tCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGc                                                                                                                                              <  3:391266/148‑1 (MQ=255)
           gggATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCAGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCaa                                                                                                                                           <  1:419561/149‑1 (MQ=255)
            ggATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAt                                                                                                                                          >  4:214363/1‑149 (MQ=255)
             gATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAtt                                                                                                                                         <  2:395514/149‑1 (MQ=255)
              aTAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAtta                                                                                                                                        <  4:372481/149‑1 (MQ=255)
              aTAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAtta                                                                                                                                        <  3:317339/149‑1 (MQ=255)
                aaCGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagt                                                                                                                                      <  7:453605/149‑1 (MQ=255)
                   gCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagttag                                                                                                                                   >  3:379854/1‑149 (MQ=255)
                    cTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGc                                                                                                                                  <  7:325058/149‑1 (MQ=255)
                     tACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCg                                                                                                                                 >  5:87741/1‑149 (MQ=255)
                       cAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCg                                                                                                                                 <  4:265646/147‑1 (MQ=255)
                          gagaAGGCGATTGATCTATTTTCCTGAAACAAGGTGATTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACAAGCAACGCCAAtta                                                                                                                                        >  2:386267/1‑137 (MQ=255)
                          gagaAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAtta                                                                                                                                        <  1:386267/137‑1 (MQ=255)
                          gagaAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAAtt                                                                                                                            >  7:67795/1‑149 (MQ=255)
                          gagaAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAAtt                                                                                                                            >  4:38690/1‑149 (MQ=255)
                               ggCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGTAATTTACGt                                                                                                                       >  6:35914/1‑149 (MQ=255)
                                gCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTc                                                                                                                      <  3:38690/149‑1 (MQ=255)
                                  gATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtta                                                                                                                                    <  3:362948/133‑1 (MQ=255)
                                  gATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtta                                                                                                                                    >  4:362948/1‑133 (MQ=255)
                                  gATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGa                                                                                                                    >  2:241996/1‑149 (MQ=255)
                                   aTTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGAt                                                                                                                   <  1:241996/149‑1 (MQ=255)
                                    ttGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATa                                                                                                                  <  3:214363/149‑1 (MQ=255)
                                       aTCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTc                                                                                                               >  2:60247/1‑149 (MQ=255)
                                       aTCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTc                                                                                                               >  2:224740/1‑149 (MQ=255)
                                        tCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCT‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCg                                                                                                              >  4:71884/1‑149 (MQ=255)
                                            ttttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTg                                                                                                          <  1:60247/149‑1 (MQ=255)
                                                  tGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCa                                                                                                    >  6:73458/1‑149 (MQ=255)
                                                  tGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCa                                                                                                    >  6:257209/1‑149 (MQ=255)
                                                        aaGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCa                                                                                               >  7:378520/1‑148 (MQ=255)
                                                           gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGc                                                                                           <  6:87741/149‑1 (MQ=255)
                                                                  ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCt                                                                                     >  3:91776/1‑148 (MQ=255)
                                                                  ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCACAAAGCGTTGCt                                                                                     >  5:278853/1‑148 (MQ=255)
                                                                   tCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCtt                                                                                    >  4:392463/1‑148 (MQ=255)
                                                                   tCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCtt                                                                                    >  5:158121/1‑148 (MQ=255)
                                                                    cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTa                                                                                   >  5:312685/1‑148 (MQ=255)
                                                                    cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTa                                                                                   >  5:312682/1‑148 (MQ=255)
                                                                     aaaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaaa                                                                                 >  2:411359/1‑149 (MQ=255)
                                                                      aaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtt                                                                                                                                     >  6:290198/1‑96 (MQ=255)
                                                                      aaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtt                                                                                                                                     <  5:290198/96‑1 (MQ=255)
                                                                      aaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaaa                                                                                 >  6:144885/1‑148 (MQ=255)
                                                                        aaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAc                                                                              <  1:224740/149‑1 (MQ=255)
                                                                        aaCTCCTGTCAAATTGCCTTTCGCCCTGAAAAATGCATACTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAc                                                                              <  3:189510/149‑1 (MQ=255)
                                                                          cTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCg                                                                            <  3:392463/149‑1 (MQ=255)
                                                                           tCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGt                                                                           >  5:183757/1‑149 (MQ=255)
                                                                             cTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                                                                          >  3:116393/1‑148 (MQ=255)
                                                                              tGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTg                                                                         >  5:199041/1‑148 (MQ=255)
                                                                                tCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaa                                                                      <  6:312685/149‑1 (MQ=255)
                                                                                tCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaa                                                                      <  6:312682/149‑1 (MQ=255)
                                                                                tCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaa                                                                       >  7:303101/1‑148 (MQ=255)
                                                                                tCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaa                                                                       >  5:250822/1‑148 (MQ=255)
                                                                                tCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGAGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaa                                                                       >  5:369472/1‑148 (MQ=255)
                                                                                  aaaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa                                                                    <  5:144885/149‑1 (MQ=255)
                                                                                     ttGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa                                                                 <  5:73458/149‑1 (MQ=255)
                                                                                        ccTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAct                                                               <  6:369472/148‑1 (MQ=255)
                                                                                        ccTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAct                                                               <  6:250822/148‑1 (MQ=255)
                                                                                        ccTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAACTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActc                                                              >  5:389000/1‑149 (MQ=255)
                                                                                          ttttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTg                                                            >  5:337498/1‑149 (MQ=255)
                                                                                           tttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg                                                           >  3:104037/1‑149 (MQ=255)
                                                                                              gCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt                                                          >  7:353724/1‑147 (MQ=255)
                                                                                               cccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCa                                                        <  8:353724/148‑1 (MQ=255)
                                                                                               cccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAg                                                       <  6:183757/149‑1 (MQ=255)
                                                                                                  tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaa                                                     >  5:376346/1‑148 (MQ=255)
                                                                                                       aaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAACTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATccc                                               <  6:389000/149‑1 (MQ=255)
                                                                                                          gCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAATTCCCCaa                                            >  6:260997/1‑149 (MQ=255)
                                                                                                           cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCaa                                            <  5:35914/148‑1 (MQ=255)
                                                                                                           cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCaa                                            <  4:91776/148‑1 (MQ=255)
                                                                                                           cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAAt                                           <  4:388367/149‑1 (MQ=255)
                                                                                                           cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAAt                                           <  2:205470/149‑1 (MQ=255)
                                                                                                            aTAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATg                                          <  8:303101/149‑1 (MQ=255)
                                                                                                              aGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGa                                         >  6:174602/1‑148 (MQ=255)
                                                                                                                   aaTATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCg                                   <  5:257209/149‑1 (MQ=255)
                                                                                                                     tatCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGAt                                 <  8:67795/149‑1 (MQ=255)
                                                                                                                     tatCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGAt                                 <  4:379854/149‑1 (MQ=255)
                                                                                                                       tCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaaa                                                    >  8:163938/1‑128 (MQ=255)
                                                                                                                       tCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGAtgc                               >  2:437356/1‑149 (MQ=255)
                                                                                                                       tCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAAAGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaaa                                                    <  7:163938/128‑1 (MQ=255)
                                                                                                                                aCTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGtt                      >  3:327258/1‑149 (MQ=255)
                                                                                                                                  tAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGtttt                    >  7:352466/1‑149 (MQ=255)
                                                                                                                                    aTTGGCGTTGCATGTAC‑TGCAAAGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAAGCCCCACTGCAAGCCGCTGCCGGCCAGttttt                   >  6:309737/1‑148 (MQ=255)
                                                                                                                                    aTTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTtc                  <  6:199041/149‑1 (MQ=255)
                                                                                                                                          gTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCatcatc            >  7:291708/1‑149 (MQ=255)
                                                                                                                                             gCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGAtg                                <  7:382246/126‑1 (MQ=255)
                                                                                                                                             gCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGAtg                                >  8:382246/1‑126 (MQ=255)
                                                                                                                                                  tACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCTGATGCCTGCCAGTTTTTCATCATCTGCGGCa     >  7:173736/1‑149 (MQ=255)
                                                                                                                                                     aTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATgg  <  6:337498/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCGCAATATCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGG  >  NC_000913/1446214‑1446510

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