Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,992,339 G→C A160G (GCT→GGT)  pgsA ← phosphatidylglycerophosphate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,992,3390GC99.2% 422.7 / ‑4.0 118A160G (GCT→GGT) pgsAphosphatidylglycerophosphate synthetase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (59/58);  total (60/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 118 total aligned reads displayed.

GTTTGCTGAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCAGCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACT  >  NC_000913/1992243‑1992434
                                                                                                |                                                                                               
ggttGCTGAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCTCAGAGTCAGTACCGCACCca                                                                                              <  2:1847378/98‑1 (MQ=255)
gTTTGCTGAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCca                                                                                              >  1:2493187/1‑100 (MQ=255)
gTTTGCTGAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCca                                                                                              <  2:2619700/100‑1 (MQ=255)
gTTTGCTGAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCca                                                                                              >  2:374119/1‑100 (MQ=255)
gTTTGCTGAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCca                                                                                              <  1:1629477/100‑1 (MQ=255)
  ttGCTGAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCcac                                                                                             <  2:135892/99‑1 (MQ=255)
   tGCTGAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACaa                                                                                           <  1:2372970/100‑1 (MQ=255)
    gCTGAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACaaa                                                                                          <  2:855982/100‑1 (MQ=255)
    gCTGAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACaaa                                                                                          >  1:2890064/1‑100 (MQ=255)
     cTGAACATTACCCCGAAACGATCAATGATAAAGCAATTCTGCACGCTCAGCGCTCAAATATTGCAACATTGACCACAAAGTCAGTACCGCACCcaca                                                                                            <  2:2355925/97‑1 (MQ=255)
     cTGAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAg                                                                                         <  2:1146852/100‑1 (MQ=255)
      tGAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGa                                                                                        >  1:1716823/1‑100 (MQ=255)
      tGAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGa                                                                                        >  1:2003317/1‑100 (MQ=255)
      tGAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGa                                                                                        >  2:326007/1‑100 (MQ=255)
      tGAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGa                                                                                        >  1:1716799/1‑100 (MQ=255)
       gAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGaa                                                                                       >  1:395864/1‑100 (MQ=255)
       gAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGaa                                                                                       >  1:490535/1‑100 (MQ=255)
         aaaTTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGaaaa                                                                                     <  2:572339/100‑1 (MQ=255)
          aaTTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAg                                                                                    <  2:3003094/100‑1 (MQ=255)
          aaTTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAg                                                                                    <  1:387294/100‑1 (MQ=255)
          aaTTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAg                                                                                    <  1:1158525/100‑1 (MQ=255)
            ttACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTg                                                                                  <  2:866796/100‑1 (MQ=255)
             tACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGc                                                                                 <  2:1983535/100‑1 (MQ=255)
             tACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGc                                                                                 <  1:892891/100‑1 (MQ=255)
               cGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCa                                                                                <  2:2325340/99‑1 (MQ=255)
                gCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAAt                                                                              <  2:2828893/100‑1 (MQ=255)
                 ccGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATa                                                                             >  2:2195290/1‑100 (MQ=255)
                 ccGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATa                                                                             >  2:2778749/1‑100 (MQ=255)
                  cGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATAc                                                                            >  1:1333126/1‑100 (MQ=255)
                   gAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATAcc                                                                           <  2:2265322/100‑1 (MQ=255)
                      aCGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGc                                                                        <  1:1511936/100‑1 (MQ=255)
                       cGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGc                                                                        >  2:598438/1‑99 (MQ=255)
                       cGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCg                                                                       <  1:221136/100‑1 (MQ=255)
                       cGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCg                                                                       <  2:1716809/100‑1 (MQ=255)
                       cGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCg                                                                       <  2:1716785/100‑1 (MQ=255)
                       cGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCg                                                                       <  1:358353/100‑1 (MQ=255)
                        gATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGt                                                                      >  1:622498/1‑100 (MQ=255)
                        gATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTAAGTACCGCACCCACAAAGAAAAGTGCAATACCAGCg                                                                       >  2:947796/1‑99 (MQ=255)
                          tCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTAc                                                                    <  2:1798164/100‑1 (MQ=255)
                           cACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGCAAAGTGCAATACCGGCGTACt                                                                   >  2:541495/1‑100 (MQ=255)
                            aCTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTc                                                                  >  1:184609/1‑100 (MQ=255)
                            aCTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTc                                                                  <  1:99156/100‑1 (MQ=255)
                            aCTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTc                                                                  >  2:2238280/1‑100 (MQ=255)
                               gATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAAc                                                               >  1:828599/1‑100 (MQ=255)
                               gATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAAc                                                               >  1:199074/1‑100 (MQ=255)
                               gATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAAc                                                               <  1:2313036/100‑1 (MQ=255)
                                aTCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAAc                                                               >  2:2470636/1‑99 (MQ=255)
                                 tCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAAccc                                                             <  1:2216063/100‑1 (MQ=255)
                                  cAAGCAAATCTGCGCGCGCAGCGCTCAAATAATGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCa                                                            >  1:2037487/1‑100 (MQ=255)
                                  cAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCa                                                            >  1:516108/1‑100 (MQ=255)
                                  cAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCa                                                            >  1:2851445/1‑100 (MQ=255)
                                   aaGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGTAATACCGGCGTACTCAAccc                                                             >  1:1990342/1‑98 (MQ=255)
                                    aGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCaaa                                                          >  2:2070608/1‑100 (MQ=255)
                                    aGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGACGTACTCAACCCaaa                                                          >  1:2988119/1‑100 (MQ=255)
                                     gCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAAt                                                         >  1:1061942/1‑100 (MQ=255)
                                     gCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAAt                                                         <  2:71262/100‑1 (MQ=255)
                                      cAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATg                                                        >  2:1845471/1‑100 (MQ=255)
                                           cTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCgg                                                   <  2:39959/100‑1 (MQ=255)
                                            tGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGa                                                  >  2:246124/1‑100 (MQ=255)
                                            tGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGa                                                  >  2:2793641/1‑100 (MQ=255)
                                             gCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGAc                                                 <  1:2027838/100‑1 (MQ=255)
                                              cACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACg                                                >  1:1236609/1‑100 (MQ=255)
                                              cACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACg                                                >  1:1619734/1‑100 (MQ=255)
                                              cACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACg                                                >  2:2228991/1‑100 (MQ=255)
                                                cgcgCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGcc                                              >  2:2354220/1‑100 (MQ=255)
                                                 gcgcAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCca                                             <  1:1031337/100‑1 (MQ=255)
                                                    cAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCAcag                                          <  2:1061932/100‑1 (MQ=255)
                                                      gcgcTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCAcagca                                        <  1:834281/100‑1 (MQ=255)
                                                       cgcTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCAcagcag                                       >  2:115696/1‑100 (MQ=255)
                                                       cgcTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCAcagcag                                       <  1:703759/100‑1 (MQ=255)
                                                       cgcTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCAcagcag                                       <  1:979731/100‑1 (MQ=255)
                                                        gcTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACagcagc                                      <  2:2765652/100‑1 (MQ=255)
                                                         cTCAAATATTGCAACATTGACCACAGAGGCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGcc                                     <  2:2851419/100‑1 (MQ=255)
                                                           cAAATATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCa                                    >  2:1041549/1‑99 (MQ=255)
                                                            aaaTATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATg                                  <  2:490530/100‑1 (MQ=255)
                                                             aaTATTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGc                                 >  2:1481617/1‑100 (MQ=255)
                                                                atTGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCaa                              >  1:2325540/1‑100 (MQ=255)
                                                                 ttGCAACATTGACCACAGAGTCAGTACCGCACACACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAAc                             >  1:2871485/1‑100 (MQ=255)
                                                                  tGCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGTCAACg                            >  2:719415/1‑100 (MQ=255)
                                                                   gCAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGc                           >  1:911260/1‑100 (MQ=255)
                                                                    cAACATTGACTACAGAGTCAGTACCGCACCCACAAAGAAAAGTGTAATACCGGCGTACTCAACCCATATATTCGGACGCCACAGCAGCCATGCCAACGcn                          <  2:1990324/100‑2 (MQ=255)
                                                                    cAACATTGACCACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGcc                          >  2:1040790/1‑100 (MQ=255)
                                                                             ccACAGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTgg                 >  2:855640/1‑100 (MQ=255)
                                                                              cacaGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTggg                >  2:2213940/1‑100 (MQ=255)
                                                                               acaGAGTCAGTACCGCAGCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGc               >  1:1379100/1‑100 (MQ=255)
                                                                               acaGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGc               <  2:831087/100‑1 (MQ=255)
                                                                               acaGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGc               <  2:2405119/100‑1 (MQ=255)
                                                                                caGAGTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCa              >  1:2637004/1‑100 (MQ=255)
                                                                                    gTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTgg          >  1:1760146/1‑100 (MQ=255)
                                                                                    gTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTgg          <  2:1259111/100‑1 (MQ=255)
                                                                                    gTCAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTgg          >  2:2921980/1‑100 (MQ=255)
                                                                                      cAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGtt        >  2:2117647/1‑100 (MQ=255)
                                                                                      cAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGtt        >  1:2443575/1‑100 (MQ=255)
                                                                                      cAGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGtt        >  2:1739087/1‑100 (MQ=255)
                                                                                       aGTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGttt       >  1:2303068/1‑100 (MQ=255)
                                                                                        gTACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGtttt      <  2:1236598/100‑1 (MQ=255)
                                                                                         tACCGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTc     <  2:2789727/100‑1 (MQ=255)
                                                                                           ccGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCAc   <  1:2070627/100‑1 (MQ=255)
                                                                                           ccGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCAc   >  1:2025013/1‑100 (MQ=255)
                                                                                            cGCACCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACt  <  1:1619467/100‑1 (MQ=255)
                                                                                                |                                                                                               
GTTTGCTGAAAATTACGCCGAAACGATCACTGATCAAGCAAATCTGCACGCGCAGCGCTCAAATATTGCAACATTGACCACAGAGTCAGTACCGCAGCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACT  >  NC_000913/1992243‑1992434

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