Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,139,669 +C coding (233/330 nt)
coding (132/159 nt)
yijI →
b4558 ←
hypothetical protein
hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,139,6671.C99.6% 804.0 / ‑3.8 254coding (231/330 nt)
coding (134/159 nt)
yijI
b4558
hypothetical protein
hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (120/133);  total (121/133)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 289 total aligned reads displayed.

GCTTCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACGCGCAT  >  NC_000913/4139575‑4139754
                                                                                             |                                                                                       
 cTTCCTGTGCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGc                                                                                        <  2:3944749/93‑1 (MQ=255)
  ttCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGcc                                                                                       >  1:181578/1‑93 (MQ=255)
    ccTGTTCGTCGGCCCCGGAAATCTGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccct                                                                                     <  2:3249047/93‑3 (MQ=255)
    ccTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccct                                                                                     <  2:3965699/93‑3 (MQ=255)
     cTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccctt                                                                                    <  1:1297981/93‑4 (MQ=255)
     cTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccctt                                                                                    <  2:4603655/93‑4 (MQ=255)
     cTGTACGTCGGCACCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccctt                                                                                    <  2:2748429/93‑4 (MQ=255)
      tGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGcccttg                                                                                   >  2:4335977/1‑89 (MQ=255)
      tGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGcccttg                                                                                   >  2:1411473/1‑89 (MQ=255)
      tGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGcccttg                                                                                   >  2:4444900/1‑89 (MQ=255)
           gtcgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTc                                                                               >  2:1019517/1‑93 (MQ=255)
            tcgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTTCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCa                                                                              <  1:3300747/93‑1 (MQ=255)
              ggCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGt                                                                            <  1:3132754/93‑1 (MQ=255)
                ccccGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTgc                                                                          <  1:862808/93‑1 (MQ=255)
                 cccGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTgcg                                                                         >  1:20291/1‑93 (MQ=255)
                 cccGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTgcg                                                                         >  2:742925/1‑93 (MQ=255)
                  ccGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTgcgc                                                                        >  2:2348722/1‑93 (MQ=255)
                    ggAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAg                                                                      <  2:2048402/93‑1 (MQ=255)
                    ggAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAg                                                                      <  2:2968303/93‑1 (MQ=255)
                     gAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAgg                                                                     >  1:4425821/1‑93 (MQ=255)
                      aaaTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGt                                                                    >  1:3211903/1‑93 (MQ=255)
                       aaTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGtt                                                                   >  2:119054/1‑93 (MQ=255)
                       aaTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGtt                                                                   >  1:3192008/1‑93 (MQ=255)
                       aaTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTCTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGtt                                                                   >  1:2470417/1‑93 (MQ=255)
                        aTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTa                                                                  >  2:1330403/1‑93 (MQ=255)
                         tCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTAt                                                                 >  2:4360252/1‑93 (MQ=255)
                         tCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTAt                                                                 <  1:1710307/93‑1 (MQ=255)
                          cAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATg                                                                <  1:2065592/93‑1 (MQ=255)
                           aGTAAGTGGAGGTGATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGc                                                               <  2:1847706/93‑1 (MQ=255)
                           aGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGc                                                               >  1:1184216/1‑93 (MQ=255)
                           aGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGc                                                               >  1:2628287/1‑93 (MQ=255)
                           aGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGc                                                               <  2:1458051/93‑1 (MQ=255)
                            gTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGcc                                                              <  1:1824176/93‑1 (MQ=255)
                              aaCTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCAc                                                            <  2:1070323/93‑1 (MQ=255)
                                cTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTc                                                          >  2:4177172/1‑93 (MQ=255)
                                 tGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCa                                                         >  2:112987/1‑93 (MQ=255)
                                  ggCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCaa                                                        <  2:363541/93‑1 (MQ=255)
                                  ggCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCaa                                                        >  2:3429477/1‑93 (MQ=255)
                                  ggCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGCGCGCAGGTTATGCCACTCaa                                                        >  1:3456764/1‑93 (MQ=255)
                                   gCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAAt                                                       >  2:3032868/1‑93 (MQ=255)
                                   gCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAAt                                                       >  2:3155934/1‑93 (MQ=255)
                                   gCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAAt                                                       <  1:3019941/93‑1 (MQ=255)
                                      gTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtg                                                    >  1:355737/1‑93 (MQ=255)
                                      gTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtg                                                    <  2:2288357/93‑1 (MQ=255)
                                      gTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtg                                                    <  1:2793490/93‑1 (MQ=255)
                                       ttATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtgt                                                   >  1:3366620/1‑93 (MQ=255)
                                       ttATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtgt                                                   >  2:1165638/1‑93 (MQ=255)
                                        tATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtgtg                                                  <  2:1652113/93‑1 (MQ=255)
                                         accgcccgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtga                                                 <  2:4092090/91‑1 (MQ=255)
                                          tCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtgat                                                <  2:2168901/93‑1 (MQ=255)
                                          tCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtgat                                                <  1:94140/93‑1 (MQ=255)
                                          tCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtgat                                                <  2:1240124/93‑1 (MQ=255)
                                            gcccgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtgatga                                              >  2:580255/1‑93 (MQ=255)
                                               cgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaaaa                                           <  2:832174/93‑1 (MQ=255)
                                                 ccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAg                                         <  1:1991939/93‑1 (MQ=255)
                                                   aGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTa                                       >  2:1259576/1‑93 (MQ=255)
                                                   aGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTa                                       <  1:2890954/93‑1 (MQ=255)
                                                    gCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAt                                      >  2:3783680/1‑93 (MQ=255)
                                                    gCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAt                                      >  1:4172466/1‑93 (MQ=255)
                                                    gCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAt                                      >  1:3057583/1‑93 (MQ=255)
                                                     cGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAtt                                     <  1:4531155/93‑1 (MQ=255)
                                                      gTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTa                                    >  2:1986590/1‑93 (MQ=255)
                                                      gTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTa                                    >  1:3161574/1‑93 (MQ=255)
                                                       tataGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTac                                   <  1:1386367/93‑1 (MQ=255)
                                                       tATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTac                                   >  1:796797/1‑93 (MQ=255)
                                                       tATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTac                                   <  1:3609900/93‑1 (MQ=255)
                                                        aTCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTaca                                  >  1:386881/1‑93 (MQ=255)
                                                         tCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTacac                                 >  1:1553987/1‑93 (MQ=255)
                                                         tCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTacac                                 >  2:80247/1‑93 (MQ=255)
                                                         tCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTacac                                 <  1:46161/93‑1 (MQ=255)
                                                          cGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTAcaca                                <  1:3501148/93‑1 (MQ=255)
                                                          cGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTAcaca                                >  2:3597785/1‑93 (MQ=255)
                                                          cGGTGCCAATCAACGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTAcaca                                <  1:3718970/93‑1 (MQ=255)
                                                            gTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCTCTCAATTTGTGATGAAAAAGTATTACACAGc                              <  1:2529600/93‑1 (MQ=255)
                                                            gTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGc                              >  1:1894639/1‑93 (MQ=255)
                                                              gCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGt                            <  2:3890092/93‑1 (MQ=255)
                                                               ccAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGtt                           <  1:4335219/93‑1 (MQ=255)
                                                                cAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGttt                          <  2:2669652/93‑1 (MQ=255)
                                                                 aaTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTc                         <  1:220136/93‑1 (MQ=255)
                                                                  aTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCa                        >  1:592520/1‑93 (MQ=255)
                                                                   tCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCaa                       <  2:4291731/93‑1 (MQ=255)
                                                                     aGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACg                     <  1:105648/93‑1 (MQ=255)
                                                                     aGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACg                     <  2:499914/93‑1 (MQ=255)
                                                                      gcgcAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGt                    <  2:3121346/93‑1 (MQ=255)
                                                                         cAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGc                 >  1:2121715/1‑93 (MQ=255)
                                                                          aaGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCt                <  1:1015603/93‑1 (MQ=255)
                                                                           aGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCtt               >  2:2915892/1‑93 (MQ=255)
                                                                            gggCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTg              >  1:2387310/1‑93 (MQ=255)
                                                                            gggCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTg              >  2:2508897/1‑93 (MQ=255)
                                                                            gggCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTg              >  2:2036551/1‑93 (MQ=255)
                                                                             ggCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGc             <  2:2703343/93‑1 (MQ=255)
                                                                              gcgcTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGcc            >  2:3299319/1‑93 (MQ=255)
                                                                              gcgcTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGcc            <  1:1483125/93‑1 (MQ=255)
                                                                                 cTTTTGNCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGa         <  1:4084869/93‑1 (MQ=255)
                                                                                 cTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGa         <  1:3592444/93‑1 (MQ=255)
                                                                                  ttttGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGAc        >  2:3458231/1‑93 (MQ=255)
                                                                                   tttGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGaag       >  2:1921920/1‑91 (MQ=255)
                                                                                     tGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGAcgcg     <  1:2747715/93‑1 (MQ=255)
                                                                                       gCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACGCGCa   <  1:1646921/93‑1 (MQ=255)
                                                                                        cGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACGCGCAt  >  2:4128050/1‑93 (MQ=255)
                                                                                             |                                                                                       
GCTTCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACGCGCAT  >  NC_000913/4139575‑4139754

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