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.7% 1170.2 / ‑5.4 341D133G (GAC→GGC) rbfAribosome‑binding factor A
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (173/167);  total (173/168)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 380 total aligned reads displayed.

CAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACC  >  NC_000913/3310332‑3310511
                                                                                          |                                                                                         
cAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGccc                                                                                         <  2:5925054/93‑1 (MQ=255)
 aGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCATCAACAAAACGCCGTTAATGTCGCGACAGCGACGACGAGGACGACTCATTAGccct                                                                                        >  1:1492098/1‑93 (MQ=255)
 aGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGccct                                                                                        >  1:4615098/1‑93 (MQ=255)
 aGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGccct                                                                                        >  1:8426689/1‑93 (MQ=255)
  gcgcATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGccctc                                                                                       >  2:906354/1‑93 (MQ=255)
   cgcATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGccctcc                                                                                      <  1:4779663/93‑1 (MQ=255)
    gcATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCcctcct                                                                                     <  1:224058/93‑1 (MQ=255)
    gcATCGATGCTGGACATAACCTGAGGCTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGAATCATTAGCcctcct                                                                                     >  2:7884482/1‑93 (MQ=255)
    ccATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCcctcct                                                                                     <  1:3259401/92‑1 (MQ=255)
     cATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCtt                                                                                    >  1:2385870/1‑93 (MQ=255)
      aTCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTg                                                                                   <  1:1122334/93‑1 (MQ=255)
       tCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGc                                                                                  >  2:3051014/1‑93 (MQ=255)
       tCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGc                                                                                  <  1:2514551/93‑1 (MQ=255)
       tCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGc                                                                                  >  1:6963952/1‑93 (MQ=255)
        cGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCt                                                                                 <  2:8646459/93‑1 (MQ=255)
         gTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTg                                                                                <  1:3098402/93‑1 (MQ=255)
          ttGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTgt                                                                               <  1:3150802/93‑1 (MQ=255)
             cTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTgtcgt                                                                            >  2:2546437/1‑93 (MQ=255)
              tGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTgtcgtc                                                                           <  1:915817/93‑1 (MQ=255)
              tGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTgtcgtc                                                                           >  1:8198916/1‑93 (MQ=255)
               ggACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTcc                                                                          <  1:1833172/93‑1 (MQ=255)
               ggACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTcc                                                                          <  1:4942809/93‑1 (MQ=255)
                 aCATACCCTGAGGTTTAGCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCgg                                                                        <  1:2794469/93‑1 (MQ=255)
                   aTACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATTGCGCGACCTCGACGACGAGGACGACTCATTAGCACTCCTTGCTGTCGTCCGGGt                                                                      >  2:1248426/1‑93 (MQ=255)
                       ccTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAac                                                                  >  2:5434531/1‑93 (MQ=255)
                       ccTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAac                                                                  >  2:528245/1‑93 (MQ=255)
                         tGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAacac                                                                >  2:8457010/1‑93 (MQ=255)
                         tGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAacac                                                                >  1:7527885/1‑93 (MQ=255)
                          gAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACacg                                                               >  1:7125906/1‑93 (MQ=255)
                           aGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACacga                                                              >  2:2079225/1‑93 (MQ=255)
                            ggTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACacgac                                                             >  1:8717240/1‑93 (MQ=255)
                            ggTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACacgac                                                             >  2:4241753/1‑93 (MQ=255)
                             gTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACacgacg                                                            >  1:6072602/1‑93 (MQ=255)
                             gTTTATATAGCAACAAAACGCCGTTAATGTCGGGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCCGTCGTCCGGGGTAACatgaag                                                            >  2:8552169/1‑91 (MQ=255)
                                tATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGttc                                                         <  2:3659134/93‑1 (MQ=255)
                                tATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGttc                                                         >  1:8492025/1‑93 (MQ=255)
                                tATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGttc                                                         <  1:2079219/93‑1 (MQ=255)
                                tATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGttc                                                         <  2:6409210/93‑1 (MQ=255)
                                 aTCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGttct                                                        <  2:130818/93‑1 (MQ=255)
                                 aTCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGttct                                                        >  1:8602266/1‑93 (MQ=255)
                                  tCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGttctt                                                       <  2:3992468/93‑1 (MQ=255)
                                  tCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGttctt                                                       <  2:8343861/93‑1 (MQ=255)
                                  tCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGttctt                                                       <  1:6014375/93‑1 (MQ=255)
                                  accAGCAACAAAACGCGCATAATGCCGAGACCGAGACCACGAGGCCGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGCCGttctt                                                       <  2:7781061/92‑1 (MQ=255)
                                   ccAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGttcttc                                                      >  2:625814/1‑93 (MQ=255)
                                   ccAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGttcttc                                                      <  1:4298664/93‑1 (MQ=255)
                                   ccAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGttcttc                                                      <  2:5139289/93‑1 (MQ=255)
                                    cAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCg                                                     >  2:7978332/1‑93 (MQ=255)
                                    cAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCg                                                     >  1:1369284/1‑93 (MQ=255)
                                     aGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGt                                                    <  2:7442421/93‑1 (MQ=255)
                                     aGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGt                                                    >  2:6871751/1‑93 (MQ=255)
                                     aGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGt                                                    <  2:4638999/93‑1 (MQ=255)
                                     aGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGt                                                    >  1:7645530/1‑93 (MQ=255)
                                      gCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTc                                                   >  1:5317951/1‑93 (MQ=255)
                                       caacaaAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCa                                                  <  2:4111660/93‑1 (MQ=255)
                                           aaaaCGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGtt                                              >  1:3875144/1‑93 (MQ=255)
                                             aaCGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTg                                            <  2:812741/93‑1 (MQ=255)
                                             aaCGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTg                                            <  1:8115146/93‑1 (MQ=255)
                                             aaCGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTg                                            >  2:3895788/1‑93 (MQ=255)
                                                 ccGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCTTGTTTGACCa                                        >  2:1550110/1‑93 (MQ=255)
                                                 ccGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCa                                        <  2:4800825/93‑1 (MQ=255)
                                                 ccGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCa                                        <  1:6557917/93‑1 (MQ=255)
                                                  cGTTAATGTCGCGGCCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCAc                                       <  1:5144207/93‑1 (MQ=255)
                                                  cGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCAc                                       >  2:6640404/1‑93 (MQ=255)
                                                  cGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCAc                                       <  2:1886611/93‑1 (MQ=255)
                                                     tAATGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCt                                    <  1:5650046/93‑1 (MQ=255)
                                                       aTGTCGCGACCGCGACGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTgg                                  <  2:8394509/93‑1 (MQ=255)
                                                            gcgaccgcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTGCGTCATGTTTGACCACGCTGGTCAcc                             >  2:7293232/1‑93 (MQ=255)
                                                             cgaccgcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCa                            <  1:7958704/93‑1 (MQ=255)
                                                              gaccgcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAg                           >  2:7070419/1‑93 (MQ=255)
                                                                ccgcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGt                         >  2:5798967/1‑93 (MQ=255)
                                                                 cgcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGtt                        >  2:2840868/1‑93 (MQ=255)
                                                                  gcgacGACGAGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGttt                       <  2:7176956/93‑1 (MQ=255)
                                                                    gacgacgaGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGa                     <  1:3051178/93‑1 (MQ=255)
                                                                      cgacgaGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACa                   >  1:6364033/1‑93 (MQ=255)
                                                                      cgacgaGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACa                   <  1:1761758/93‑1 (MQ=255)
                                                                       gacgaGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACAt                  >  2:2361609/1‑93 (MQ=255)
                                                                       gacgaGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACAt                  <  2:7728032/93‑1 (MQ=255)
                                                                        acgaGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATg                 <  1:1222454/93‑1 (MQ=255)
                                                                        acgaGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATg                 <  1:8246176/93‑1 (MQ=255)
                                                                        acgaGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATg                 <  1:7396527/93‑1 (MQ=255)
                                                                          gaGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATgcg               >  2:3078219/1‑93 (MQ=255)
                                                                           aGGACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATgcgc              >  2:6553775/1‑93 (MQ=255)
                                                                            ggACGACTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCa             >  1:1665152/1‑93 (MQ=255)
                                                                             gacgacTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCAt            >  2:3545933/1‑93 (MQ=255)
                                                                             gacgacTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTCCGCCTTGTTTGACCACCCGGGTAACAGGGTTTGACATGCCCac            >  2:4402677/1‑92 (MQ=255)
                                                                                gacTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATccc         <  2:4915200/93‑1 (MQ=255)
                                                                                gacTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATccc         >  2:4490969/1‑93 (MQ=255)
                                                                                gacTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATccc         <  2:3302473/93‑1 (MQ=255)
                                                                                cacTCATTAGCCCTCCTTGCTGTCGGCCGGGTTAACACGACGGTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATccc         <  2:2311688/92‑1 (MQ=255)
                                                                                 acTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCt        <  1:4113492/93‑1 (MQ=255)
                                                                                 acTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCt        >  1:3655949/1‑93 (MQ=255)
                                                                                  cTCATTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCtt       <  1:6570790/93‑1 (MQ=255)
                                                                                     aTTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCaa    <  2:2188259/93‑1 (MQ=255)
                                                                                     aTTAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCaa    >  2:6005236/1‑93 (MQ=255)
                                                                                      ttAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAAc   >  1:3995061/1‑93 (MQ=255)
                                                                                      ttAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAAc   <  1:7326078/93‑1 (MQ=255)
                                                                                      ttAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAAc   <  1:2601615/93‑1 (MQ=255)
                                                                                      ttAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAAc   >  1:6407677/1‑93 (MQ=255)
                                                                                       tAGCCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAAcc  <  1:6721489/93‑1 (MQ=255)
                                                                                          |                                                                                         
CAGCGCATCGTTGCTGGACATACCCTGAGGTTTATCCAGCAACAAAACGCCGTTAATGTCGCGACCGCGACGACGAGGACGACTCATTAGTCCTCCTTGCTGTCGTCCGGGTTAACACGACGTTCTTCGTCATGTTTGACCACGCTGGTCACCAGGTTTGACATGCGCATCCCTTCAACC  >  NC_000913/3310332‑3310511

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