Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,310,422 T→C D133G (GAC→GGC)  rbfA ← ribosome‑binding factor A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,310,4220TC99.6% 944.6 / ‑4.7 277D133G (GAC→GGC) rbfAribosome‑binding factor A
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (157/117);  total (159/118)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.90e-01

Only 100 of 293 total aligned reads displayed.

AGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACC  >  NC_000913/3310333‑3310511
                                                                                         |                                                                                         
aGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGccct                                                                                        >  1:2599054/1‑93 (MQ=255)
aGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGccct                                                                                        >  1:5800700/1‑93 (MQ=255)
  ggcATCCTTTCTGGACATTCCCTTAGGTTTAATCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGccctcc                                                                                      <  1:646171/92‑1 (MQ=255)
   gcATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCcctcct                                                                                     >  2:2268920/1‑93 (MQ=255)
   gcATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCcctcct                                                                                     >  2:6220727/1‑93 (MQ=255)
   gcATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCcctcct                                                                                     <  2:1159326/93‑1 (MQ=255)
   gcATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCcctcct                                                                                     <  2:1014596/93‑1 (MQ=255)
   gcATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCcctcct                                                                                     <  1:4333604/93‑1 (MQ=255)
    cATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCtt                                                                                    >  1:4943924/1‑93 (MQ=255)
    cATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCtt                                                                                    >  2:1512099/1‑93 (MQ=255)
    cATCGTTGCTGGACATACCCTGAGGTTTATCCAACAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCtt                                                                                    >  1:3494413/1‑93 (MQ=255)
     ctCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTg                                                                                   <  1:2438282/92‑1 (MQ=255)
     aTCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTg                                                                                   >  2:5792724/1‑93 (MQ=255)
      tCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGc                                                                                  >  2:2379682/1‑93 (MQ=255)
          tGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTgtc                                                                              <  1:4979128/93‑1 (MQ=255)
           gCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTgtcg                                                                             >  1:1592158/1‑93 (MQ=255)
             tgcaCATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTgtcgtc                                                                           <  2:6353985/90‑1 (MQ=255)
             tGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTgtcgtc                                                                           <  2:1190484/93‑1 (MQ=255)
              ggACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTcc                                                                          <  2:73416/93‑1 (MQ=255)
              ggACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTACTTGCTGTCGTcc                                                                          >  2:5111854/1‑93 (MQ=255)
                aCATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCgg                                                                        >  2:2210918/1‑93 (MQ=255)
                aCATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCgg                                                                        >  1:4772353/1‑93 (MQ=255)
                  aTACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGCGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGt                                                                      >  2:616535/1‑93 (MQ=255)
                    aCCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTa                                                                    >  2:3055137/1‑93 (MQ=255)
                    aCCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCATTGCTGTCGTCCGGGTTa                                                                    >  2:3416410/1‑93 (MQ=255)
                     cccTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTaa                                                                   >  1:3394454/1‑93 (MQ=255)
                     cccTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTaa                                                                   >  1:5691853/1‑93 (MQ=255)
                      ccTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAac                                                                  <  2:3616154/93‑1 (MQ=255)
                       cTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAaca                                                                 <  1:3334665/93‑1 (MQ=255)
                       cTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAaca                                                                 >  2:461957/1‑93 (MQ=255)
                         gAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCTGGTTAACacg                                                               <  1:3964274/93‑1 (MQ=255)
                         gAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACacg                                                               >  2:5668737/1‑93 (MQ=255)
                         gAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACacg                                                               >  1:5484085/1‑93 (MQ=255)
                          aGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACacga                                                              <  1:94749/93‑1 (MQ=255)
                           ggTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACacgac                                                             <  2:4452125/93‑1 (MQ=255)
                           ggTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACacgac                                                             <  2:1637127/93‑1 (MQ=255)
                               tATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGttc                                                         >  1:6344655/1‑93 (MQ=255)
                                agccAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGttct                                                        <  2:5379536/91‑1 (MQ=255)
                                aTCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGttct                                                        >  2:4661308/1‑93 (MQ=255)
                                aTCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGttct                                                        >  2:5191341/1‑93 (MQ=255)
                                 tCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGttctt                                                       <  2:1853894/93‑1 (MQ=255)
                                   cAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCg                                                     <  1:5565784/93‑1 (MQ=255)
                                    aGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGt                                                    >  1:4662688/1‑93 (MQ=255)
                                      caacaaAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCa                                                  <  2:4807528/93‑1 (MQ=255)
                                      caacaaAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCa                                                  <  1:5322984/93‑1 (MQ=255)
                                       aacaaAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCAt                                                 <  2:3562709/93‑1 (MQ=255)
                                        acaaAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATg                                                <  1:751014/93‑1 (MQ=255)
                                          aaaaCGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGtt                                              >  1:2834771/1‑93 (MQ=255)
                                          aaaaCGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGtt                                              <  1:5097399/93‑1 (MQ=255)
                                          aaaaCGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGtt                                              >  2:3708495/1‑93 (MQ=255)
                                          aaaaCGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGtt                                              <  1:4029585/93‑1 (MQ=255)
                                           aaaCGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGttt                                             <  2:1768581/93‑1 (MQ=255)
                                           aaaCGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGttt                                             >  2:1583950/1‑93 (MQ=255)
                                             aCGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGa                                           >  2:2564308/1‑93 (MQ=255)
                                               gCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGAcc                                         >  1:4286861/1‑93 (MQ=255)
                                                 cGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCAc                                       <  2:2664785/93‑1 (MQ=255)
                                                   ttAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGc                                     >  2:4919989/1‑93 (MQ=255)
                                                   ttAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCCTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGc                                     <  1:5245918/93‑1 (MQ=255)
                                                     aaTGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTg                                   <  1:5199231/93‑1 (MQ=255)
                                                      aTGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTgg                                  <  1:142339/93‑1 (MQ=255)
                                                       tGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGt                                 <  2:4242414/93‑1 (MQ=255)
                                                        gTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTc                                <  1:1859519/93‑1 (MQ=255)
                                                          cgcgaccgcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCAc                              <  2:2855516/93‑1 (MQ=255)
                                                           gcgaccgcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCAcc                             <  1:4844318/93‑1 (MQ=255)
                                                           gcgaccgcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCAcc                             <  2:5863027/93‑1 (MQ=255)
                                                             caccgcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAg                           <  2:2083345/92‑1 (MQ=255)
                                                              accgcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAgg                          >  1:2039544/1‑93 (MQ=255)
                                                              accgcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAgg                          <  2:628136/93‑1 (MQ=255)
                                                               ccgcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGt                         <  1:4576125/93‑1 (MQ=255)
                                                               ccgcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGt                         <  1:4071773/93‑1 (MQ=255)
                                                                cgcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGtt                        <  2:6061306/93‑1 (MQ=255)
                                                                 gcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGttt                       <  2:4553441/93‑1 (MQ=255)
                                                                 gcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGttt                       >  1:1051836/1‑93 (MQ=255)
                                                                 gcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGttt                       >  1:2702023/1‑93 (MQ=255)
                                                                 gcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGttt                       >  1:1004758/1‑93 (MQ=255)
                                                                 gcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGttt                       <  1:3725946/93‑1 (MQ=255)
                                                                    acgacgaGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGAc                    >  2:1321485/1‑93 (MQ=255)
                                                                     cgacgaGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACa                   >  2:4096339/1‑93 (MQ=255)
                                                                     cgacgaGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACa                   <  1:2859893/93‑1 (MQ=255)
                                                                      gacgaGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACAt                  >  1:6253436/1‑93 (MQ=255)
                                                                       acgaGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATg                 >  2:4004743/1‑93 (MQ=255)
                                                                       acgaGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATg                 >  2:2895935/1‑93 (MQ=255)
                                                                         gaGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATgcg               >  1:1964540/1‑93 (MQ=255)
                                                                         gaGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATgcg               >  1:123026/1‑93 (MQ=255)
                                                                          aGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATgcgc              >  2:5620888/1‑93 (MQ=255)
                                                                           ggACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATgcgcc             >  1:5918045/1‑92 (MQ=255)
                                                                           ggACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCa             >  1:3021436/1‑93 (MQ=255)
                                                                           ggACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCa             >  2:2481176/1‑93 (MQ=255)
                                                                              cgacTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATcc          >  1:1686083/1‑93 (MQ=255)
                                                                              cgacTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATcc          <  1:1320016/93‑1 (MQ=255)
                                                                              cgacTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATcc          <  1:4522823/93‑1 (MQ=255)
                                                                                acTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCt        <  1:2352039/93‑1 (MQ=255)
                                                                                acTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCt        <  2:5974830/93‑1 (MQ=255)
                                                                                   cATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCa     >  2:2784622/1‑93 (MQ=255)
                                                                                    cttAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCaa    <  1:3584566/92‑1 (MQ=255)
                                                                                    aTTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCaa    >  2:4331751/1‑93 (MQ=255)
                                                                                     ttAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAAc   >  2:3788503/1‑93 (MQ=255)
                                                                                     ttAGCCCTCCTTGCTGTCGACCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAAc   >  1:6238705/1‑93 (MQ=255)
                                                                                      tAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAAcc  <  1:3152589/93‑1 (MQ=255)
                                                                                      tAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAAcc  <  1:5128623/93‑1 (MQ=255)
                                                                                         |                                                                                         
AGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACC  >  NC_000913/3310333‑3310511

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