Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 103,614 A→T 5.2% K154* (AAG→TAG)  ftsQ → divisome assembly protein, membrane anchored protein involved in growth of wall at septum

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913103,6140AT 5.2% 323.0 / 3.3 105K154* (AAG→TAG) ftsQdivisome assembly protein, membrane anchored protein involved in growth of wall at septum
Reads supporting (aligned to +/- strand):  ref base A (51/48);  new base T (4/2);  total (55/50)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

CGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGTTCCTGGCAGTTGACGCT  >  NC_000913/103465‑103762
                                                                                                                                                     |                                                                                                                                                    
cGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCagca                                                                                                                                                      <  1:1708289/150‑1 (MQ=255)
      ttAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAgg                                                                                                                                                <  2:1769743/150‑1 (MQ=255)
        aaGCAGGTGAGCGTCAGAAAGCAGTTGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTg                                                                                                                                              >  2:85948/1‑150 (MQ=255)
          gCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCt                                                                                                                                            <  1:784411/150‑1 (MQ=255)
            aGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTc                                                                                                                                          <  2:1237593/150‑1 (MQ=255)
             ggTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGt                                                                                                                                               >  1:503565/1‑144 (MQ=255)
             ggTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGt                                                                                                                                               <  2:503565/144‑1 (MQ=255)
                gAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCTGCTTGCTGGTGcgttctgg                                                                                                                                      <  1:1497490/150‑8 (MQ=255)
                 cgCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATg                                                                                                                                     <  2:138802/149‑1 (MQ=255)
                 aGCGTCAGAAATCAGTGGGCTTATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGGGGAATGATCAACATATGGTAGACTCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATg                                                                                                                                     <  2:450800/150‑1 (MQ=255)
                  gCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGATTATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGc                                                                                                                                    <  2:969484/150‑1 (MQ=255)
                  gCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCATCATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGc                                                                                                                                    <  1:78099/150‑1 (MQ=255)
                  gCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGc                                                                                                                                    <  2:1010943/150‑1 (MQ=255)
                  gCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGc                                                                                                                                    <  2:478416/150‑1 (MQ=255)
                  gCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGc                                                                                                                                    >  1:316646/1‑150 (MQ=255)
                       aGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTAt                                                                                                                               >  2:1372331/1‑150 (MQ=255)
                        gAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATg                                                                                                                              >  2:555705/1‑150 (MQ=255)
                           aGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCTCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGcc                                                                                                                           <  1:522201/150‑1 (MQ=255)
                            gCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGccc                                                                                                                          >  1:1328260/1‑150 (MQ=255)
                             cAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCg                                                                                                                         <  1:655507/150‑1 (MQ=255)
                              aGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTTCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCgg                                                                                                                        <  1:1729506/150‑1 (MQ=255)
                                 ggCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAg                                                                                                                     <  2:909889/150‑1 (MQ=255)
                                   ccTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAATGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATTCTGTATGGCCCGGAAGGc                                                                                                                   >  1:128980/1‑150 (MQ=255)
                                   ccTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGc                                                                                                                   >  1:1173937/1‑150 (MQ=255)
                                       atgaATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAgcg                                                                                                               >  1:196163/1‑150 (MQ=255)
                                       atgaATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAgcg                                                                                                               <  1:1383170/150‑1 (MQ=255)
                                        tgaATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAgcgc                                                                                                              >  1:180007/1‑150 (MQ=255)
                                        tgaATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAgcgc                                                                                                              <  2:338362/150‑1 (MQ=255)
                                           aTTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCaa                                                                                                           >  2:1073988/1‑150 (MQ=255)
                                              gAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGa                                                                                                        <  2:1561535/150‑1 (MQ=255)
                                                 gATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgt                                                                                                     >  2:1181822/1‑150 (MQ=255)
                                                  aTTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAgtg                                                                                                    <  2:1368891/150‑1 (MQ=255)
                                                       tCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCa                                                                                               <  2:642243/150‑1 (MQ=255)
                                                          ggTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCTAGCAGGTGATTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAggg                                                                                            >  2:1319385/1‑150 (MQ=255)
                                                          ggTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAggg                                                                                            <  2:1328260/150‑1 (MQ=255)
                                                           gTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGc                                                                                           >  1:1699962/1‑150 (MQ=255)
                                                             tGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGGAGGCAGCGCCAATGAAGTGTTGCAGGGCTa                                                                                         >  2:1087429/1‑150 (MQ=255)
                                                              gAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTAt                                                                                        >  2:1011034/1‑150 (MQ=255)
                                                                  atGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATcgcg                                                                                    >  1:1280114/1‑150 (MQ=255)
                                                                   tgtgCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGa                                                                                   <  1:492524/150‑1 (MQ=255)
                                                                    gtgCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGaa                                                                                  >  2:1484130/1‑150 (MQ=255)
                                                                     tgCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGaaa                                                                                 <  1:1657853/150‑1 (MQ=255)
                                                                        cGATTGCGCTGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATgg                                                                              <  1:105914/150‑1 (MQ=255)
                                                                           ttGCGCGGTGGAATGAGCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCTGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGc                                                                           <  1:510809/150‑1 (MQ=255)
                                                                                   tGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTg                                                                   >  2:1456413/1‑150 (MQ=255)
                                                                                   tGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTg                                                                   >  2:647287/1‑150 (MQ=255)
                                                                                    ggAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTgg                                                                  <  1:8126/150‑1 (MQ=255)
                                                                                     gAATGATCAACATATTGTAGCCGCGGAATGAAATACCTTCAGCGTGCCGCCAGAACTCACCAGCAAGCAGGTGCGTCCAATGCTGGATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGc                                                                 <  1:768100/150‑1 (MQ=255)
                                                                                     gAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGc                                                                 <  2:1206834/150‑1 (MQ=255)
                                                                                     gAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGc                                                                 <  1:776896/150‑1 (MQ=255)
                                                                                      aaTGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCa                                                                <  2:1021493/150‑1 (MQ=255)
                                                                                        tGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCaaa                                                              >  1:498812/1‑150 (MQ=255)
                                                                                        tGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCaaa                                                              >  1:1575865/1‑150 (MQ=255)
                                                                                            cAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGAc                                                          >  1:1374025/1‑150 (MQ=255)
                                                                                              aCATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAg                                                        <  2:21662/150‑1 (MQ=255)
                                                                                              aCATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAg                                                        <  2:1206583/150‑1 (MQ=255)
                                                                                               cATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGa                                                       <  1:734296/150‑1 (MQ=255)
                                                                                                 tatGGTAGACGCGGAAGGAAATTCCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCATAGGACAGAtt                                                     >  1:1740693/1‑150 (MQ=255)
                                                                                                 tatGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGAtt                                                     >  1:1751567/1‑150 (MQ=255)
                                                                                                 tatGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGAtt                                                     >  1:1733708/1‑150 (MQ=255)
                                                                                                      tAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTActc                                                >  1:405752/1‑150 (MQ=255)
                                                                                                       tgACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTActct                                               <  2:1121954/149‑1 (MQ=255)
                                                                                                       aGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTActct                                               >  2:1543818/1‑150 (MQ=255)
                                                                                                        gACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCGTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCCGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTg                                              <  1:924563/150‑1 (MQ=255)
                                                                                                           gcgGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTgaag                                           <  1:1372331/150‑1 (MQ=255)
                                                                                                             ggaaggaaATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGGTTCGAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTgaagga                                         <  2:631475/150‑1 (MQ=255)
                                                                                                                aggaaATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGACgc                                      >  1:925543/1‑150 (MQ=255)
                                                                                                                aggaaATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAgc                                      >  1:1395908/1‑150 (MQ=255)
                                                                                                                  gaaATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAgcgg                                    >  2:951596/1‑150 (MQ=255)
                                                                                                                      tACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGAt                                >  2:1389332/1‑150 (MQ=255)
                                                                                                                       aCCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATg                               >  2:942488/1‑150 (MQ=255)
                                                                                                                       aCCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATg                               >  2:409386/1‑150 (MQ=255)
                                                                                                                       aCCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTGCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGTCTATCGCGAAATGGGGCAGATGTTGGCAAATGACAGATTTACTCTGAAGGAAGCGGCGATg                               <  2:1646341/150‑1 (MQ=255)
                                                                                                                         cTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGAc                             >  1:1805411/1‑150 (MQ=255)
                                                                                                                         cTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGAGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGAc                             >  2:127012/1‑150 (MQ=255)
                                                                                                                                 ggccgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCgcggc                     <  1:1069586/149‑1 (MQ=255)
                                                                                                                                  gccgccTCAACGCACCAGCAAGCAGGTGCTTCCAATGCGGTATGGCCCGGAAGGCAGCGCCAATGAAGTGGTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCgcggcg                    <  1:450763/150‑1 (MQ=255)
                                                                                                                                  gccgccAGAACGCACCAGCAAGCAGCTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCgcggcg                    <  1:991576/150‑1 (MQ=255)
                                                                                                                                   ccgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGTGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAATCGGCGATGACCGCGCGGCGt                   >  2:1687686/1‑150 (MQ=255)
                                                                                                                                   ccgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGt                   >  1:68733/1‑150 (MQ=255)
                                                                                                                                    cgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGtt                  >  2:1306593/1‑150 (MQ=255)
                                                                                                                                    cgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGtt                  >  1:740284/1‑150 (MQ=255)
                                                                                                                                    cgccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGtt                  >  1:1544543/1‑150 (MQ=255)
                                                                                                                                     gccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGTTc                 >  1:1410084/1‑150 (MQ=255)
                                                                                                                                     gccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGTTc                 >  1:1408446/1‑150 (MQ=255)
                                                                                                                                     gccAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGTTc                 >  2:602850/1‑150 (MQ=255)
                                                                                                                                      ccAGAACGCACCATCAAGCAGGTGCTTCCAATGCTTTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGATCGCGCGGCGTTcc                <  2:340974/150‑1 (MQ=255)
                                                                                                                                       ttggaagCACCTGCTTGCTGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGTTCCt               >  1:1639194/7‑150 (MQ=255)
                                                                                                                                         ggaagCACCTGCTTGCTGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATTAAGTGTTCCAGGGCTATCGCGAAATTCGGCAGATGCTGGCAAAGGACAGATTTACTAGGAAGGAAGCGGAGATGACCGCGCGGCGTTCCTgg             >  2:1102668/5‑150 (MQ=255)
                                                                                                                                         ggaagCACCTGCTTGCTGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGTTCCTgg             >  1:24311/5‑150 (MQ=255)
                                                                                                                                         gAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGTTCCTgg             <  1:1073988/150‑1 (MQ=255)
                                                                                                                                          aaCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGTTCCTGGc            <  2:1173937/150‑1 (MQ=255)
                                                                                                                                          aaCGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGTTCCTGGc            <  1:1681018/150‑1 (MQ=255)
                                                                                                                                            cGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGTTCCTGGCAg          >  1:1500775/1‑150 (MQ=255)
                                                                                                                                              cACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGTTCCTGGCAGtt        >  1:842885/1‑150 (MQ=255)
                                                                                                                                                ccAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGTTCCTGGCAGTTgc      >  2:1275716/1‑149 (MQ=255)
                                                                                                                                                ccAGCAAGCAGGTGCTTCCAATGATGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGAGGCGTTCCTGGCAGTTGa      >  1:550402/1‑150 (MQ=255)
                                                                                                                                                  agcaagcaGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGATCCTGGCAGTTGACg    >  2:859147/1‑150 (MQ=255)
                                                                                                                                                   gcaagcaGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGTTCCTGGCAGTTGACGc   >  2:528354/1‑150 (MQ=255)
                                                                                                                                                    caagcaGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGTTCCTGGCAGTTGACGCt  >  1:1099423/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGTGGATTAAGCAGGTGAGCGTCAGAAAGCAGTGGCCTGATGAATTGAAGATTCATCTGGTTGAATATGTGCCGATTGCGCGGTGGAATGATCAACATATGGTAGACGCGGAAGGAAATACCTTCAGCGTGCCGCCAGAACGCACCAGCAAGCAGGTGCTTCCAATGCTGTATGGCCCGGAAGGCAGCGCCAATGAAGTGTTGCAGGGCTATCGCGAAATGGGGCAGATGCTGGCAAAGGACAGATTTACTCTGAAGGAAGCGGCGATGACCGCGCGGCGTTCCTGGCAGTTGACGCT  >  NC_000913/103465‑103762

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: