Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,004,160 (G)7→8 coding (214/711 nt) ycbF → putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,004,1531.G93.8% 321.4 / 14.3 96coding (207/711 nt)ycbFputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base . (5/1);  new base G (56/34);  total (61/35)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 106 total aligned reads displayed.

CTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑GG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAA  >  NC_000913/1004054‑1004250
                                                                                                    |                                                                                                   
cTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGcc                                                                                                      >  1:981644/1‑100 (MQ=255)
 tttCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑g                                                                                                    >  1:1379483/1‑100 (MQ=255)
 tttCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGAAAAATCAACCGCCCAACGCGTTGGGCC‑g                                                                                                    >  2:726342/1‑100 (MQ=255)
   tCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑gg‑g                                                                                                 <  1:1284355/100‑1 (MQ=255)
     cGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑gg‑ggg                                                                                               <  1:101238/100‑1 (MQ=255)
      ggCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑gg‑gggg                                                                                              >  1:1205982/1‑100 (MQ=255)
       gCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑gg‑gggg‑g                                                                                            <  1:323108/100‑1 (MQ=255)
          gacagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑gg‑gggg‑ggaa                                                                                         >  1:321258/1‑97 (MQ=255)
          gacagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑gg‑gggg‑ggaa                                                                                         >  1:2626929/1‑97 (MQ=255)
             agagAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCg                                                                                       <  2:2452387/100‑1 (MQ=255)
             agagAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCg                                                                                       <  1:2615518/100‑1 (MQ=255)
              gagaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGt                                                                                      >  1:932432/1‑100 (MQ=255)
              gagaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGt                                                                                      >  1:193690/1‑100 (MQ=255)
                gaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCg                                                                                    <  1:2652694/100‑1 (MQ=255)
                  aTCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGAc                                                                                  <  1:812294/100‑1 (MQ=255)
                    cGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAg                                                                                <  1:2416408/100‑1 (MQ=255)
                    cGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTCATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAg                                                                                <  2:1847816/100‑1 (MQ=255)
                      atatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGt                                                                               <  2:347993/99‑1 (MQ=255)
                         tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCa                                                                           >  1:1250837/1‑100 (MQ=255)
                          cTATTTTACTGGCGAATACCAGTCAGGAATCCTGGCTTATTAACGGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAg                                                                          <  2:561266/100‑1 (MQ=255)
                          cTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAg                                                                          >  1:527094/1‑100 (MQ=255)
                            aTTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCa                                                                        >  1:49515/1‑100 (MQ=255)
                            aTTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCa                                                                        >  2:2834004/1‑100 (MQ=255)
                               ttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCa                                                                     <  2:2424745/100‑1 (MQ=255)
                                  cTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTa                                                                  <  2:2713257/100‑1 (MQ=255)
                                  cTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTa                                                                  >  1:877631/1‑100 (MQ=255)
                                  cTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTa                                                                  >  1:2733494/1‑100 (MQ=255)
                                  cTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTa                                                                  >  1:440268/1‑100 (MQ=255)
                                    gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACt                                                                >  2:1390231/1‑100 (MQ=255)
                                    gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACt                                                                >  1:2013868/1‑100 (MQ=255)
                                    gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACt                                                                >  2:2533151/1‑100 (MQ=255)
                                    gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACt                                                                >  2:1041998/1‑100 (MQ=255)
                                    gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGCAAAATCAACCGCCCAACGCGTTGCGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACt                                                                >  2:1042327/1‑100 (MQ=255)
                                      cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTgg                                                              <  2:2782104/100‑1 (MQ=255)
                                        aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGcc                                                            >  1:1830311/1‑100 (MQ=255)
                                        aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGcc                                                            <  2:876442/100‑1 (MQ=255)
                                        aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGcc                                                            >  1:1233350/1‑100 (MQ=255)
                                          tACCAGTCAGGGATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTCGGGCC‑GGTGGGG‑GGAGCGTCGACAGTGCCAGCACCATTACTGGCCGc                                                          >  2:925956/1‑100 (MQ=255)
                                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCt                                                         >  2:1133228/1‑100 (MQ=255)
                                             cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑GG‑GGGGTGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                                       >  1:2171776/1‑100 (MQ=255)
                                              aGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg                                                      >  1:2260634/1‑100 (MQ=255)
                                              aGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg                                                      >  2:922412/1‑100 (MQ=255)
                                                tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc                                                    >  1:1803023/1‑100 (MQ=255)
                                                tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc                                                    <  2:508294/100‑1 (MQ=255)
                                                 cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa                                                   >  1:345539/1‑100 (MQ=255)
                                                 cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa                                                   <  2:1205978/100‑1 (MQ=255)
                                                   ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCGAGCACCATTACTGGCCGCTCCGCCACt                                                 >  2:544762/1‑100 (MQ=255)
                                                   ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACt                                                 >  2:1409912/1‑100 (MQ=255)
                                                   ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACt                                                 >  1:946455/1‑100 (MQ=255)
                                                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAt                                              <  2:444117/100‑1 (MQ=255)
                                                         cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtct                                           <  2:255037/100‑1 (MQ=255)
                                                         cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtct                                           <  2:255911/100‑1 (MQ=255)
                                                           ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTcc                                         >  2:1731751/1‑100 (MQ=255)
                                                            gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCt                                        <  2:1233346/100‑1 (MQ=255)
                                                             cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                                       <  1:1112958/100‑1 (MQ=255)
                                                             cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                                       >  2:858614/1‑100 (MQ=255)
                                                             cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                                       >  2:2435826/1‑100 (MQ=255)
                                                             cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                                       >  1:819427/1‑100 (MQ=255)
                                                               tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGaa                                     >  1:2447744/1‑100 (MQ=255)
                                                                 ttaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGc                                   >  2:753409/1‑100 (MQ=255)
                                                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                                  >  2:1535654/1‑100 (MQ=255)
                                                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                                  <  1:1390236/100‑1 (MQ=255)
                                                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑GG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGccc                                 >  1:1277022/1‑100 (MQ=255)
                                                                   aaCAGTAAAATCAACCGCCCAACGCGTTGGGCC‑GG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcccn                                <  2:2200958/100‑2 (MQ=255)
                                                                     cAGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCgg                               >  1:2232909/1‑100 (MQ=255)
                                                                      aGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGt                              <  1:2207771/100‑1 (MQ=255)
                                                                      aGTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGt                              >  1:2882066/1‑100 (MQ=255)
                                                                       gTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                             <  1:698807/100‑1 (MQ=255)
                                                                       gTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                             >  1:1703149/1‑100 (MQ=255)
                                                                       gTAAAATCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                             <  2:506223/100‑1 (MQ=255)
                                                                         aaaaTCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACg                           <  1:1409917/100‑1 (MQ=255)
                                                                         aaaaTCAACCCCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACg                           <  1:225010/100‑1 (MQ=255)
                                                                          aaaTCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGa                          >  1:747037/1‑100 (MQ=255)
                                                                          aaaTCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGa                          >  1:1157114/1‑100 (MQ=255)
                                                                          aaaTCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGa                          >  1:1674988/1‑100 (MQ=255)
                                                                            aTCAACCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACa                        <  2:14435/100‑1 (MQ=255)
                                                                             cccccggCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAg                       >  1:2536107/7‑100 (MQ=255)
                                                                              ccccggCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                      >  2:2441680/6‑100 (MQ=255)
                                                                              ccccggCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTAGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                      >  2:52384/6‑100 (MQ=255)
                                                                               aaCCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc                     <  2:1803010/100‑1 (MQ=255)
                                                                               aaCCGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc                     >  2:591324/1‑100 (MQ=255)
                                                                                 nggCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCAc                   >  1:854888/3‑100 (MQ=255)
                                                                                 cggCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTCTTCTCCTGAAGCCCGGTACGACAGGCAc                   >  1:687309/3‑100 (MQ=255)
                                                                                  cGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTTCGCCACTTATTCTCCTGACGCCCGGTACGACAGTCACg                  <  2:628662/100‑1 (MQ=255)
                                                                                  cGCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACg                  >  2:2514399/1‑100 (MQ=255)
                                                                                   gCCCAACTCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGt                 >  1:2707598/1‑100 (MQ=255)
                                                                                   gCCCAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGt                 >  1:2519240/1‑100 (MQ=255)
                                                                                     ccAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTg               >  2:2735452/1‑100 (MQ=255)
                                                                                     ccAACGCGTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTg               <  1:364255/100‑1 (MQ=255)
                                                                                         cgcgTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCt           >  2:1051759/1‑100 (MQ=255)
                                                                                         cgcgTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCt           <  2:240349/100‑1 (MQ=255)
                                                                                           cgTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTg         >  2:827773/1‑100 (MQ=255)
                                                                                           cgTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTg         >  2:828009/1‑100 (MQ=255)
                                                                                            gTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGc        >  1:2424889/1‑100 (MQ=255)
                                                                                            gTTGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGc        >  1:1402463/1‑100 (MQ=255)
                                                                                             ttGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCt       <  2:2863051/100‑1 (MQ=255)
                                                                                              tGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTg      <  1:2148071/100‑1 (MQ=255)
                                                                                              tGGGCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTg      <  2:1393796/100‑1 (MQ=255)
                                                                                               gggCCGGG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTga     >  2:9601/1‑100 (MQ=255)
                                                                                                 gCC‑GG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGaa  >  1:1649617/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
CTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCC‑GG‑GGGG‑GAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAA  >  NC_000913/1004054‑1004250

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