Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,440,503 A→G V411A (GTC→GCC)  prlA ← putative ATPase subunit of translocase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,440,5030AG87.6% 756.2 / 91.2 298V411A (GTC→GCC) prlAputative ATPase subunit of translocase
Reads supporting (aligned to +/- strand):  ref base A (19/18);  new base G (118/143);  total (137/161)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 317 total aligned reads displayed.

GTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATA  >  NC_000913/3440412‑3440595
                                                                                           |                                                                                            
 tAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGgca                                                                                            >  2:5762576/1‑93 (MQ=255)
   gCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGgcaac                                                                                          <  2:4020292/93‑1 (MQ=255)
   gCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGacaac                                                                                          >  2:2272152/1‑93 (MQ=255)
     cTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGcaacaa                                                                                        >  2:719111/1‑93 (MQ=255)
      tttCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCaacaac                                                                                       >  1:1269020/1‑93 (MQ=255)
        tCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGa                                                                                     <  1:712238/93‑1 (MQ=255)
        tCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGa                                                                                     <  2:5194702/93‑1 (MQ=255)
          aGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATa                                                                                   >  1:827728/1‑93 (MQ=255)
          aGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATa                                                                                   <  1:4156455/93‑1 (MQ=255)
          aGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATa                                                                                   >  2:6195419/1‑93 (MQ=255)
          aGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATa                                                                                   >  2:812038/1‑93 (MQ=255)
           ggTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATaa                                                                                  >  1:4294989/1‑93 (MQ=255)
           ggTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATaa                                                                                  <  1:658824/93‑1 (MQ=255)
             ttCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGc                                                                                >  1:6819416/1‑93 (MQ=255)
             ttCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGc                                                                                >  2:6250320/1‑93 (MQ=255)
             ttCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGc                                                                                <  2:779443/93‑1 (MQ=255)
              tCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCa                                                                               >  2:2410433/1‑93 (MQ=255)
                gCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGt                                                                             <  2:4985900/93‑1 (MQ=255)
                  cttcttCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGa                                                                           >  2:3254026/1‑93 (MQ=255)
                   ttcttcAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAg                                                                          >  1:3500043/1‑93 (MQ=255)
                   ttcttcAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAg                                                                          <  2:2314211/93‑1 (MQ=255)
                   ttcttcAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAg                                                                          >  1:2094902/1‑93 (MQ=255)
                   ttcttcAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAg                                                                          <  2:1007294/93‑1 (MQ=255)
                    tcttcAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAgg                                                                         <  2:7426545/93‑1 (MQ=255)
                     cttcAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGt                                                                        <  1:144227/93‑1 (MQ=255)
                      ttcAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTc                                                                       >  2:5833566/1‑93 (MQ=255)
                       tcAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTcc                                                                      >  2:609775/1‑93 (MQ=255)
                       tcAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTcc                                                                      >  1:2212578/1‑93 (MQ=255)
                         aaTGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCa                                                                    >  1:1026033/1‑93 (MQ=255)
                         aaTGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCa                                                                    <  1:738508/93‑1 (MQ=255)
                         aaTGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCa                                                                    >  1:6308481/1‑93 (MQ=255)
                          aTGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCAc                                                                   >  2:462843/1‑93 (MQ=255)
                          aTGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCAc                                                                   >  1:2496742/1‑93 (MQ=255)
                           tGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCAcc                                                                  >  2:1631349/1‑93 (MQ=255)
                           tGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCAcc                                                                  <  1:5002565/93‑1 (MQ=255)
                           tGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCAcc                                                                  >  2:6588735/1‑93 (MQ=255)
                              aGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGaa                                                               >  2:4354353/1‑93 (MQ=255)
                              aGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGaa                                                               >  2:7475199/1‑93 (MQ=255)
                                  tCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAg                                                           <  1:3910200/93‑1 (MQ=255)
                                    aTACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGaa                                                         <  1:4488446/93‑1 (MQ=255)
                                    aTACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGaa                                                         <  1:2311400/93‑1 (MQ=255)
                                      actgactgGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACg                                                       >  1:4392104/1‑93 (MQ=255)
                                      actgactgGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACg                                                       <  2:2255794/93‑1 (MQ=255)
                                      actgactgGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACg                                                       <  2:6941303/93‑1 (MQ=255)
                                       ctgactgGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACgg                                                      <  2:5556839/93‑1 (MQ=255)
                                        tgactgGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGt                                                     <  2:6639239/93‑1 (MQ=255)
                                         gactgGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTa                                                    <  2:4210830/93‑1 (MQ=255)
                                           ctgGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACt                                                  <  2:1468492/93‑1 (MQ=255)
                                                catcatCAGAGGTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCAt                                             <  1:3961996/93‑1 (MQ=255)
                                                 atcatcAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCAtt                                            <  2:872874/93‑1 (MQ=255)
                                                 atcatcAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCAtt                                            >  2:1353201/1‑93 (MQ=255)
                                                 atcatcAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCAtt                                            <  1:3098316/93‑1 (MQ=255)
                                                  tcatcaGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTg                                           >  2:5375247/1‑93 (MQ=255)
                                                  tcatcaGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTg                                           <  1:2008070/93‑1 (MQ=255)
                                                   catcaGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGc                                          <  1:4250979/93‑1 (MQ=255)
                                                   catcaGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGc                                          <  1:5818378/93‑1 (MQ=255)
                                                    atcaGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCa                                         <  1:2675448/93‑1 (MQ=255)
                                                     tcaGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCAt                                        >  1:1356/1‑93 (MQ=255)
                                                     tcaGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCAt                                        >  1:4010771/1‑93 (MQ=255)
                                                      caGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATc                                       <  2:4788618/93‑1 (MQ=255)
                                                       agagTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCa                                      >  1:2060894/1‑93 (MQ=255)
                                                        gagTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCAc                                     >  1:5748880/1‑93 (MQ=255)
                                                        gagTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCAc                                     <  2:593993/93‑1 (MQ=255)
                                                         agTTTGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACg                                    >  1:5162637/1‑93 (MQ=255)
                                                           tttGCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCa                                  <  2:1518275/93‑1 (MQ=255)
                                                              gCACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGa                               <  1:7038975/93‑1 (MQ=255)
                                                               cACTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTTATTGCATCACGCATGaa                              <  2:198549/93‑1 (MQ=255)
                                                                 cTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACt                            >  1:5923653/1‑93 (MQ=255)
                                                                 cTTGAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACt                            <  1:2639465/93‑1 (MQ=255)
                                                                  ttGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTc                           >  1:5237290/1‑93 (MQ=255)
                                                                    gAGCCATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCg                         <  2:1776113/93‑1 (MQ=255)
                                                                       ccATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGa                      <  2:1142331/93‑1 (MQ=255)
                                                                        cATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGAt                     >  1:3730750/1‑93 (MQ=255)
                                                                        cATAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGAt                     <  1:2941734/93‑1 (MQ=255)
                                                                         aTAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATc                    >  1:4963758/1‑93 (MQ=255)
                                                                         aTAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATc                    >  2:1815077/1‑93 (MQ=255)
                                                                          tAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCa                   <  1:5553745/93‑1 (MQ=255)
                                                                          tAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCa                   >  2:3806346/1‑93 (MQ=255)
                                                                          tAAAGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCa                   <  1:615975/93‑1 (MQ=255)
                                                                           aaaGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAg                  <  1:4907271/93‑1 (MQ=255)
                                                                            aaGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAgg                 >  2:651736/1‑93 (MQ=255)
                                                                            aaGTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAgg                 <  2:4539157/93‑1 (MQ=255)
                                                                              gTCCATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCa               >  1:4537644/1‑93 (MQ=255)
                                                                                ccATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGa             <  2:4129040/93‑1 (MQ=255)
                                                                                ccATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGa             <  1:7073998/93‑1 (MQ=255)
                                                                                ccATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGa             <  2:724560/93‑1 (MQ=255)
                                                                                ccATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGa             >  1:6052828/1‑93 (MQ=255)
                                                                                 cATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGAt            >  2:2990262/1‑93 (MQ=255)
                                                                                 cATAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGAt            >  2:2847399/1‑93 (MQ=255)
                                                                                  aTAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATa           <  1:3638198/93‑1 (MQ=255)
                                                                                   tAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATaa          >  2:6317312/1‑93 (MQ=255)
                                                                                   tAATCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATaa          >  2:360166/1‑93 (MQ=255)
                                                                                    aaTCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATaaa         <  2:3231887/93‑1 (MQ=255)
                                                                                    aaTCACGACAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATaaa         <  2:3449881/93‑1 (MQ=255)
                                                                                     aTCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAg        >  1:1674904/1‑93 (MQ=255)
                                                                                      tCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAgg       >  2:2621614/1‑93 (MQ=255)
                                                                                      tCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAgg       <  2:1684543/93‑1 (MQ=255)
                                                                                      tCACGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAgg       <  1:1580998/93‑1 (MQ=255)
                                                                                      tCACGACAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAgg       <  1:3148006/93‑1 (MQ=255)
                                                                                        aCGGCAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTa     <  1:3269436/93‑1 (MQ=255)
                                                                                           |                                                                                            
GTAGCCTTTCAGGTTCGCCTTCTTCAATGCAGACTCATACTGACTGGACATCATCAGAGTTTGCACTTGAGCCATAAAGTCCATAATCACGACAACAACGATAAGCAGTGAGGTCCCACCGAAGTAGAACGGTACTTTCATTGCATCACGCATGAACTCCGGGATCAGGCAGATAAAGGTAATA  >  NC_000913/3440412‑3440595

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