Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,322,698 A→G 16.0% intergenic (‑193/+35) secG ← / ← glmM preprotein translocase membrane subunit/phosphoglucosamine mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,322,6980AG16.0% 480.3 / 76.3 210intergenic (‑193/+35)secG/glmMpreprotein translocase membrane subunit/phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base A (89/84);  new base G (13/20);  total (103/107)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

Only 100 of 214 total aligned reads displayed.

CTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGC  >  NC_000913/3322599‑3322797
                                                                                                   |                                                                                                   
cTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGa                                                                                                     >  1:310800/1‑100 (MQ=255)
 tAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGAc                                                                                                    <  2:3789868/100‑1 (MQ=255)
    ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAAt                                                                                                 >  1:323032/1‑100 (MQ=255)
    ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAAt                                                                                                 >  2:4053803/1‑100 (MQ=255)
     cAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATc                                                                                                >  2:2098265/1‑100 (MQ=255)
      aaaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATTTGTTCACTGAAAAAAAAGGCGACAATcg                                                                                               <  2:629034/100‑1 (MQ=255)
       aaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgc                                                                                              <  1:706988/100‑1 (MQ=255)
       aaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgc                                                                                              <  2:2728116/100‑1 (MQ=255)
       aaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgc                                                                                              >  1:2940855/1‑100 (MQ=255)
       aaGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgc                                                                                              >  1:995320/1‑100 (MQ=255)
           gAATGACCTTCGCAAGGGCAATTTAATCGCATTGTATCAATTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCgccgcc                                                                                          <  2:1989808/100‑1 (MQ=255)
           gAAGGACCAGCGCAAGAGCAATGTCATAGGATTGTAGCAAAGAGAGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCgccgcc                                                                                          <  2:1486954/100‑1 (MQ=255)
               gACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCtttt                                                                                      >  1:3053982/1‑100 (MQ=255)
                aCCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttt                                                                                     <  2:3324947/100‑1 (MQ=39)
                aCCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCttttt                                                                                     <  1:3052876/100‑1 (MQ=255)
                 ccTTCGCAAGTGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTCAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                    <  2:2765521/100‑1 (MQ=255)
                 ccTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCtttttt                                                                                    <  2:253052/100‑2 (MQ=38)
                 ccTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                    >  2:3911524/1‑100 (MQ=255)
                 ccTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                                    <  1:1992629/100‑1 (MQ=255)
                  cTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttc                                                                                   <  2:713101/100‑3 (MQ=38)
                   ttCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttct                                                                                  <  1:4182178/100‑4 (MQ=38)
                   ttCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttct                                                                                  <  1:782722/100‑4 (MQ=38)
                    tCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttctc                                                                                 <  1:733647/100‑5 (MQ=37)
                    tCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttctc                                                                                 <  1:3907446/100‑5 (MQ=37)
                     cgtaaGGGCAATTTTATTGCATTGTATCAATTGCGAAAAAAAACAGCAGCACCTCTGTGTTCGCTGAGAAAAAAGGCTGCGATTGTCGTCttttttctca                                                                                <  2:1447269/97‑6 (MQ=255)
                     cGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttctca                                                                                <  1:3217518/100‑6 (MQ=37)
                     cGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCa                                                                                >  1:3190774/1‑100 (MQ=255)
                      gCAAGGGCAATTTTATCGCATTGTATCAACTGTGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttctcag                                                                               <  1:3333205/100‑7 (MQ=25)
                      gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttctcag                                                                               <  2:67483/100‑7 (MQ=37)
                       cAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGt                                                                              <  2:59529/100‑1 (MQ=255)
                       cAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGt                                                                              >  2:1341860/1‑100 (MQ=255)
                         aGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTa                                                                            <  1:3372834/100‑1 (MQ=255)
                         aGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTa                                                                            <  1:1379741/100‑1 (MQ=255)
                          gggCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTAt                                                                           >  1:3667403/1‑100 (MQ=255)
                               aTTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAAc                                                                      >  1:1929060/1‑100 (MQ=255)
                                ttttATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACg                                                                     <  1:3070189/100‑1 (MQ=255)
                                   tATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCtt                                                                  <  2:241532/100‑1 (MQ=255)
                                    aTCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCttt                                                                 <  1:1061356/100‑1 (MQ=255)
                                     tCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTa                                                                >  2:2292613/1‑100 (MQ=255)
                                        cATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAAc                                                             <  1:2905425/100‑1 (MQ=255)
                                        cATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTAATCTAACGCTTTAAAc                                                             >  2:2452907/1‑100 (MQ=255)
                                         aTTGTATCAACTGCGGAAAAAAACAGAAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTGAGCCAGTTATCTAACGCTTTAAACg                                                            >  1:1594967/1‑100 (MQ=255)
                                           tGTATCAATTGCGGAAAAAAACAGCATAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGc                                                          <  1:2491842/100‑1 (MQ=255)
                                                cAACTGCGGCAAAAAACAGCAAAACCGATGTGTTCGCTTAGAAAAAAGGCGACAATCGCCGCCTTTTTATCCAGTTATCTAACGCTTTAAACGGCTTTTa                                                     >  2:3232151/1‑100 (MQ=255)
                                                cAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTa                                                     <  2:3969915/100‑1 (MQ=255)
                                                  aCTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACt                                                   <  1:2346953/100‑1 (MQ=255)
                                                   cTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTg                                                  <  1:3943065/100‑1 (MQ=255)
                                                   cTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTg                                                  <  1:1543436/100‑1 (MQ=255)
                                                   cTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTg                                                  <  1:16490/100‑1 (MQ=255)
                                                    tGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGTCAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTg                                                  <  1:1046149/99‑1 (MQ=255)
                                                     gCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCa                                                >  2:1720979/1‑100 (MQ=255)
                                                        gAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATc                                              >  2:3111109/1‑99 (MQ=255)
                                                        gAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCg                                             >  2:3024938/1‑100 (MQ=255)
                                                        gAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCg                                             >  1:2093797/1‑100 (MQ=255)
                                                         aaaaaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCgg                                            <  1:143962/100‑1 (MQ=255)
                                                         aaaaaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCgg                                            <  1:4069882/100‑1 (MQ=255)
                                                          aaaaaaCCGCAAAACGGATGTGTGCGCTGAGAAAAAAGGCGGCAATCGCCGCCTTTTTAGCCAGTGATCTAACGCTTTAAACGGCTTTTACTGCATCGGc                                           <  2:2892073/100‑1 (MQ=255)
                                                             aaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGAt                                        >  2:730186/1‑100 (MQ=255)
                                                               aCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGc                                      <  1:4034086/100‑1 (MQ=255)
                                                                 aGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCgg                                    <  1:2980660/100‑1 (MQ=255)
                                                                  gCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCgg                                    >  1:3819920/1‑99 (MQ=255)
                                                                  gCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgt                                   >  1:221849/1‑100 (MQ=255)
                                                                  gCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgt                                   >  1:470468/1‑100 (MQ=255)
                                                                   cAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgtg                                  <  1:744304/100‑1 (MQ=255)
                                                                      agccGATGTGTTCGCTGAGAAAAATGGCGACAATCGCCGCCTGTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGc                               <  2:3996274/98‑1 (MQ=255)
                                                                      aaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGc                               >  1:178460/1‑100 (MQ=255)
                                                                        ccGATGTGTTCGCTGAGAAAAAAGGCGCCAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCa                              <  1:2085835/99‑1 (MQ=255)
                                                                             gtgtTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGTGGGGTGCAAATTCa                        >  1:2281554/1‑100 (MQ=255)
                                                                             gtgtTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCa                        <  2:4089017/100‑1 (MQ=255)
                                                                              tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                       <  2:2987487/100‑1 (MQ=255)
                                                                              tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                       <  2:4150400/100‑1 (MQ=255)
                                                                                 tCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCa                    <  1:4046630/100‑1 (MQ=255)
                                                                                 tCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCACCGGCGATGCGGTGTGCAAATTCAGTCa                    >  2:3518975/1‑100 (MQ=255)
                                                                                   gCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCAcc                  <  2:834112/100‑1 (MQ=255)
                                                                                   gCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCAcc                  <  1:2301765/100‑1 (MQ=255)
                                                                                    cTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATATAACGCGTTAAACGGCTTCTACTGCAACCGCGATGCGGTGCGCAAATTCAGTcgcat                 >  2:3708940/1‑96 (MQ=255)
                                                                                     tGAGAAAAAAGGCGACAATCGCCGCCTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc               >  2:756195/1‑100 (MQ=255)
                                                                                      ggctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc               >  2:2124277/5‑100 (MQ=37)
                                                                                      ggctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc               >  2:534658/5‑100 (MQ=37)
                                                                                      gagaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc               <  2:1914451/100‑1 (MQ=255)
                                                                                       gctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcg              >  1:577326/4‑100 (MQ=38)
                                                                                       gctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcg              >  1:930302/4‑100 (MQ=38)
                                                                                       agaAAAAAGGCGACAATCGCCGTCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgc               >  2:432989/1‑99 (MQ=255)
                                                                                        ctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcgc             >  2:162695/3‑100 (MQ=38)
                                                                                        ctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcgc             >  2:490998/3‑100 (MQ=38)
                                                                                        ctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTACCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcgc             >  2:2239496/3‑100 (MQ=37)
                                                                                        caaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcgc             >  1:2851612/2‑100 (MQ=38)
                                                                                          aaaaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCtt           >  2:1322730/1‑100 (MQ=255)
                                                                                          aaaaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCtt           >  1:2577689/1‑100 (MQ=255)
                                                                                            aaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCg         <  1:1166159/100‑1 (MQ=255)
                                                                                            aaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCg         <  1:2667306/100‑1 (MQ=255)
                                                                                            aaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCg         <  2:1527400/100‑1 (MQ=255)
                                                                                             aaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCCCCTGCGCTTCGt        >  2:1094413/1‑100 (MQ=255)
                                                                                             aaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGt        <  2:1667309/100‑1 (MQ=255)
                                                                                               ggCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCt      >  1:1224186/1‑100 (MQ=255)
                                                                                                 cGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTc    >  2:1603857/1‑100 (MQ=255)
                                                                                                 cGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTc    >  1:916604/1‑100 (MQ=255)
                                                                                                 cGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTc    >  1:2405202/1‑100 (MQ=255)
                                                                                                  gACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCg   <  1:3861168/100‑1 (MQ=255)
                                                                                                  gACAATCGCCGCCTTTTTAGCCAGTGATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCg   <  2:3896623/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGC  >  NC_000913/3322599‑3322797

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