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,7790AC98.3% 803.2 / 1.9 232L67R (CTG→CGG) mreCcell wall structural complex MreBCD transmembrane component MreC
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base C (115/112);  total (118/114)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 234 total aligned reads displayed.

CTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGA  >  NC_000913/3399682‑3399875
                                                                                                 |                                                                                                
 tGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCcgcg                                                                                               <  1:4805602/100‑1 (MQ=255)
  gTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGt                                                                                              >  2:5083973/1‑100 (MQ=255)
  gTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGt                                                                                              >  2:3964930/1‑100 (MQ=255)
   tttGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTc                                                                                             <  1:3394747/100‑1 (MQ=255)
      gTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGc                                                                                          <  1:1390408/100‑1 (MQ=255)
      gTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGc                                                                                          <  1:2841833/100‑1 (MQ=255)
      gTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGc                                                                                          <  2:4795842/100‑1 (MQ=255)
      gTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGc                                                                                          <  1:226148/100‑1 (MQ=255)
       tATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCg                                                                                         >  1:5405049/1‑100 (MQ=255)
        aTTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGa                                                                                        >  2:2625/1‑100 (MQ=255)
          tctCCAAGCATTAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATCGGTCACGCGAGGCGCGCGTCTGCGATa                                                                                      <  1:4911723/98‑1 (MQ=255)
          tGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATa                                                                                      >  1:4239055/1‑100 (MQ=255)
          tGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATa                                                                                      <  1:3644066/100‑1 (MQ=255)
           gTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATAc                                                                                     <  2:2895298/100‑1 (MQ=255)
           gTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATAc                                                                                     >  1:3846398/1‑100 (MQ=255)
            tCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACg                                                                                    <  1:5357084/100‑1 (MQ=255)
            tCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATAcn                                                                                    <  1:347437/100‑2 (MQ=255)
             ccAAGAATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGATTTCACGTTCTAATTGCTCACGCGAGGCCCGCGTCTGCGATaaga                                                                                   >  2:4842710/1‑97 (MQ=255)
                 gCATCAGCAGTTCACTGTTGTTCAGCAACGGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATc                                                                               <  1:4509943/100‑1 (MQ=255)
                  cATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATcc                                                                              >  2:3421941/1‑100 (MQ=255)
                   aTCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCa                                                                             >  1:629740/1‑100 (MQ=255)
                    acagcagTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAg                                                                            <  2:3631236/99‑1 (MQ=255)
                     ccgcagTTCACTGTTTTTCCGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGc                                                                           <  2:72570/98‑1 (MQ=255)
                     cagcagTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGc                                                                           <  2:1722204/100‑1 (MQ=255)
                     cagcagTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGc                                                                           <  1:3053361/100‑1 (MQ=255)
                      agcagTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCa                                                                          >  1:1529430/1‑100 (MQ=255)
                       gcagTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCaa                                                                         >  2:4205289/1‑100 (MQ=255)
                        cagTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAAt                                                                        <  1:5251432/100‑1 (MQ=255)
                        cagTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAAt                                                                        <  1:1320642/100‑1 (MQ=255)
                          gTTCTCTGTTTTTTAGCAACACGCCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCTCGCGTCTGCGATACGCCATCCAGCAATTc                                                                      <  2:1618470/100‑1 (MQ=255)
                          gTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTc                                                                      >  1:3990978/1‑100 (MQ=255)
                          gTTCACTGTTTCTCAGCAGCAGTTCCAGACGTAACGCCCGGTTTTCAAGTTCTATTTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTc                                                                      <  2:3742265/100‑1 (MQ=255)
                           ttaaCTGTGTTTCAGCAACAATCCCTGACGTAACGCCCGATTTTCAAGTCCTCAATTGTCACTCGAGGCCCTTGTCTGCGATACGCCATCCAGCAATACa                                                                     <  2:898479/97‑1 (MQ=255)
                           ttCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCa                                                                     >  1:2026043/1‑100 (MQ=255)
                            tCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCACGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCAc                                                                    >  2:2782966/1‑100 (MQ=255)
                             cACTGTTTTTCAGCAACAGTTCGTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACg                                                                   >  1:4666617/1‑100 (MQ=255)
                              aCTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACga                                                                  >  2:2019564/1‑100 (MQ=255)
                               cTGTTTTTCAGCAACAGTTCCTGCCGTAACGCCCGGTCTTCAAGTTCTCATTGGTCACGCGAGGCCCCCGCCTGCGATCCGCCATCCAGCACTTCACgaa                                                                 >  2:3007230/1‑99 (MQ=255)
                                  tttttCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACgaggag                                                              >  2:327476/1‑100 (MQ=255)
                                   ttttCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGCGGAGc                                                             >  2:3793239/1‑100 (MQ=255)
                                      tCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCAtt                                                          >  2:5344083/1‑100 (MQ=255)
                                      tCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCAtt                                                          <  2:3334427/100‑1 (MQ=255)
                                        aGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTgg                                                        <  2:3187122/100‑1 (MQ=255)
                                        aGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTgg                                                        <  1:1965166/100‑1 (MQ=255)
                                        aGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTgg                                                        <  1:3164993/100‑1 (MQ=255)
                                          cAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGaa                                                      <  2:84684/100‑1 (MQ=255)
                                           aaCAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGaaa                                                     <  2:5371572/100‑1 (MQ=255)
                                              aGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACaa                                                  >  1:71049/1‑100 (MQ=255)
                                               gTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACaaa                                                 >  1:4100720/1‑100 (MQ=255)
                                                 tCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGt                                               >  2:1948346/1‑100 (MQ=255)
                                                 tCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGt                                               <  1:1501428/100‑1 (MQ=255)
                                                  ccTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAACATTGGAAACAAAGTa                                              >  1:1018153/1‑100 (MQ=255)
                                                   cTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAg                                             >  2:4571065/1‑100 (MQ=255)
                                                   cTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAg                                             <  2:1194523/100‑1 (MQ=255)
                                                   cTGACGTAACGCCCGGATTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAg                                             >  2:5329418/1‑100 (MQ=255)
                                                    tGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGa                                            <  2:3986007/100‑1 (MQ=255)
                                                    tGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGa                                            <  2:1059869/100‑1 (MQ=255)
                                                     gACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGaa                                           <  1:57256/100‑1 (MQ=255)
                                                      aCGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGaaa                                          >  2:2338608/1‑100 (MQ=255)
                                                        gTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAgg                                        >  2:1892130/1‑100 (MQ=255)
                                                        gTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAgg                                        >  1:4966336/1‑100 (MQ=255)
                                                         tAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGa                                       >  1:245742/1‑100 (MQ=255)
                                                         tAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGa                                       >  2:4013296/1‑100 (MQ=255)
                                                         tAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGa                                       >  1:1935794/1‑100 (MQ=255)
                                                            cGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTg                                    <  2:1310602/100‑1 (MQ=255)
                                                            cGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTg                                    >  1:4170516/1‑100 (MQ=255)
                                                            cGCCCGGATTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTg                                    >  2:4332846/1‑100 (MQ=255)
                                                             gCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGa                                   >  2:1124238/1‑100 (MQ=255)
                                                              nccGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGAc                                  >  1:4438037/2‑100 (MQ=255)
                                                              cccGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGAc                                  <  2:70282/100‑1 (MQ=255)
                                                              cccGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGAc                                  >  1:50812/1‑100 (MQ=255)
                                                              cccGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGAc                                  >  1:51169/1‑100 (MQ=255)
                                                                  gTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGAcggcg                              <  1:224332/100‑1 (MQ=255)
                                                                  gTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGAcggcg                              <  1:1486664/100‑1 (MQ=255)
                                                                  gTTTTCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGAcggcg                              <  1:4703412/100‑1 (MQ=255)
                                                                   ttttCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGAcggcgg                             <  1:3341639/100‑1 (MQ=255)
                                                                   ttttCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGAcggcgg                             <  1:1683236/100‑1 (MQ=255)
                                                                   ttttCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGAcggcgg                             <  1:3303894/100‑1 (MQ=255)
                                                                   ttttCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGAcggcgg                             <  1:1950545/100‑1 (MQ=255)
                                                                   ttttCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGAcggcgg                             <  2:1393749/100‑1 (MQ=255)
                                                                      tCAAGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTAt                          >  2:5018509/1‑100 (MQ=255)
                                                                        aaGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCACCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATcc                        >  1:4977507/1‑100 (MQ=255)
                                                                         aGTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCa                       <  1:5417019/100‑1 (MQ=255)
                                                                          gTTCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCat                      <  1:3716745/100‑1 (MQ=255)
                                                                           ttCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCata                     >  1:2357131/1‑100 (MQ=255)
                                                                           ttCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCata                     <  2:5052650/100‑1 (MQ=255)
                                                                           atCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCata                     <  2:3044004/99‑1 (MQ=255)
                                                                            tCTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCatat                    >  2:765586/1‑100 (MQ=255)
                                                                             cTAATTGGTCACGCGAGGCCCGCGTCTGGGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCAtata                   <  1:748930/100‑1 (MQ=255)
                                                                             cTAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCAtata                   <  2:414651/100‑1 (MQ=255)
                                                                              tAATTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATaa                  >  2:2233912/1‑100 (MQ=255)
                                                                               aaTTTGTCACGCGAGGCCCGCGTCAGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAg                 <  1:2175990/100‑1 (MQ=255)
                                                                               aaTTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAg                 >  1:3240411/1‑100 (MQ=255)
                                                                               aaTTGGTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAg                 >  2:1776318/1‑100 (MQ=255)
                                                                                   ggTCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACg             <  2:1568789/100‑1 (MQ=255)
                                                                                     tCACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGa           >  1:1807673/1‑100 (MQ=255)
                                                                                      cACGCGAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGAt          <  1:916257/100‑1 (MQ=255)
                                                                                        cgcgAGGCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGAttt        >  1:2414656/1‑100 (MQ=255)
                                                                                           gAGGCCCGCGTCTGCGATACGCCATCCAGCAATCCGCGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGAc     <  2:2377993/100‑1 (MQ=255)
                                                                                              gCCCGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGa  >  1:274504/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGA  >  NC_000913/3399682‑3399875

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