Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,399,779 A→C L67R (CTG→CGG)  mreC ← cell wall structural complex MreBCD transmembrane component MreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,399,7790AC99.5% 671.4 / NA 191L67R (CTG→CGG) mreCcell wall structural complex MreBCD transmembrane component MreC
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (100/90);  total (100/91)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 191 total aligned reads displayed.

CTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTG  >  NC_000913/3399682‑3399874
                                                                                                 |                                                                                               
cTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGc                                                                                               >  1:2646768/1‑100 (MQ=255)
  gTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGt                                                                                             >  1:3119772/1‑100 (MQ=255)
   tttGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTc                                                                                            >  2:2770237/1‑100 (MQ=255)
    ttGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTCATTGGTCACGCGAGGCCCGCGTCt                                                                                           >  2:3484555/1‑100 (MQ=255)
      gTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGc                                                                                         >  1:2482343/1‑100 (MQ=255)
       tATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCg                                                                                        >  1:570479/1‑100 (MQ=255)
       tATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCg                                                                                        >  1:401408/1‑100 (MQ=255)
       tATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCg                                                                                        >  2:2845542/1‑100 (MQ=255)
       tATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCg                                                                                        <  1:3422438/100‑1 (MQ=255)
           gTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTTACGCGAGGCCCGCGTCTGCGATAc                                                                                    >  2:768479/1‑100 (MQ=255)
           gTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATAc                                                                                    >  1:3590323/1‑100 (MQ=255)
           gTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATAc                                                                                    <  2:891646/100‑1 (MQ=255)
           gTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATAc                                                                                    <  1:557898/100‑1 (MQ=255)
           gTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATAc                                                                                    <  1:1375178/100‑1 (MQ=255)
             ccAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGc                                                                                  <  2:354350/100‑1 (MQ=255)
              cAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGcc                                                                                 >  2:1358287/1‑100 (MQ=255)
              cAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGcc                                                                                 >  1:3933317/1‑100 (MQ=255)
              cAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGcc                                                                                 >  2:1334254/1‑100 (MQ=255)
              cAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGcc                                                                                 >  1:1232613/1‑100 (MQ=255)
                aGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCAt                                                                               <  1:1001703/100‑1 (MQ=255)
                   aTCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCa                                                                            >  1:2503069/1‑100 (MQ=255)
                   aTCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCa                                                                            <  1:605525/100‑1 (MQ=255)
                   aTCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCa                                                                            >  2:1285793/1‑100 (MQ=255)
                    tCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAg                                                                           <  1:3367558/100‑1 (MQ=255)
                    tCAGCAGTTCACTGGTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAg                                                                           <  2:1502857/100‑1 (MQ=255)
                     cagcagTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGc                                                                          >  2:1024190/1‑100 (MQ=255)
                     cagcagTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGc                                                                          <  2:2637309/100‑1 (MQ=255)
                     cagcagTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGc                                                                          <  1:3018658/100‑1 (MQ=255)
                     cagcagTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGc                                                                          <  2:3312615/100‑1 (MQ=255)
                      agcagTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCa                                                                         >  2:3285785/1‑100 (MQ=255)
                       gcagTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCaa                                                                        >  1:99102/1‑100 (MQ=255)
                       gcagTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCaa                                                                        >  1:2273374/1‑100 (MQ=255)
                        cagTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAAt                                                                       <  1:458917/100‑1 (MQ=255)
                          gTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTc                                                                     >  2:1260899/1‑100 (MQ=255)
                           ttCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCa                                                                    >  1:2934790/1‑100 (MQ=255)
                            tCACTGTTTTTCAGCAACAGTCCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCAc                                                                   >  2:3162422/1‑100 (MQ=255)
                             cACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACg                                                                  >  1:2893031/1‑100 (MQ=255)
                               cTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACgag                                                                >  2:1803134/1‑100 (MQ=255)
                                 gTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACgagga                                                              >  1:3858156/1‑100 (MQ=255)
                                  tttttCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACgaggag                                                             >  2:487031/1‑100 (MQ=255)
                                   ttttCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGc                                                            <  1:2143477/100‑1 (MQ=255)
                                     ttCAGCAACAGTTCCTGACGTAACGCCCGTTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCAt                                                          <  2:1628667/100‑1 (MQ=255)
                                     ttCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCAt                                                          <  1:3052420/100‑1 (MQ=255)
                                      tCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCAtt                                                         <  2:2322647/100‑1 (MQ=255)
                                        aGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTgg                                                       >  2:1856033/1‑100 (MQ=255)
                                         gCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGa                                                      >  1:1469311/1‑100 (MQ=255)
                                         gCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGa                                                      >  2:56753/1‑100 (MQ=255)
                                         gCAACAGCTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGa                                                      >  2:1735084/1‑100 (MQ=255)
                                          cAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGaa                                                     >  2:1454977/1‑100 (MQ=255)
                                          cAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGaa                                                     <  2:2652211/100‑1 (MQ=255)
                                            aCAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAAc                                                   >  1:299284/1‑100 (MQ=255)
                                             cAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACa                                                  >  1:2673912/1‑100 (MQ=255)
                                              aGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAAGTGGTCACGCGAGGCCCGCGTCTGCAATACGCCATCCAGCAAATCACGAGGAGCATTGGAAACaa                                                 <  2:439018/100‑1 (MQ=255)
                                                ttCCTGACGTAACGCCCGGATTTCAAGTTCTAATTGGTCACGCGAAGCCCGCGTCTGCGATACGCCATCCAGCAAGTCACGAGGAGCATTGGAAACAAAg                                               >  2:1893105/1‑100 (MQ=255)
                                                ctCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAg                                               <  1:1763798/99‑1 (MQ=255)
                                                  ccTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTa                                             >  2:1748590/1‑100 (MQ=255)
                                                  ccTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTa                                             >  1:2819594/1‑100 (MQ=255)
                                                   cTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGCGGAGCATTGGAAACAAAGTAg                                            >  2:86338/1‑100 (MQ=255)
                                                   cTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAg                                            >  1:1025895/1‑100 (MQ=255)
                                                   cTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAg                                            >  1:314163/1‑100 (MQ=255)
                                                   cTGACGGAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAg                                            <  2:2444198/100‑1 (MQ=255)
                                                    tGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGa                                           >  1:740639/1‑100 (MQ=255)
                                                     gACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGaa                                          >  2:112036/1‑100 (MQ=255)
                                                      aCGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGaan                                         <  1:3241889/100‑2 (MQ=255)
                                                       cGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAg                                        <  1:57743/100‑1 (MQ=255)
                                                          aaCGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCAGCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGAc                                     <  1:3475937/100‑1 (MQ=255)
                                                             gCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGa                                  <  2:3230105/100‑1 (MQ=255)
                                                             gCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGa                                  <  1:3532635/100‑1 (MQ=255)
                                                             gCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGa                                  >  2:359627/1‑100 (MQ=255)
                                                              ctcGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGCTACGCGATCCAGCAATTCACGAGGAGTATTGGATACAAAGTAGAAAGGACTGAc                                 <  2:3985417/98‑1 (MQ=255)
                                                              cccGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGAc                                 >  1:2436839/1‑100 (MQ=255)
                                                                   ttttCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGAcggcgg                            <  2:2399442/100‑1 (MQ=255)
                                                                   ttttCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGAcggcgg                            <  1:1914058/100‑1 (MQ=255)
                                                                     ttCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTa                          <  2:2799023/100‑1 (MQ=255)
                                                                      tCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTAt                         <  1:380073/100‑1 (MQ=255)
                                                                        aaGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATcc                       >  1:2762074/1‑100 (MQ=255)
                                                                         aGTTCTACTTGGTCCCGCGACGCCCGCGTCAGCGATACGCCATCCCCCAATTCACCAGCGGCATCGCTAACACAGTAGCAAGGACTGACGGCGGTATCCa                      >  2:893966/1‑100 (MQ=255)
                                                                         aGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCa                      >  1:2908361/1‑100 (MQ=255)
                                                                          gTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCat                     <  2:2735726/100‑1 (MQ=255)
                                                                          gTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCat                     <  1:447204/100‑1 (MQ=255)
                                                                           ttCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCata                    <  1:1680481/100‑1 (MQ=255)
                                                                             cTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCAtata                  <  1:2447189/100‑1 (MQ=255)
                                                                             cTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCAtata                  >  1:2942301/1‑100 (MQ=255)
                                                                              tacttGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATaa                 >  2:2278031/4‑100 (MQ=255)
                                                                               aaTTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAg                <  2:968028/100‑1 (MQ=255)
                                                                               aaTTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAg                <  2:2358704/100‑1 (MQ=255)
                                                                                aTTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGt               >  2:1339361/1‑100 (MQ=255)
                                                                                 ttGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTa              <  1:1533491/100‑1 (MQ=255)
                                                                                  tGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTAc             >  1:1830898/1‑100 (MQ=255)
                                                                                   ggTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACg            <  2:2472633/100‑1 (MQ=255)
                                                                                    gTCACGCGAGGCCCACGTCTGCGATAAGCCATCCAGCAATTCACGAGAAGCATTTGAAACAAAGTAGAAAGGACTGACGGCGGTATTCATATAATTACgg           >  2:3733663/1‑100 (MQ=255)
                                                                                       aCGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGAtt        <  2:3912594/100‑1 (MQ=255)
                                                                                       aCGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGAtt        >  1:786459/1‑100 (MQ=255)
                                                                                        cgcgAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGAttt       >  1:2577896/1‑100 (MQ=255)
                                                                                         gcgAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTg      <  2:2106726/100‑1 (MQ=255)
                                                                                          cgAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGa     >  2:2345217/1‑100 (MQ=255)
                                                                                           gAGGCCCGCTTCTGCGATACGCCATCCAGCAATTCACGAGGAACATTGGAATCAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGAc    <  2:433908/100‑1 (MQ=255)
                                                                                           gAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGAc    <  1:158688/100‑1 (MQ=255)
                                                                                           gAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGAc    <  2:1295671/100‑1 (MQ=255)
                                                                                            aGGCCCGCGGCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACt   <  1:511461/100‑1 (MQ=255)
                                                                                                 |                                                                                               
CTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTG  >  NC_000913/3399682‑3399874

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: ‑