Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,926,900 T→G 21.4% D647A (GAC→GCC)  ptrB ← protease II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,926,9000TG21.4% 331.4 / 98.4 160D647A (GAC→GCC) ptrBprotease II
Reads supporting (aligned to +/- strand):  ref base T (67/58);  new base G (15/19);  total (83/77)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

Only 100 of 160 total aligned reads displayed.

ATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAA  >  NC_000913/1926802‑1926998
                                                                                                  |                                                                                                  
aTGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTc                                                                                                   >  2:3165026/1‑100 (MQ=255)
 tGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCg                                                                                                  <  2:828023/100‑1 (MQ=255)
 tGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCg                                                                                                  >  2:414461/1‑100 (MQ=255)
  gTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                                                 >  1:98876/1‑100 (MQ=255)
  gTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                                                 <  2:3614100/100‑1 (MQ=255)
  gTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                                                 <  2:3597335/100‑1 (MQ=255)
  gTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                                                 >  1:2332432/1‑100 (MQ=255)
   ttCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTCAAGCGACCAGATTTACCGACATGGACTCAGACCATGGCGGt                                                                                                <  1:3080468/100‑1 (MQ=255)
   ttCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGt                                                                                                <  2:1546523/100‑1 (MQ=38)
   ttCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGt                                                                                                <  2:1654051/100‑1 (MQ=38)
    tCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTCCCGCCATGGCCTGAGTCCATGTCGGt                                                                                                >  1:2429337/1‑99 (MQ=255)
    tCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                                                               >  1:1852818/1‑100 (MQ=255)
    tCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGtn                                                                                               <  1:1226116/100‑2 (MQ=37)
     cccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTAc                                                                                              >  2:3397049/1‑100 (MQ=255)
     cccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaa                                                                                              <  1:817959/100‑2 (MQ=37)
     cccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaa                                                                                              <  2:3279350/100‑2 (MQ=37)
      ccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACa                                                                                             <  2:461241/100‑1 (MQ=255)
      ccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaa                                                                                             <  1:1014626/100‑3 (MQ=37)
      ccTCCGTCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCTTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                                                               <  1:1710916/98‑1 (MQ=255)
        tGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaatc                                                                                           <  2:1125792/100‑5 (MQ=25)
         gcgcCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACagag                                                                                          >  2:2034507/1‑100 (MQ=255)
          cgcCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGc                                                                                         <  2:145903/100‑1 (MQ=255)
          cgcCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaatctg                                                                                         <  1:2087724/100‑7 (MQ=25)
           gcCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCAGGTCGGTACAGAGCa                                                                                        >  2:1906250/1‑100 (MQ=255)
           gcCAGCGCGACCAGAAAAGCATATACCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaatctgg                                                                                        <  2:1374484/100‑8 (MQ=25)
            ccAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCaa                                                                                       >  1:532261/1‑100 (MQ=255)
            ccAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCaa                                                                                       >  1:1904497/1‑100 (MQ=255)
             cAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaatctggtc                                                                                      <  1:4269559/100‑10 (MQ=21)
             cAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaatctggtc                                                                                      <  1:1858161/100‑10 (MQ=21)
               gcgcGACCAGAAAAGCATATTCCATCGCTACGCCTTAGTACGAGTTAAAGCAACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaa                                                                                    >  1:2382811/1‑100 (MQ=255)
                 gtgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaaa                                                                                  <  2:4397843/98‑1 (MQ=255)
                 gcgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaaa                                                                                  <  2:3155551/100‑1 (MQ=255)
                  cgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACGGCCATGACCTGAGTCCATGTCGGTACAGAGCAATAAAAg                                                                                 >  1:3090915/1‑100 (MQ=255)
                  cgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAg                                                                                 >  1:3223100/1‑100 (MQ=255)
                   gACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGa                                                                                >  2:2265545/1‑100 (MQ=255)
                      cAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATgg                                                                             <  2:1620539/100‑1 (MQ=255)
                       aGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGt                                                                            <  2:1242734/100‑1 (MQ=255)
                       aGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGt                                                                            >  1:2303774/1‑100 (MQ=255)
                       aGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGt                                                                            >  1:4013278/1‑100 (MQ=255)
                            aGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGtacggc                                                                       >  2:1728919/1‑95 (MQ=255)
                            aGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCg                                                                       >  2:1189190/1‑100 (MQ=255)
                             gCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCgg                                                                      <  1:74069/100‑1 (MQ=255)
                             gCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCgg                                                                      >  2:4366679/1‑100 (MQ=255)
                              cATATTCCATCGCTACGCCTTCGTACGATTTAAGACGACCAGATTTACCGCCATGGCCTGAGTCCATGCCGGTAAAGAACAATAAAAGCCGCTCATCGGt                                                                     >  1:4391296/1‑100 (MQ=255)
                                 atTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACCGAGCAATAAAAGATGGTAATCGGtttt                                                                  >  1:4188980/1‑100 (MQ=255)
                                  ttacATCGCTACGCGTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTc                                                                 <  1:3549581/97‑1 (MQ=255)
                                  ttCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTc                                                                 >  2:375489/1‑100 (MQ=255)
                                   tCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCGATAAAAGATGGTCATCGGTTTTc                                                                 <  1:541652/99‑1 (MQ=255)
                                   tCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCa                                                                <  1:1001494/100‑1 (MQ=255)
                                   tCCATCGCTAAGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGTTACAGAGCAATAAAAGATGGTCATCGGTTTTCa                                                                <  1:2892560/100‑1 (MQ=255)
                                    ccATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAg                                                               >  1:4294148/1‑100 (MQ=255)
                                     cATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGc                                                              >  1:2856067/1‑100 (MQ=255)
                                       tCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTc                                                            <  1:2921580/100‑1 (MQ=255)
                                       tCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTc                                                            <  1:3165008/100‑1 (MQ=255)
                                          cTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTcgcg                                                         >  2:2250451/1‑100 (MQ=255)
                                                 ttCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTa                                                  >  2:4225196/1‑100 (MQ=255)
                                                    gTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCg                                               <  2:2552680/100‑1 (MQ=255)
                                                      ncGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAACGATGGTCATCGGCTTTCAGCTCGCGCAATTTAGCGAc                                             >  1:1226103/2‑100 (MQ=255)
                                                        gATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGAccc                                           >  1:2467751/1‑100 (MQ=255)
                                                         aTTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCa                                          >  1:1418886/1‑100 (MQ=255)
                                                          tttAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAt                                         >  2:4224635/1‑100 (MQ=255)
                                                            tAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAttt                                       >  2:2134325/1‑100 (MQ=255)
                                                               aGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGc                                    <  1:230596/100‑1 (MQ=255)
                                                               aGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGc                                    <  2:1475253/100‑1 (MQ=255)
                                                               aGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGc                                    >  2:4077247/1‑100 (MQ=255)
                                                               aGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCAGCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGc                                    <  1:2362509/100‑1 (MQ=255)
                                                                  gACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg                                 >  2:3530943/1‑100 (MQ=255)
                                                                  gACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg                                 <  2:3877157/100‑1 (MQ=255)
                                                                   aCCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg                                 <  2:2753876/99‑1 (MQ=255)
                                                                   aCCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGt                                >  1:1518401/1‑100 (MQ=255)
                                                                     cAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTc                              <  1:3624687/100‑1 (MQ=255)
                                                                      aGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTcc                             >  2:3796116/1‑100 (MQ=255)
                                                                       tctgtACCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTccc                            >  1:3327067/5‑100 (MQ=25)
                                                                       gATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTccc                            <  2:3068410/100‑1 (MQ=255)
                                                                       gATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTcc                             >  1:1349070/1‑99 (MQ=255)
                                                                         tttACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCaa                          >  2:1054586/1‑100 (MQ=255)
                                                                         tgtACCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCAGGTTcac                            >  2:3018956/3‑96 (MQ=25)
                                                                          gtACCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAat                         >  1:3060838/2‑100 (MQ=37)
                                                                           tACCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAata                        >  2:3038106/1‑100 (MQ=38)
                                                                            aCCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAatat                       >  2:545684/1‑100 (MQ=38)
                                                                            aCCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAatat                       >  2:3496408/1‑100 (MQ=38)
                                                                            aCCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAatat                       >  2:1695112/1‑100 (MQ=38)
                                                                             ccGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATAtt                      <  1:1650476/100‑1 (MQ=255)
                                                                             ccGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATAtt                      >  1:1871828/1‑100 (MQ=38)
                                                                              cGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTg                     <  2:4318941/100‑1 (MQ=255)
                                                                              cGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTg                     <  1:1853513/100‑1 (MQ=255)
                                                                               gCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGc                    <  1:4206091/100‑1 (MQ=255)
                                                                               gCCATGGCCTGAGTCCACGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGc                    <  2:3082661/100‑1 (MQ=255)
                                                                                 cATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCan                  <  2:4004494/100‑2 (MQ=255)
                                                                                 cATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCAc                  <  1:399087/100‑1 (MQ=255)
                                                                                  aTGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACTCATTTTGCCGGTTCCCAATATTGCAcc                 <  1:690010/100‑1 (MQ=255)
                                                                                   tGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCt                <  2:1309259/100‑1 (MQ=255)
                                                                                    ggCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTg               <  1:4224604/100‑1 (MQ=255)
                                                                                     gCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTga              >  2:2794232/1‑100 (MQ=255)
                                                                                      ccTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTgag             <  1:2034500/100‑1 (MQ=255)
                                                                                      ccTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTgag             >  1:2694967/1‑100 (MQ=255)
                                                                                      actcaggccATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTgag             <  1:1670423/93‑1 (MQ=255)
                                                                                           gTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCg        >  2:4266443/1‑100 (MQ=255)
                                                                                             ccATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTg      >  2:4245120/1‑100 (MQ=255)
                                                                                              cATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTGTGCCGGTGCCCAATATTGCACCTGAGAAGCGTg      <  1:2594425/99‑1 (MQ=255)
                                                                                                  |                                                                                                  
ATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAA  >  NC_000913/1926802‑1926998

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