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.5% 653.4 / ‑4.3 201coding (231/330 nt)
coding (134/159 nt)
yijI
b4558
hypothetical protein
hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (102/98);  total (102/99)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 220 total aligned reads displayed.

GCTTCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACGCGCAT  >  NC_000913/4139575‑4139754
                                                                                             |                                                                                        
gcTTCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTg                                                                                           >  1:4680583/1‑93 (MQ=255)
 cTTCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGc                                                                                         <  2:3757353/93‑1 (MQ=255)
   tCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccc                                                                                       >  2:4646561/1‑92 (MQ=255)
   tCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccc                                                                                       >  2:1372980/1‑92 (MQ=255)
     cTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAATGGCGCTTTTTGCGTTGccctt                                                                                     >  1:1176269/1‑90 (MQ=255)
     cTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccctt                                                                                     >  2:2692818/1‑90 (MQ=255)
         tcgtcgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCgg                                                                                 >  2:4731073/1‑93 (MQ=255)
         tcgtcgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCgg                                                                                 >  1:3028267/1‑93 (MQ=255)
          cgtcgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGt                                                                                <  2:2118864/93‑1 (MQ=255)
          cgtcgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGt                                                                                >  1:307531/1‑93 (MQ=255)
           gtcgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTc                                                                               >  2:3658714/1‑93 (MQ=255)
            tcgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCa                                                                              >  2:1240581/1‑93 (MQ=255)
               gCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTg                                                                           >  1:3926008/1‑93 (MQ=255)
                gctgggAAAGCAGAAACTGGCAGTTTACGCCCGCCAGTGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTAAGTgc                                                                          <  2:627491/89‑1 (MQ=255)
                 cccGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTgcg                                                                         >  2:1455512/1‑93 (MQ=255)
                  ccGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTgcgc                                                                        >  1:1330552/1‑93 (MQ=255)
                  ccGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTgcgc                                                                        >  2:3298696/1‑93 (MQ=255)
                   cGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCa                                                                       >  1:4201531/1‑93 (MQ=255)
                   cGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCa                                                                       <  1:717846/93‑1 (MQ=255)
                    ggAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAg                                                                      >  2:757091/1‑93 (MQ=255)
                     gAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAgg                                                                     >  1:1984133/1‑93 (MQ=255)
                      aaaTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGt                                                                    <  2:3946493/93‑1 (MQ=255)
                       aaTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGtt                                                                   >  1:2327104/1‑93 (MQ=255)
                        aTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTa                                                                  <  1:2356554/93‑1 (MQ=255)
                        aTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTa                                                                  <  1:1127091/93‑1 (MQ=255)
                        aTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTa                                                                  <  2:1619549/93‑1 (MQ=255)
                         tCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTAt                                                                 <  2:3102418/93‑1 (MQ=255)
                          cAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATg                                                                <  1:3591604/93‑1 (MQ=255)
                          cAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATg                                                                >  2:871860/1‑93 (MQ=255)
                            gTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGcc                                                              >  2:3096749/1‑93 (MQ=255)
                              aaCTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCAc                                                            <  1:1614383/93‑1 (MQ=255)
                                cTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTc                                                          <  1:2623644/93‑1 (MQ=255)
                                 tGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCa                                                         >  2:2455003/1‑93 (MQ=255)
                                  ggCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCaa                                                        <  1:2800041/93‑1 (MQ=255)
                                    cAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAAtt                                                      <  1:1859765/93‑1 (MQ=255)
                                    cAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAAtt                                                      >  2:2467951/1‑93 (MQ=255)
                                    cAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAAtt                                                      <  2:2669089/93‑1 (MQ=255)
                                      gTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtg                                                    <  2:2965938/93‑1 (MQ=255)
                                       ttATCGCCCGCCAGAGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtgt                                                   >  2:2358623/1‑93 (MQ=255)
                                        tATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtgtg                                                  >  2:2298171/1‑93 (MQ=255)
                                        tATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtgtg                                                  <  1:412475/93‑1 (MQ=255)
                                         ttCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtga                                                 <  1:1250286/92‑1 (MQ=255)
                                         ttCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtga                                                 <  2:3404783/92‑1 (MQ=255)
                                          tCGGCCGGCAGCGTATCGGTGCCAATCACCGCAAGGGCGCTTTTGGCGTTGCCGTTGCGGTCAGTGCGCAGGTTATGCCTCTCAATTTGtgat                                                <  2:3552227/93‑1 (MQ=255)
                                          tCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGTGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTGATGCCACTCAATTTGtgat                                                >  2:1143509/1‑93 (MQ=255)
                                          tCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtgat                                                >  2:2258013/1‑93 (MQ=255)
                                              ccgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaaa                                            >  1:1439725/1‑93 (MQ=255)
                                                gccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaaaaa                                          >  2:2199511/1‑93 (MQ=255)
                                                gccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaaaaa                                          <  2:2341919/93‑1 (MQ=255)
                                                 ccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAg                                         <  1:3719575/93‑1 (MQ=255)
                                                 ccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAg                                         >  2:2022599/1‑93 (MQ=255)
                                                 ccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAg                                         >  2:3195269/1‑93 (MQ=255)
                                                  cAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGt                                        >  2:1042421/1‑93 (MQ=255)
                                                   aGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTa                                       >  2:3337905/1‑93 (MQ=255)
                                                    gCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAt                                      >  1:4725692/1‑93 (MQ=255)
                                                    gCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAt                                      >  1:2931977/1‑93 (MQ=255)
                                                    gCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAt                                      <  2:339736/93‑1 (MQ=255)
                                                     cGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAtt                                     <  2:4358940/93‑1 (MQ=255)
                                                     cGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAtt                                     <  1:4462712/93‑1 (MQ=255)
                                                      gTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTa                                    >  1:4599835/1‑93 (MQ=255)
                                                       tATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTac                                   >  2:2483362/1‑93 (MQ=255)
                                                       tATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTac                                   >  2:4470545/1‑93 (MQ=255)
                                                        aTCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATGTGTGATGAAAAAGTATTaca                                  <  2:163707/93‑1 (MQ=255)
                                                         tCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTacac                                 >  2:1662332/1‑93 (MQ=255)
                                                         tCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTacac                                 <  1:2309898/93‑1 (MQ=255)
                                                         tCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTacac                                 >  2:1099132/1‑93 (MQ=255)
                                                         tCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTacac                                 >  2:4102611/1‑93 (MQ=255)
                                                           ggTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAg                               <  1:3692645/93‑1 (MQ=255)
                                                            gTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGc                              <  1:1664045/93‑1 (MQ=255)
                                                             tGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCg                             >  1:3984569/1‑93 (MQ=255)
                                                             tGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCg                             <  2:4119908/93‑1 (MQ=255)
                                                             tGCCAATCAGCGCAAGGGCGCTCTGGCCGTTGCCCTTGCTGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCg                             <  1:2476746/93‑1 (MQ=255)
                                                              gCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGt                            >  2:3047268/1‑93 (MQ=255)
                                                                cAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTTAATTTGTGATGAAAAAGTATTACACAGCGttt                          >  1:3384265/1‑93 (MQ=255)
                                                                 aaTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTc                         <  1:1812492/93‑1 (MQ=255)
                                                                 aaTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTc                         <  1:3127931/93‑1 (MQ=255)
                                                                 aaTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTc                         >  2:3174987/1‑93 (MQ=255)
                                                                 aaTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTc                         <  1:3158845/93‑1 (MQ=255)
                                                                  aTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCa                        >  2:2895589/1‑93 (MQ=255)
                                                                  aTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCa                        >  1:2573475/1‑93 (MQ=255)
                                                                  aTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCa                        <  1:2222017/93‑1 (MQ=255)
                                                                   tCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCaa                       <  2:1790407/93‑1 (MQ=255)
                                                                   tCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCaa                       >  2:3360981/1‑93 (MQ=255)
                                                                   tCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCaa                       >  1:638469/1‑93 (MQ=255)
                                                                   tCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCaa                       >  2:1378669/1‑93 (MQ=255)
                                                                    cAGCGCAAGGGTGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAAc                      >  1:3205850/1‑93 (MQ=255)
                                                                      gcgcAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGt                    >  1:806017/1‑93 (MQ=255)
                                                                       cgcAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTg                   <  2:553232/93‑1 (MQ=255)
                                                                       cgcAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTg                   <  1:2705984/93‑1 (MQ=255)
                                                                         cAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGc                 >  1:4331725/1‑93 (MQ=255)
                                                                         cAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGc                 <  1:1738902/93‑1 (MQ=255)
                                                                          aaGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCt                <  1:1638931/93‑1 (MQ=255)
                                                                          aaGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCt                >  2:2396135/1‑93 (MQ=255)
                                                                          aaGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCt                <  2:3733524/93‑1 (MQ=255)
                                                                             ggCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGc             >  2:1513678/1‑93 (MQ=255)
                                                                             ggCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGc             >  1:2221384/1‑93 (MQ=255)
                                                                                  ttttGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGAc        <  2:4037248/93‑1 (MQ=255)
                                                                                   tttGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGAcg       <  1:170158/93‑1 (MQ=255)
                                                                                       gCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACGCGCa   <  1:2918187/93‑1 (MQ=255)
                                                                                       gCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACGCGCa   <  1:353522/93‑1 (MQ=255)
                                                                                             |                                                                                        
GCTTCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACGCGCAT  >  NC_000913/4139575‑4139754

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