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.

CTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAA  >  NC_000913/1004054‑1004250
                                                                                                    |                                                                                                   
cTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGcc                                                                                                      >  1:981644/1‑100 (MQ=255)
 tttCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCg                                                                                                    >  1:1379483/1‑100 (MQ=255)
 tttCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGAAAAATCAACCGCCCAACGCGTTGGGCCg                                                                                                    >  2:726342/1‑100 (MQ=255)
   tCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggg                                                                                                 <  1:1284355/100‑1 (MQ=255)
     cGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggggg                                                                                               <  1:101238/100‑1 (MQ=255)
      ggCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgggggg                                                                                              >  1:1205982/1‑100 (MQ=255)
       gCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggggggg                                                                                            <  1:323108/100‑1 (MQ=255)
          gacagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggggggggaa                                                                                         >  1:321258/1‑97 (MQ=255)
          gacagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggggggggaa                                                                                         >  1:2626929/1‑97 (MQ=255)
             agagAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCg                                                                                       <  2:2452387/100‑1 (MQ=255)
             agagAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCg                                                                                       <  1:2615518/100‑1 (MQ=255)
              gagaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGt                                                                                      >  1:932432/1‑100 (MQ=255)
              gagaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGt                                                                                      >  1:193690/1‑100 (MQ=255)
                gaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCg                                                                                    <  1:2652694/100‑1 (MQ=255)
                  aTCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGAc                                                                                  <  1:812294/100‑1 (MQ=255)
                    cGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAg                                                                                <  1:2416408/100‑1 (MQ=255)
                    cGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTCATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAg                                                                                <  2:1847816/100‑1 (MQ=255)
                      atatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGt                                                                               <  2:347993/99‑1 (MQ=255)
                         tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCa                                                                           >  1:1250837/1‑100 (MQ=255)
                          cTATTTTACTGGCGAATACCAGTCAGGAATCCTGGCTTATTAACGGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAg                                                                          <  2:561266/100‑1 (MQ=255)
                          cTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAg                                                                          >  1:527094/1‑100 (MQ=255)
                            aTTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCa                                                                        >  1:49515/1‑100 (MQ=255)
                            aTTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCa                                                                        >  2:2834004/1‑100 (MQ=255)
                               ttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCa                                                                     <  2:2424745/100‑1 (MQ=255)
                                  cTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTa                                                                  <  2:2713257/100‑1 (MQ=255)
                                  cTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTa                                                                  >  1:877631/1‑100 (MQ=255)
                                  cTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTa                                                                  >  1:2733494/1‑100 (MQ=255)
                                  cTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTa                                                                  >  1:440268/1‑100 (MQ=255)
                                    gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACt                                                                >  2:1390231/1‑100 (MQ=255)
                                    gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACt                                                                >  1:2013868/1‑100 (MQ=255)
                                    gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACt                                                                >  2:2533151/1‑100 (MQ=255)
                                    gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACt                                                                >  2:1041998/1‑100 (MQ=255)
                                    gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGCAAAATCAACCGCCCAACGCGTTGCGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACt                                                                >  2:1042327/1‑100 (MQ=255)
                                      cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTgg                                                              <  2:2782104/100‑1 (MQ=255)
                                        aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                                                            >  1:1830311/1‑100 (MQ=255)
                                        aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                                                            <  2:876442/100‑1 (MQ=255)
                                        aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                                                            >  1:1233350/1‑100 (MQ=255)
                                          tACCAGTCAGGGATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTCGGGCCGGTGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGc                                                          >  2:925956/1‑100 (MQ=255)
                                           aCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCt                                                         >  2:1133228/1‑100 (MQ=255)
                                             cAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGTGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                                       >  1:2171776/1‑100 (MQ=255)
                                              aGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg                                                      >  1:2260634/1‑100 (MQ=255)
                                              aGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg                                                      >  2:922412/1‑100 (MQ=255)
                                                tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc                                                    >  1:1803023/1‑100 (MQ=255)
                                                tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc                                                    <  2:508294/100‑1 (MQ=255)
                                                 cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa                                                   >  1:345539/1‑100 (MQ=255)
                                                 cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa                                                   <  2:1205978/100‑1 (MQ=255)
                                                   ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCGAGCACCATTACTGGCCGCTCCGCCACt                                                 >  2:544762/1‑100 (MQ=255)
                                                   ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACt                                                 >  2:1409912/1‑100 (MQ=255)
                                                   ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACt                                                 >  1:946455/1‑100 (MQ=255)
                                                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAt                                              <  2:444117/100‑1 (MQ=255)
                                                         cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtct                                           <  2:255037/100‑1 (MQ=255)
                                                         cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtct                                           <  2:255911/100‑1 (MQ=255)
                                                           ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTcc                                         >  2:1731751/1‑100 (MQ=255)
                                                            gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCt                                        <  2:1233346/100‑1 (MQ=255)
                                                             cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                                       <  1:1112958/100‑1 (MQ=255)
                                                             cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                                       >  2:858614/1‑100 (MQ=255)
                                                             cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                                       >  2:2435826/1‑100 (MQ=255)
                                                             cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                                       >  1:819427/1‑100 (MQ=255)
                                                               tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGaa                                     >  1:2447744/1‑100 (MQ=255)
                                                                 ttaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGc                                   >  2:753409/1‑100 (MQ=255)
                                                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                                  >  2:1535654/1‑100 (MQ=255)
                                                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                                  <  1:1390236/100‑1 (MQ=255)
                                                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGccc                                 >  1:1277022/1‑100 (MQ=255)
                                                                   aaCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcccn                                <  2:2200958/100‑2 (MQ=255)
                                                                     cAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCgg                               >  1:2232909/1‑100 (MQ=255)
                                                                      aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGt                              <  1:2207771/100‑1 (MQ=255)
                                                                      aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGt                              >  1:2882066/1‑100 (MQ=255)
                                                                       gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                             <  1:698807/100‑1 (MQ=255)
                                                                       gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                             >  1:1703149/1‑100 (MQ=255)
                                                                       gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                             <  2:506223/100‑1 (MQ=255)
                                                                         aaaaTCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACg                           <  1:1409917/100‑1 (MQ=255)
                                                                         aaaaTCAACCCCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACg                           <  1:225010/100‑1 (MQ=255)
                                                                          aaaTCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGa                          >  1:747037/1‑100 (MQ=255)
                                                                          aaaTCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGa                          >  1:1157114/1‑100 (MQ=255)
                                                                          aaaTCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGa                          >  1:1674988/1‑100 (MQ=255)
                                                                            aTCAACCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACa                        <  2:14435/100‑1 (MQ=255)
                                                                             cccccggCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAg                       >  1:2536107/7‑100 (MQ=255)
                                                                              ccccggCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                      >  2:2441680/6‑100 (MQ=255)
                                                                              ccccggCCCAACGCGTTGGGCCGGGGGGGGAAGCGTAGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                      >  2:52384/6‑100 (MQ=255)
                                                                               aaCCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc                     <  2:1803010/100‑1 (MQ=255)
                                                                               aaCCGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc                     >  2:591324/1‑100 (MQ=255)
                                                                                 nggCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCAc                   >  1:854888/3‑100 (MQ=255)
                                                                                 cggCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTCTTCTCCTGAAGCCCGGTACGACAGGCAc                   >  1:687309/3‑100 (MQ=255)
                                                                                  cGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTTCGCCACTTATTCTCCTGACGCCCGGTACGACAGTCACg                  <  2:628662/100‑1 (MQ=255)
                                                                                  cGCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACg                  >  2:2514399/1‑100 (MQ=255)
                                                                                   gCCCAACTCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGt                 >  1:2707598/1‑100 (MQ=255)
                                                                                   gCCCAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGt                 >  1:2519240/1‑100 (MQ=255)
                                                                                     ccAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTg               >  2:2735452/1‑100 (MQ=255)
                                                                                     ccAACGCGTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTg               <  1:364255/100‑1 (MQ=255)
                                                                                         cgcgTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCt           >  2:1051759/1‑100 (MQ=255)
                                                                                         cgcgTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCt           <  2:240349/100‑1 (MQ=255)
                                                                                           cgTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTg         >  2:827773/1‑100 (MQ=255)
                                                                                           cgTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTg         >  2:828009/1‑100 (MQ=255)
                                                                                            gTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGc        >  1:2424889/1‑100 (MQ=255)
                                                                                            gTTGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGc        >  1:1402463/1‑100 (MQ=255)
                                                                                             ttGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCt       <  2:2863051/100‑1 (MQ=255)
                                                                                              tGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTg      <  1:2148071/100‑1 (MQ=255)
                                                                                              tGGGCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTg      <  2:1393796/100‑1 (MQ=255)
                                                                                               gggCCGGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTga     >  2:9601/1‑100 (MQ=255)
                                                                                                 gCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGaa  >  1:1649617/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
CTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAA  >  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: