Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,751,328 A→G M176T (ATG→ACG)  b1673 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,751,3280AG98.5% 429.5 / ‑2.6 137M176T (ATG→ACG) b1673hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (72/63);  total (73/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

Only 100 of 141 total aligned reads displayed.

TCGTCCCTTCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACATGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTGGTCAATCGACAAATTTCTTCCGTC  >  NC_000913/1751239‑1751417
                                                                                         |                                                                                         
tCGTCCCTTCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGc                                                                                        >  2:2981006/1‑93 (MQ=255)
 cGTCCCTTCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGcc                                                                                       <  2:3422284/93‑1 (MQ=255)
 cGTCCCTTCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGcc                                                                                       <  1:1577093/93‑1 (MQ=255)
 cGTCCCTTCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGcc                                                                                       <  1:3945038/93‑1 (MQ=255)
 cGTCCCTCTAAACGCAATCGCTTGCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGcc                                                                                       <  1:3945039/93‑1 (MQ=255)
  gTCCCTTCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCa                                                                                      >  1:2112330/1‑93 (MQ=255)
   tCCCTTCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCag                                                                                     >  2:166169/1‑93 (MQ=255)
   tCCCTTCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCag                                                                                     >  1:3771703/1‑93 (MQ=255)
    cccTTCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCaga                                                                                    >  2:3185514/1‑93 (MQ=255)
    cccTTCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCaga                                                                                    <  2:3945070/93‑1 (MQ=255)
    cccTTCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCaga                                                                                    >  2:3185513/1‑93 (MQ=255)
     ccTTCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCagag                                                                                   >  2:3871910/1‑93 (MQ=255)
     ccTTCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCagag                                                                                   <  1:3397711/93‑1 (MQ=255)
      cTTCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCagaga                                                                                  <  1:579574/93‑1 (MQ=255)
      cTTCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCagaga                                                                                  <  2:1106250/93‑1 (MQ=255)
        tCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCagagaga                                                                                >  2:2927160/1‑93 (MQ=255)
        tCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCagagaga                                                                                <  1:881161/93‑1 (MQ=255)
        tCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCagagaga                                                                                <  1:2590037/93‑1 (MQ=255)
         cAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCagagagag                                                                               <  2:2931237/93‑1 (MQ=255)
          aaCCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAgg                                                                              <  2:3066878/93‑1 (MQ=255)
           aCCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCTCACCGCCGCTGTGGTTACGGCTATTCAACATTCCATAGAGAGGa                                                                             >  2:3775317/1‑93 (MQ=255)
           aCCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGa                                                                             >  2:3775467/1‑93 (MQ=255)
           aCCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGa                                                                             >  2:41251/1‑93 (MQ=255)
            ccGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGaa                                                                            <  1:2402409/93‑1 (MQ=255)
            ccGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGaa                                                                            <  1:2406406/93‑1 (MQ=255)
             cGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAAc                                                                           >  1:1749649/1‑93 (MQ=255)
                aaTCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAg                                                                        >  2:2526581/1‑93 (MQ=255)
                 aTCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAgg                                                                       >  2:139429/1‑93 (MQ=255)
                 aTCGCTTTCAGGTTTTTCGAACCCATTAGTGCGCCAGTTCCCGCCCCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAgg                                                                       <  2:2640586/93‑1 (MQ=255)
                  tCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGt                                                                      <  2:43246/93‑1 (MQ=255)
                  tCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGt                                                                      >  2:368440/1‑93 (MQ=255)
                  tCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGt                                                                      >  1:3315585/1‑93 (MQ=255)
                   tgCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGtt                                                                     <  2:4087721/92‑1 (MQ=255)
                    gCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGttt                                                                    >  1:433173/1‑93 (MQ=255)
                     cTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGtttt                                                                   >  1:512276/1‑93 (MQ=255)
                      tttCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTc                                                                  >  2:3883751/1‑93 (MQ=255)
                        tCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTccc                                                                <  2:4074910/93‑1 (MQ=255)
                        tCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTccc                                                                >  2:3358466/1‑93 (MQ=255)
                        tCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTccc                                                                >  2:3411570/1‑93 (MQ=255)
                         cAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTcccc                                                               >  1:4728087/1‑93 (MQ=255)
                           ggTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAg                                                             >  1:1568658/1‑93 (MQ=255)
                           ggTTTTCCGACCCCATTATGGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAg                                                             <  2:1281010/93‑1 (MQ=255)
                             tttttCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGcc                                                           <  2:4627680/93‑1 (MQ=255)
                              ttttCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGCGAGAGGAACAAGGTTTTCACCAGCCt                                                          >  2:2538721/1‑93 (MQ=255)
                              ttttCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCt                                                          <  1:3452941/93‑1 (MQ=255)
                               tttCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTg                                                         >  1:4179700/1‑93 (MQ=255)
                                    aaCCCATTATGGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACATGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCaa                                                    <  1:2348003/93‑1 (MQ=255)
                                     aCCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAAt                                                   <  1:2263967/93‑1 (MQ=255)
                                     aCCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAAt                                                   >  1:168821/1‑93 (MQ=255)
                                      cccATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATa                                                  >  2:773867/1‑93 (MQ=255)
                                      cccATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATa                                                  <  1:676302/93‑1 (MQ=255)
                                       ccATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAg                                                 >  1:198646/1‑93 (MQ=255)
                                        cATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAgc                                                <  2:1501193/93‑1 (MQ=255)
                                        cATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAgc                                                <  1:618192/93‑1 (MQ=255)
                                         attattGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAgcc                                               >  1:3362625/1‑93 (MQ=255)
                                         attattGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAgcc                                               >  2:3181159/1‑93 (MQ=255)
                                         attattGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAgcc                                               >  2:4259223/1‑93 (MQ=255)
                                          ttattGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTACCCAGCCTGACCAATAgccg                                              <  2:3339849/93‑1 (MQ=255)
                                                 gcCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCcacaca                                       >  1:2991519/1‑93 (MQ=255)
                                                 gcCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCcacaca                                       >  1:2991520/1‑93 (MQ=255)
                                                 gcCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCcacaca                                       >  2:1740841/1‑93 (MQ=255)
                                                 gcCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCcacaca                                       <  1:1113172/93‑1 (MQ=255)
                                                  ccAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAg                                      <  2:2222009/93‑1 (MQ=255)
                                                   cAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAgg                                     >  1:4656207/1‑93 (MQ=255)
                                                    aGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGt                                    >  1:908358/1‑93 (MQ=255)
                                                      ttCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACCAGGTTTTCCCCAGCCTGACCAATAGCCGCCACCCAggatc                                  >  2:3954227/1‑90 (MQ=255)
                                                       tCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGCGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTc                                 <  2:2696314/93‑1 (MQ=255)
                                                       tCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTc                                 <  2:2671307/93‑1 (MQ=255)
                                                       tCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTc                                 <  2:2968074/93‑1 (MQ=255)
                                                       tCCCGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTc                                 <  1:1458400/93‑1 (MQ=255)
                                                        cccGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTcc                                <  1:3626810/93‑1 (MQ=255)
                                                         ccGCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCg                               <  2:3359779/93‑1 (MQ=255)
                                                           gCACCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGa                             <  1:1136663/93‑1 (MQ=255)
                                                             aCCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACt                           >  2:2878559/1‑93 (MQ=255)
                                                             aCCGCCGCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACt                           >  2:1000089/1‑93 (MQ=255)
                                                              ccgccgCTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACta                          <  2:2776291/93‑2 (MQ=255)
                                                               cgccgcTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACtag                         >  1:1836007/1‑91 (MQ=255)
                                                                   gcTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCaa                     >  2:4442770/1‑93 (MQ=255)
                                                                   gcTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCaa                     >  2:4442769/1‑93 (MQ=255)
                                                                   gcTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCaa                     >  2:1958349/1‑93 (MQ=255)
                                                                   gcTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCaa                     <  1:2018498/93‑1 (MQ=255)
                                                                    cTGTGGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAAt                    <  2:1664491/93‑1 (MQ=255)
                                                                       tgGTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGa                 >  1:3916965/1‑93 (MQ=255)
                                                                         gTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACa               >  2:4602950/1‑93 (MQ=255)
                                                                         gTTACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACa               >  2:445319/1‑93 (MQ=255)
                                                                          ttACGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACaa              >  1:1640585/1‑93 (MQ=255)
                                                                            aCGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAAt            >  2:1065376/1‑93 (MQ=255)
                                                                             cGGCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAAtt           <  1:2843431/93‑1 (MQ=255)
                                                                              ggCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAAttt          <  1:1988658/93‑1 (MQ=255)
                                                                              ggCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAAttt          <  2:3903359/93‑1 (MQ=255)
                                                                              ggCTATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAAttt          <  2:1138291/93‑1 (MQ=255)
                                                                                cTATTCAACGTGCGAGTGAGATGAACAAGGTTTGCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATttct        <  2:2276297/93‑1 (MQ=255)
                                                                                 tATTCAACGTGCCAGCGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGACGCCACACAGGTTTCCGGACTGGTCAATCGACAAATttctc       >  2:2862983/1‑92 (MQ=255)
                                                                                 tATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATttctt       >  2:2641646/1‑93 (MQ=255)
                                                                                 tATTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATttctt       >  2:3652139/1‑93 (MQ=255)
                                                                                  aTTCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATttcttc      <  1:1711035/93‑1 (MQ=255)
                                                                                   ttCAACGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTcc     >  2:549380/1‑93 (MQ=255)
                                                                                    tCAACGTGCCAGAGAGAGGAACAAGGTCTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCcg    >  1:390612/1‑93 (MQ=255)
                                                                                      aaCGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCcgtc  >  1:4395403/1‑93 (MQ=255)
                                                                                      aaCGTGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTAGTCAATCGACAAATTTCTTCcgtc  >  2:2311528/1‑93 (MQ=255)
                                                                                         |                                                                                         
TCGTCCCTTCAACCGCAATCGCTTTCAGGTTTTTCGAACCCATTATTGCGCCAGTTCCCGCACCGCCGCTGTGGTTACGGCTATTCAACATGCCAGAGAGAGGAACAAGGTTTTCCCCAGCCTGACCAATAGCCGCCACACAGGTTTCCGGACTGGTCAATCGACAAATTTCTTCCGTC  >  NC_000913/1751239‑1751417

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