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.3% 948.5 / ‑1.1 270coding (231/330 nt)
coding (134/159 nt)
yijI
b4558
hypothetical protein
hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (136/131);  total (137/133)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.34e-01

Only 100 of 301 total aligned reads displayed.

GCTTCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACGCGCA  >  NC_000913/4139575‑4139753
                                                                                             |                                                                                      
 cTTCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGc                                                                                       <  1:3846883/93‑1 (MQ=255)
 cTTCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGc                                                                                       <  2:2857784/93‑1 (MQ=255)
    ccTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccct                                                                                    <  1:1073963/93‑3 (MQ=255)
    ccTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccct                                                                                    <  1:263348/93‑3 (MQ=255)
    ccTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccct                                                                                    <  1:5868820/93‑3 (MQ=255)
     cTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccctt                                                                                   >  2:3989693/1‑90 (MQ=255)
        ttCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCg                                                                                 >  1:1406406/1‑93 (MQ=255)
        ttCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCg                                                                                 <  2:4441730/93‑1 (MQ=255)
         tagtcgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCgg                                                                                >  1:4057603/3‑93 (MQ=255)
            tcgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCa                                                                             >  2:5871214/1‑93 (MQ=255)
              ggCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGt                                                                           <  1:2979323/93‑1 (MQ=255)
              ggCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGt                                                                           >  1:3601890/1‑93 (MQ=255)
              ggCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGt                                                                           >  2:3472377/1‑93 (MQ=255)
              ggCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGt                                                                           <  2:306353/93‑1 (MQ=255)
               gCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTg                                                                          >  1:2453462/1‑93 (MQ=255)
               gCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTg                                                                          <  2:736953/93‑1 (MQ=255)
               gCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTg                                                                          <  2:6046631/93‑1 (MQ=255)
               gCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTg                                                                          >  1:708526/1‑93 (MQ=255)
                ccccGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTgc                                                                         <  1:3193082/93‑1 (MQ=255)
                ccccGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTgc                                                                         >  1:1049539/1‑93 (MQ=255)
                     gAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAgg                                                                    <  2:2534888/93‑1 (MQ=255)
                     gAAATCAGTAACTGGCAGTTATCGCCCGCCAGAGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTTGCGTTGCCCTTGCGGTCAGTGCGCAgt                                                                    >  1:1461685/1‑92 (MQ=255)
                       aaTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGtt                                                                  <  2:2358842/93‑1 (MQ=255)
                       aaTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGtt                                                                  <  2:2938358/93‑1 (MQ=255)
                       aaTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGtt                                                                  <  1:651555/93‑1 (MQ=255)
                        aTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTa                                                                 >  1:6350065/1‑93 (MQ=255)
                        aTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTa                                                                 <  2:3661011/93‑1 (MQ=255)
                        aTCAGTAACTGGCAGTGATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTa                                                                 <  2:4864975/93‑1 (MQ=255)
                            gTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGcc                                                             <  1:5555/93‑1 (MQ=255)
                               aCTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACt                                                          <  1:2548769/93‑1 (MQ=255)
                                 tGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCa                                                        <  2:1032015/93‑1 (MQ=255)
                                 tGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCa                                                        >  2:5129914/1‑93 (MQ=255)
                                 tGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCa                                                        >  1:4456725/1‑93 (MQ=255)
                                 tGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCa                                                        <  2:1500475/93‑1 (MQ=255)
                                     aGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAAttt                                                    >  2:2791127/1‑93 (MQ=255)
                                      gTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtg                                                   >  2:330042/1‑93 (MQ=255)
                                        tATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtgtg                                                 <  2:2810760/93‑1 (MQ=255)
                                         aTCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtga                                                <  1:4004015/93‑1 (MQ=255)
                                          tCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtgat                                               >  2:3608687/1‑93 (MQ=255)
                                          tCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtgat                                               <  1:2313444/93‑1 (MQ=255)
                                           cgcccgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtgatg                                              >  2:4968470/1‑93 (MQ=255)
                                           cgcccgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtgatg                                              >  1:4552992/1‑93 (MQ=255)
                                            gcccgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtgatga                                             <  2:1553796/93‑1 (MQ=255)
                                              ccgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaaa                                           >  1:3295081/1‑93 (MQ=255)
                                                gccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaaaaa                                         >  1:3717326/1‑93 (MQ=255)
                                                  cAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGt                                       <  1:1598607/93‑1 (MQ=255)
                                                  cAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGt                                       <  1:822495/93‑1 (MQ=255)
                                                   aGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTa                                      <  2:5191721/93‑1 (MQ=255)
                                                   aGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTa                                      <  2:2127046/93‑1 (MQ=255)
                                                    gCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAt                                     >  2:460531/1‑93 (MQ=255)
                                                    gCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAt                                     >  2:3255003/1‑93 (MQ=255)
                                                    gCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAt                                     >  1:188991/1‑93 (MQ=255)
                                                     cGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAtt                                    <  1:5030925/93‑1 (MQ=255)
                                                     cGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAtt                                    >  1:4793150/1‑93 (MQ=255)
                                                      gTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTa                                   <  2:1146857/93‑1 (MQ=255)
                                                      gTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTa                                   <  1:432696/93‑1 (MQ=255)
                                                      gTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTa                                   >  2:5348304/1‑93 (MQ=255)
                                                       tATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTac                                  >  2:3807380/1‑93 (MQ=255)
                                                       tATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTac                                  >  1:5415811/1‑93 (MQ=255)
                                                       tATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTac                                  <  1:2734518/93‑1 (MQ=255)
                                                         tCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTacac                                <  2:1417430/93‑1 (MQ=255)
                                                          cGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTAcaca                               <  1:5223365/93‑1 (MQ=255)
                                                           ggTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGTTCAGTGCGCAGCTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAg                              <  1:3952242/93‑1 (MQ=255)
                                                            gTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGc                             >  2:1788414/1‑93 (MQ=255)
                                                            gTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGc                             <  1:6053424/93‑1 (MQ=255)
                                                            gTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGc                             >  2:6090228/1‑93 (MQ=255)
                                                            gTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGc                             <  2:6153159/93‑1 (MQ=255)
                                                              gCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGt                           <  2:1940235/93‑1 (MQ=255)
                                                              gCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGt                           >  2:1668394/1‑93 (MQ=255)
                                                              gCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGt                           <  2:4515363/93‑1 (MQ=255)
                                                               ccAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGtt                          >  2:3450566/1‑93 (MQ=255)
                                                                cAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGttt                         <  2:4317358/93‑1 (MQ=255)
                                                                cAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGttt                         <  1:931208/93‑1 (MQ=255)
                                                                 aaTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTc                        <  1:1272660/93‑1 (MQ=255)
                                                                 aaTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTc                        >  1:184739/1‑93 (MQ=255)
                                                                 aaTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTc                        >  2:4692353/1‑93 (MQ=255)
                                                                  aTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCa                       >  2:2041616/1‑93 (MQ=255)
                                                                   tCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAGGTATTACACAGCGTTTCaa                      <  1:6303104/93‑1 (MQ=255)
                                                                      gcgcAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGt                   <  1:2779155/93‑1 (MQ=255)
                                                                       cgcAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTg                  <  1:5832371/93‑1 (MQ=255)
                                                                       cgcAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTg                  <  1:2007156/93‑1 (MQ=255)
                                                                        gcAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTgg                 >  2:3185249/1‑93 (MQ=255)
                                                                        gcAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTgg                 >  2:1851029/1‑93 (MQ=255)
                                                                        gcAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTgg                 >  1:5918975/1‑93 (MQ=255)
                                                                          aaGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCt               >  1:2884750/1‑93 (MQ=255)
                                                                          aaGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCt               <  2:5443340/93‑1 (MQ=255)
                                                                           aGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCtt              >  1:2180935/1‑93 (MQ=255)
                                                                           aGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCtt              >  2:6338645/1‑93 (MQ=255)
                                                                           aGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCtt              <  2:1217382/93‑1 (MQ=255)
                                                                            gggCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTg             >  2:8361/1‑93 (MQ=255)
                                                                            gggCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTg             >  1:3377463/1‑93 (MQ=255)
                                                                            gggCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTg             >  1:4646990/1‑93 (MQ=255)
                                                                             ggCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGc            >  1:191577/1‑93 (MQ=255)
                                                                              gcgcTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGcc           <  1:565641/93‑1 (MQ=255)
                                                                              gcgcTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGcc           <  1:1184215/93‑1 (MQ=255)
                                                                                gcTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCgg         >  1:4125518/1‑93 (MQ=255)
                                                                                gcTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCgg         >  2:5973431/1‑93 (MQ=255)
                                                                                   tttGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGAcg      <  1:4899693/93‑1 (MQ=255)
                                                                                     tGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGAcgcg    <  2:3359083/93‑1 (MQ=255)
                                                                                       gCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACGCGCa  >  2:5705493/1‑93 (MQ=255)
                                                                                             |                                                                                      
GCTTCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACGCGCA  >  NC_000913/4139575‑4139753

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