Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,647,924 G→A P795L (CCG→CTG)  b2520 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,647,9240GA85.4% 653.1 / 68.6 277P795L (CCG→CTG) b2520hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/17);  new base A (128/100);  total (150/127)
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.70e-01

Only 100 of 281 total aligned reads displayed.

TTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCT  >  NC_000913/2647832‑2648016
                                                                                            |                                                                                            
   aaTACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGcc                                                                                           <  1:2795731/93‑1 (MQ=255)
     tACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCt                                                                                         <  2:596111/93‑1 (MQ=255)
     tACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCt                                                                                         >  2:5667577/1‑93 (MQ=255)
     tACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCt                                                                                         >  2:2667311/1‑93 (MQ=255)
      aCAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCtt                                                                                        <  1:4793559/93‑1 (MQ=255)
           tCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTg                                                                                   <  2:5341254/93‑1 (MQ=255)
             aTGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGa                                                                                 >  2:1645876/1‑93 (MQ=255)
             aTGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGa                                                                                 >  2:2387104/1‑93 (MQ=255)
              tGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGAc                                                                                >  2:6711950/1‑93 (MQ=255)
               gACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACt                                                                               >  1:3030906/1‑93 (MQ=255)
                aCGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTc                                                                              <  2:2419552/93‑1 (MQ=255)
                  gATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGa                                                                            >  2:5417658/1‑93 (MQ=255)
                   aTTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGAc                                                                           >  1:3257931/1‑93 (MQ=255)
                   aTTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGAc                                                                           >  1:629928/1‑93 (MQ=255)
                   aTTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGAc                                                                           >  1:1794223/1‑93 (MQ=255)
                   aTTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGAc                                                                           <  2:6794348/93‑1 (MQ=255)
                   aTTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGAc                                                                           <  1:3794981/93‑1 (MQ=255)
                     tCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGAcaa                                                                         >  2:904035/1‑91 (MQ=255)
                      ccAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCAt                                                                        <  1:6396019/93‑1 (MQ=255)
                       cAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATc                                                                       >  2:181409/1‑93 (MQ=255)
                       cAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATc                                                                       <  1:4600530/93‑1 (MQ=255)
                       cAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATc                                                                       <  1:5263824/93‑1 (MQ=255)
                       cAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATc                                                                       <  2:4621336/93‑1 (MQ=255)
                         ggTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGc                                                                     <  2:1030029/93‑1 (MQ=255)
                               aTCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAAcc                                                               <  1:3199322/93‑1 (MQ=255)
                                tCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCt                                                              <  1:7151253/93‑1 (MQ=255)
                                  gACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCttt                                                            >  2:7077988/1‑93 (MQ=255)
                                  gACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCttt                                                            >  1:606558/1‑93 (MQ=255)
                                  gACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCttt                                                            >  1:6543846/1‑93 (MQ=255)
                                  gACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCttt                                                            <  2:2720874/93‑1 (MQ=255)
                                   aCCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTa                                                           >  2:6080153/1‑93 (MQ=255)
                                     cGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTAcc                                                         <  2:7279203/93‑1 (MQ=255)
                                     cGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTAcc                                                         >  2:4910635/1‑93 (MQ=255)
                                      ggAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTAccc                                                        >  2:5793485/1‑93 (MQ=255)
                                      ggAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTAccc                                                        <  2:3470246/93‑1 (MQ=255)
                                       gAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACccg                                                       <  2:6321805/93‑1 (MQ=255)
                                          tCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACccgccg                                                    >  1:6846458/1‑93 (MQ=255)
                                          tCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACccgccg                                                    >  2:7464764/1‑93 (MQ=255)
                                           cGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGt                                                   >  2:3164887/1‑93 (MQ=255)
                                           cGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGt                                                   >  2:839099/1‑93 (MQ=255)
                                            gTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGtt                                                  >  1:5414003/1‑93 (MQ=255)
                                            gTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGtt                                                  >  1:6920458/1‑93 (MQ=255)
                                            gTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGtt                                                  >  1:5061521/1‑93 (MQ=255)
                                            gTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCACCGCATACCCTTTCCCCGCCGtt                                                  <  2:763373/93‑1 (MQ=255)
                                             tCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTg                                                 >  1:4381702/1‑93 (MQ=255)
                                               aGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGgc                                               >  1:513323/1‑93 (MQ=255)
                                               aGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGgc                                               >  2:791784/1‑93 (MQ=255)
                                               aGACCCACCCCTTGAGCCCGAACATCAATCTCTTCCCCCCACAGCAGCCCTTCACTGGACTCGACCACCGCATACCCTTTCCCCGCCGTTGgc                                               <  1:1166766/93‑1 (MQ=255)
                                                gATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGgcg                                              >  1:6355237/1‑93 (MQ=255)
                                                 aTCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGgcgc                                             >  1:729348/1‑93 (MQ=255)
                                                  tCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGcgcg                                            >  2:144256/1‑93 (MQ=255)
                                                   ccAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCgcgg                                           >  2:3377873/1‑93 (MQ=255)
                                                   ccAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCgcgg                                           >  1:4118747/1‑93 (MQ=255)
                                                    cAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCgcggc                                          >  2:3129618/1‑93 (MQ=255)
                                                     aGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCgcggcg                                         >  2:4776225/1‑93 (MQ=255)
                                                     aGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCgcggcg                                         >  1:1101525/1‑93 (MQ=255)
                                                      gCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATACCCTTTCCCCGCCGTTGGCGCGGCGa                                        <  2:7191722/93‑1 (MQ=255)
                                                      ccccTTGACCCCGAACATAAACCTCTTCCCTCCGCAGCAGCCCTTCCCTGGCCTCGACCATCGCATAACCTTTCCCCGCCGTTGGCGCGGCGa                                        <  2:4980013/92‑1 (MQ=255)
                                                       cccTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGat                                       >  1:1515027/1‑93 (MQ=255)
                                                       cccTTGAGCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGat                                       >  1:4555640/1‑93 (MQ=255)
                                                       cccTGGACCCCAAACATCAATCTCTTCCAACCACAGCACCCCTTCACGGGACTCGCCCATCGCATACCCTTTCCCCGCCGTTGGCGCGGCGat                                       <  1:4322938/93‑1 (MQ=255)
                                                        ccTTGACCCCAAACATCAATCTCTTCCCCCCACAGCACCCCTTCCCGGGCCCCCCCCATCGCAAACCCTTCCCCCGCCGTTGGCGCGGCGata                                      <  1:3419931/93‑1 (MQ=255)
                                                         ctggAGCCCGAACATCAATCTCTCGCCACCACAGCACCCCTTCCCGGGACTCGACCAGCGCATACCCTTTCCCCGCCGTTGGCGCGGCGatat                                     <  1:5979032/90‑1 (MQ=255)
                                                              gCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAAc                                >  1:3575596/1‑93 (MQ=255)
                                                              gCCCGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAAc                                >  2:4399786/1‑93 (MQ=255)
                                                              gCCCGAACATCAACCTCTTGCCACCACAGCCGCCCTTCACTGACCTCGCCCATCGCATAACCTTTCCCCGCCGTTGGCGCGGCGATATGCAAc                                <  2:2074031/93‑1 (MQ=255)
                                                               cccGAACATCAATCTCTTCCCCCCACGGCGCCCCTTCACTGGCCTCGACCATCGCATACCTTTTCCCCGCCGTTGGCGCGGCGATATGCAACt                               <  1:2206267/93‑1 (MQ=255)
                                                                 cGAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTa                             >  1:4876203/1‑93 (MQ=255)
                                                                 cGAACATCAATCCCTTGCCACCACAGCAGCCCTTCCCTGGACCCGACCACCGCATAACTTTTTCCCGCCGTTGGCGCGGCGATATGCAACTTa                             <  2:3639854/93‑1 (MQ=255)
                                                                  gAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTaa                            <  2:6585984/93‑1 (MQ=255)
                                                                  gAACATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTaa                            >  1:5565891/1‑93 (MQ=255)
                                                                   aaCATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGCCCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAAt                           <  2:2923706/93‑1 (MQ=255)
                                                                   aaCATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATACCCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAAt                           <  2:5638142/93‑1 (MQ=255)
                                                                    acacCATTCTCTTCCCACCACAGCCCCCCTTCACTGGCCTCGACCATCGCATACCTTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATg                          <  1:5870560/93‑1 (MQ=255)
                                                                    aCATCAATCTCTTCCCACCACAGCAGCCCTTCCCGGGCCTCGCCCATCGCATACCCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATg                          <  1:6689163/93‑1 (MQ=255)
                                                                     cATCAATCTCTTGCCACCACAGCAGCCCTTCCCTGGACTCGCCCATCGCATACCCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATgg                         <  2:6941988/93‑1 (MQ=255)
                                                                     cATCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATgg                         >  1:3151628/1‑93 (MQ=255)
                                                                      aTCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGgt                        >  2:1165178/1‑93 (MQ=255)
                                                                      aTCAATCTCTCCCCACCACGGCACCCCTCCCCGGGACCCGACCATCCCATACCCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGgt                        <  1:4453198/93‑1 (MQ=255)
                                                                       tCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGgtg                       >  2:2985441/1‑93 (MQ=255)
                                                                       tCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGgtg                       >  2:4065137/1‑93 (MQ=255)
                                                                       tCAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGgtg                       >  1:230422/1‑93 (MQ=255)
                                                                        cAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGgtgt                      >  1:3054544/1‑93 (MQ=255)
                                                                        cAATCTCTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATACCCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGgtgt                      <  2:3877767/93‑1 (MQ=255)
                                                                        cAATCTCTTGCCACCACAGCACCCCTTCACTGGCCTCGACCATCGCATACCCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGgtgt                      <  2:3294221/93‑1 (MQ=255)
                                                                         aaTCTCTTGCCCCCACAGCAGCCCTTCACTGGACTCGCCCATCGCATACCCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTc                     <  2:5528849/93‑1 (MQ=255)
                                                                         aaTCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTc                     >  1:6956306/1‑93 (MQ=255)
                                                                           tctcTTGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGc                   >  1:5299277/1‑93 (MQ=255)
                                                                              ctggCCACCACAGCAGCCCTTCCCTGGACTCGACCATCGCATACCCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAg                <  2:4695091/90‑1 (MQ=255)
                                                                              ctTCCCACCCCAGCCCCCCTTCACTGGCCTCGACCACCGCATACCCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAg                <  1:1724886/93‑1 (MQ=255)
                                                                               ttGCCACCACAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAgg               >  2:3906888/1‑93 (MQ=255)
                                                                               tcgCCCCCACAGCACCCCTCCACGGGACTCGCCCATCGCATACCTTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAgg               <  2:1223834/91‑1 (MQ=255)
                                                                                   caccacAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGa           >  2:6117366/1‑93 (MQ=255)
                                                                                    cccacAGCCCCCCTCCCCTGGCCTCCCCCACCGCATACCCTTCCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGAt          <  1:1289376/92‑1 (MQ=255)
                                                                                    cccacAGCACCCCTTCACGGGACCCGCCCATCGCATACCCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGAt          <  1:1912781/92‑1 (MQ=255)
                                                                                     ccacAGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATa         >  2:6054721/1‑93 (MQ=255)
                                                                                      cacaGCCCCCCTTCACGGGCCTCGCCCACCGCATACCCTTCCCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATaa        <  1:6210288/93‑1 (MQ=255)
                                                                                      cacaGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATaa        >  2:3727841/1‑93 (MQ=255)
                                                                                       acaGCAGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAAc       >  2:2458442/1‑93 (MQ=255)
                                                                                            ggCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCt  >  1:7526557/1‑93 (MQ=255)
                                                                                            |                                                                                            
TTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCT  >  NC_000913/2647832‑2648016

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