Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 379,246 (G)10→9 intergenic (‑120/+47) yaiN ← / ← yaiO putative alpha helix chain/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913379,2370G.97.5% 504.4 / 7.5 122intergenic (‑111/+56)yaiN/yaiOputative alpha helix chain/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base . (73/43);  total (76/46)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-01

Only 100 of 156 total aligned reads displayed.

GTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCA  >  NC_000913/379145‑379327
                                                                                            |                                                                                           
gTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAg                                                                                             >  2:2736661/1‑93 (MQ=255)
 tATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgg                                                                                            <  2:3394266/93‑1 (MQ=255)
 tATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgg                                                                                            >  2:2067327/1‑93 (MQ=255)
 tATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgg                                                                                            >  1:2362047/1‑93 (MQ=255)
  aTTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggg                                                                                           >  2:2521773/1‑93 (MQ=255)
  aTTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggg                                                                                           >  2:4018445/1‑93 (MQ=255)
  aTTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggg                                                                                           >  2:2030238/1‑93 (MQ=255)
   ttCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggg                                                                                          >  2:187525/1‑93 (MQ=255)
      tataTGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggg                                                                                       <  2:1366448/93‑1 (MQ=255)
      tataTGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggg                                                                                       <  1:3879004/93‑1 (MQ=255)
      tataTGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggg                                                                                       >  2:2591540/1‑93 (MQ=255)
      tataTGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggg                                                                                       <  1:1800618/93‑1 (MQ=255)
       atatGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggggg                                                                                      <  1:1541491/93‑1 (MQ=255)
       atatGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGA‑‑GTGGTCGGGAAACAAAGAGGAAAgggggggggg                                                                                    >  1:4344914/1‑93 (MQ=255)
         atGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggggct                                                                                    >  2:1030172/1‑91 (MQ=255)
          tGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggggggct                                                                                   >  1:2129770/1‑91 (MQ=255)
          tGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggggggct                                                                                   <  2:1386928/93‑3 (MQ=255)
           gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggggggcta                                                                                  >  1:3678822/1‑90 (MQ=255)
           gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggggggcta                                                                                  <  2:1522557/93‑4 (MQ=255)
             cAATGCATACCCCCCTATAGTATATACCACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAAt                                                                               <  1:3673873/93‑1 (MQ=255)
             cAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAAt                                                                               <  1:1959530/93‑1 (MQ=255)
              aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATc                                                                              >  1:1977910/1‑93 (MQ=255)
              aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATc                                                                              >  1:114457/1‑93 (MQ=255)
                  cATCCCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCa                                                                          <  1:4384771/93‑1 (MQ=255)
                    tACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAgg                                                                        >  1:4020554/1‑93 (MQ=255)
                    tACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAgg                                                                        >  1:2858744/1‑93 (MQ=255)
                      ccccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGa                                                                      >  1:2080539/1‑93 (MQ=255)
                       cccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGaa                                                                     <  1:2795838/93‑1 (MQ=255)
                       cccccTATAGTAAATAACACTGGAGAATAAAATTTATCGGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGaa                                                                     <  1:1004037/93‑1 (MQ=255)
                         cccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCAATTGggaagggaag                                                                    >  2:1943808/1‑93 (MQ=255)
                          ccTATAGTATATAACACTGGAGAATAGAAGGTATCCGGTGAATGTGATCGGAAAACAAAGAGCAAAGGGGGGCTGCTAATCGGCAGGGAAGgc                                                                  <  2:549221/93‑1 (MQ=255)
                          ccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgc                                                                  >  2:1857888/1‑93 (MQ=255)
                           cTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGGAGgcc                                                                 >  1:4137978/1‑93 (MQ=255)
                             ataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCCAATCGGCAGGGGAGggcgc                                                               >  2:1424448/1‑93 (MQ=255)
                              taGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgcccac                                                              >  1:3555379/1‑90 (MQ=255)
                               aGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGccc                                                             >  1:4465590/1‑93 (MQ=255)
                               aGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGcccc                                                            >  1:2441511/1‑93 (MQ=255)
                                gTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCcccccc                                                            >  2:1055154/1‑93 (MQ=255)
                                gTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGcccc                                                            >  2:3101548/1‑93 (MQ=255)
                                gTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGcccc                                                            >  1:3249743/1‑93 (MQ=255)
                                 tatataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGGCGCCCCg                                                           >  2:1701517/1‑93 (MQ=255)
                                 tatataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCg                                                           >  1:1878764/1‑93 (MQ=255)
                                 tatataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCg                                                           >  2:3277109/1‑93 (MQ=255)
                                 tatataACACGGGAAAAAAAATTTTATGGGGTAAATGTGGTCGAAAAAAAAAGAGAAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCg                                                           <  2:4039762/93‑1 (MQ=255)
                                  atataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCgg                                                          >  1:363428/1‑93 (MQ=255)
                                  ataaaACACTGGAAAAAAAATTTTATCCGGTGAATGTGGCGGAAAAAAAAAGAGGAA‑‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGa                                                         <  2:3300456/93‑1 (MQ=255)
                                   tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGa                                                         >  1:2112252/1‑93 (MQ=255)
                                    ataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGAt                                                        <  1:404388/93‑1 (MQ=255)
                                    ataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCCGAt                                                        >  1:1883909/1‑93 (MQ=255)
                                     taACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATa                                                       >  1:567333/1‑93 (MQ=255)
                                      aaCACTGAAGAATAAAATTTATCCGGTGAATGGCGTCGTAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAg                                                      <  2:1216630/93‑1 (MQ=255)
                                         acgggAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgg                                                   <  2:3609392/90‑1 (MQ=255)
                                           tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggg                                                 >  1:3884969/1‑92 (MQ=255)
                                              agaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGggcggc                                              >  1:4556010/1‑93 (MQ=255)
                                               gaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCa                                             >  2:4245396/1‑93 (MQ=255)
                                                aaTAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGGAg                                            >  1:1025252/1‑93 (MQ=255)
                                                 aTAAATTTTTTCCGGTGAATGTCGGCGAAAAACAAAGGGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGaa                                          <  2:832564/93‑1 (MQ=255)
                                                 aTAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGa                                           >  2:215763/1‑93 (MQ=255)
                                                  tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGaa                                          >  2:4650596/1‑93 (MQ=255)
                                                  aaaaaTTTATCGGGTGAATGTGGTCGAAAAACAAAGGGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGaa                                          <  1:115002/92‑1 (MQ=255)
                                                    aaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAgg                                        >  1:1838669/1‑93 (MQ=255)
                                                    aaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCCGGTACGGGGAGGCAGaaag                                        >  2:3994623/1‑91 (MQ=255)
                                                     aaTTTATCCGGTGAATGTGGTCGAAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGa                                       <  1:3850021/93‑1 (MQ=255)
                                                       tttATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAAt                                     >  1:741328/1‑93 (MQ=255)
                                                       tttATCCGGGGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCggaaggaaa                                     >  2:4675926/1‑92 (MQ=255)
                                                         tATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGGAGAAGGAATCa                                   >  2:3031852/1‑93 (MQ=255)
                                                           tCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGa                                 >  2:589961/1‑93 (MQ=255)
                                                            ccGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGa                                 >  2:4436216/1‑93 (MQ=255)
                                                              ggTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAttt                             <  2:1497539/93‑1 (MQ=255)
                                                               gTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAttt                             >  1:2248295/1‑93 (MQ=255)
                                                                ggAATGTGTTGGGAAAGAAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTc                           <  2:943828/92‑1 (MQ=255)
                                                                 gAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTcc                          >  2:3174105/1‑93 (MQ=255)
                                                                  aaTGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCa                         >  1:4623829/1‑93 (MQ=255)
                                                                  aaTGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAg                        >  2:433764/1‑93 (MQ=255)
                                                                  aaTGTGGGGGGAAAAAAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCa                         <  2:4063397/93‑1 (MQ=255)
                                                                   aTGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGGTAGCGGGCGGCAGAAGGAATCAGAATTTCCa                         >  1:4115752/1‑93 (MQ=255)
                                                                    tgtcGTCGCAAAAAAAAGAGGAAAGGGGGTGGGCTAATCGGCAGGGAAGGCCGACCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAgg                       <  1:4405441/93‑1 (MQ=255)
                                                                      tgGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGGATCAGAATTTCCAgggc                     >  1:1331591/1‑91 (MQ=255)
                                                                      tgGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTc                     <  1:3534550/93‑1 (MQ=255)
                                                                      tgGACGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTc                     <  1:2845778/93‑1 (MQ=255)
                                                                         tCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGa                  >  1:1849639/1‑93 (MQ=255)
                                                                          cGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGAc                 <  1:1510259/93‑1 (MQ=255)
                                                                           ggAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACg                >  1:1989875/1‑93 (MQ=255)
                                                                            gAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCCGGTCAGACgg               >  2:1124190/1‑93 (MQ=255)
                                                                            gAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACgg               <  2:2222272/93‑1 (MQ=255)
                                                                            gAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAATTCCAGGTCAGACgg               >  1:3321271/1‑93 (MQ=255)
                                                                            gAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACggg              <  1:4010814/93‑1 (MQ=255)
                                                                              aaaCAAAGAGGAAAGGGGGGGGGTTAATCGGCAGGGAAGGCCCCCCCGGATAGCGGGGGGCACAGGGAATAAGAATTTTCAGGTGCGACGGGc             >  2:342818/1‑93 (MQ=255)
                                                                              aaaCAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGc             <  1:958402/93‑1 (MQ=255)
                                                                              aaaCAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGc             >  2:2167476/1‑93 (MQ=255)
                                                                              aaaCAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCCCCCCGGATAGCCGGCGGGAGAAGGAATCAGAATTTTCAGGGCAGACGGGc             >  1:1952922/1‑93 (MQ=255)
                                                                               accAAAGAGGAGAGGGGGGGTCCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGCCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt            <  2:2697436/91‑1 (MQ=255)
                                                                                aCAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCCCCCCGGATAGCGGGCGGCGGAAGGAATCAGAATTTCCCGGGTCAACGGGCtc           >  2:74177/1‑92 (MQ=255)
                                                                                   aaGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCaa        <  2:2135385/93‑1 (MQ=255)
                                                                                    agagGAAAGGGGGGGGGCTAATCGGCAGGGAAGGGCGCCCCGGGTAGCGGGCGGCAGGAGGAATGAGGATTTCCAAGTCAGAAGGGCTGGAAg       >  1:3806871/1‑93 (MQ=255)
                                                                                    agagGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg       >  1:3172720/1‑93 (MQ=255)
                                                                                    agagGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg       >  2:3173242/1‑93 (MQ=255)
                                                                                      agGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGtt     <  1:204724/93‑1 (MQ=255)
                                                                                      agGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGtt     <  1:1275671/93‑1 (MQ=255)
                                                                                         gagggggggggggCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGc   <  1:3189895/90‑1 (MQ=255)
                                                                                            |                                                                                           
GTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCA  >  NC_000913/379145‑379327

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