Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,862,881 T→C L208P (CTA→CCA)  ygbM → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,862,8810TC99.4% 537.4 / NA 171L208P (CTA→CCA) ygbMhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (90/80);  minor base G (1/0);  total (91/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

ATACTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACTACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGATCG  >  NC_000913/2862791‑2862970
                                                                                          |                                                                                         
 tACTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGaccacc                                                                                        >  2:728856/1‑93 (MQ=255)
  aCTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCt                                                                                       <  2:924523/93‑1 (MQ=255)
   cTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTg                                                                                      >  1:1469868/1‑93 (MQ=255)
    ttttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGa                                                                                     >  1:4401890/1‑93 (MQ=255)
      ttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATc                                                                                   >  2:2993075/1‑93 (MQ=255)
        cATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGa                                                                                 <  1:1345211/93‑1 (MQ=255)
          tGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACa                                                                               <  1:2791959/93‑1 (MQ=255)
           gcgcAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACAt                                                                              <  2:2559407/93‑1 (MQ=255)
           gcgcAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACAt                                                                              <  2:1676640/93‑1 (MQ=255)
             gcAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGa                                                                            <  1:99224/93‑1 (MQ=255)
             gcAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGa                                                                            >  1:3022316/1‑93 (MQ=255)
               aaaaaGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAAc                                                                          >  1:1059926/1‑93 (MQ=255)
                aaaaGTAGATGGCACCCTGACGCATTTACTCAGCTCCTCCGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAAcc                                                                         <  1:1933813/93‑1 (MQ=255)
                aaaaGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAAcc                                                                         >  2:1765310/1‑93 (MQ=255)
                  aaGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCgg                                                                       >  1:2432383/1‑93 (MQ=255)
                   aGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGa                                                                      >  1:511966/1‑93 (MQ=255)
                   aGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGa                                                                      <  2:4246942/93‑1 (MQ=255)
                   aGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGa                                                                      >  2:2291471/1‑93 (MQ=255)
                    gTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGAc                                                                     >  1:2097716/1‑93 (MQ=255)
                    gTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGAc                                                                     >  1:3365817/1‑93 (MQ=255)
                      aGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGa                                                                   <  2:269404/93‑1 (MQ=255)
                       gATTGCAACCTGACTCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGAt                                                                  <  2:1586683/93‑1 (MQ=255)
                       gATTGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGAt                                                                  <  1:1025286/93‑1 (MQ=255)
                       gATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGAt                                                                  <  2:2108659/93‑1 (MQ=255)
                        aTGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATg                                                                 <  2:804976/93‑1 (MQ=255)
                        aTGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATg                                                                 >  2:4014832/1‑93 (MQ=255)
                        aTGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATg                                                                 >  1:1674685/1‑93 (MQ=255)
                         tGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATgg                                                                >  2:3246293/1‑93 (MQ=255)
                         tGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATgg                                                                <  1:1263008/93‑1 (MQ=255)
                         tGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATgg                                                                >  1:3027215/1‑93 (MQ=255)
                         gggCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGCACCAATTGCCGGACCACCTGATCGACATGAACCGGACGATgg                                                                >  2:1634299/2‑93 (MQ=255)
                          ggCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGga                                                               <  2:3793322/93‑1 (MQ=255)
                          ggCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGga                                                               >  1:2131053/1‑93 (MQ=255)
                           gCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGgag                                                              >  2:1979902/1‑93 (MQ=255)
                            cAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACGACCTGATCGACATGAACCGGACGATGgaga                                                             >  1:406921/1‑93 (MQ=255)
                             aaCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCTACATTAACCGGACGATGGAGaa                                                            >  1:1237867/1‑93 (MQ=255)
                              aCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGaaa                                                           <  1:749601/93‑1 (MQ=255)
                              aCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGaaa                                                           >  2:4469047/1‑93 (MQ=255)
                              aCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGaaa                                                           >  1:3811940/1‑93 (MQ=255)
                               ccTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAAt                                                          >  1:1904518/1‑93 (MQ=255)
                                  gACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCaa                                                       >  1:1594124/1‑93 (MQ=255)
                                  gACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCaa                                                       <  1:1132190/93‑1 (MQ=255)
                                  gACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCaa                                                       >  1:492688/1‑93 (MQ=255)
                                   aCGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAAc                                                      <  2:3870892/93‑1 (MQ=255)
                                    cGCATTTTATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACt                                                     >  2:2626884/1‑93 (MQ=255)
                                     gCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTa                                                    <  1:1358579/93‑1 (MQ=255)
                                       aTTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATc                                                  <  1:1807447/93‑1 (MQ=255)
                                        tttAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATcc                                                 <  2:2082729/93‑1 (MQ=255)
                                         ttAATCTGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCg                                                >  2:3899520/1‑93 (MQ=255)
                                          tAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGt                                               <  2:2739170/93‑1 (MQ=255)
                                          tAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGt                                               <  1:3412980/93‑1 (MQ=255)
                                          tAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGt                                               <  1:3661766/93‑1 (MQ=255)
                                           aaTCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTg                                              <  2:1213379/93‑1 (MQ=255)
                                             tCCGCGGCTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGc                                            <  2:4341096/93‑1 (MQ=255)
                                             tCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGc                                            <  2:1835664/93‑1 (MQ=255)
                                             tCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATACACTATCCGTGGc                                            >  1:2668761/1‑93 (MQ=255)
                                                gcgACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGt                                         >  2:1324166/1‑93 (MQ=255)
                                                  gACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTc                                       >  1:3735548/1‑93 (MQ=255)
                                                  gACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTc                                       >  2:2215044/1‑93 (MQ=255)
                                                  gACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAAATATCCGTGGCTGTTc                                       >  1:267113/1‑93 (MQ=255)
                                                   aCTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTcc                                      <  2:4396623/93‑1 (MQ=255)
                                                     tACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGc                                    >  2:2519193/1‑93 (MQ=255)
                                                      aCGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGcc                                   <  2:4295512/93‑1 (MQ=255)
                                                        gCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTg                                 <  1:2159947/93‑1 (MQ=255)
                                                         ccGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGt                                >  1:1556069/1‑93 (MQ=255)
                                                          cGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGtt                               <  2:4525080/93‑1 (MQ=255)
                                                          cGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGtt                               >  1:1877104/1‑93 (MQ=255)
                                                           gggAAATATGCGCACGTACAAATTGCCGGACCGCCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTc                              <  1:3250104/93‑1 (MQ=255)
                                                           gggAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTc                              <  1:4208027/93‑1 (MQ=255)
                                                           gggAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTc                              >  2:499667/1‑93 (MQ=255)
                                                           gggAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTc                              >  1:2399010/1‑93 (MQ=255)
                                                            ggAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCg                             >  2:4640419/1‑93 (MQ=255)
                                                             gAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGa                            >  1:3878774/1‑93 (MQ=255)
                                                             gAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGa                            >  1:3725257/1‑93 (MQ=255)
                                                                 tatGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAg                        >  1:747564/1‑93 (MQ=255)
                                                                 tatGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAg                        <  2:3361460/93‑1 (MQ=255)
                                                                  atGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAgg                       >  2:2671032/1‑93 (MQ=255)
                                                                   tGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAggg                      >  1:7985/1‑92 (MQ=255)
                                                                   tGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGt                      <  2:143776/93‑1 (MQ=255)
                                                                   tGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGt                      >  2:3444732/1‑93 (MQ=255)
                                                                   tGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGt                      <  1:4042968/93‑1 (MQ=255)
                                                                      gcACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTggg                   >  1:3185004/1‑93 (MQ=255)
                                                                       cACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGa                  <  1:3289715/93‑1 (MQ=255)
                                                                         cGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGata                >  1:3317039/1‑93 (MQ=255)
                                                                          gTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGatat               >  2:3337946/1‑93 (MQ=255)
                                                                           tACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATc              <  1:849591/93‑1 (MQ=255)
                                                                           tACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATc              <  1:3135146/93‑1 (MQ=255)
                                                                             cAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAg            <  1:3639060/93‑1 (MQ=255)
                                                                             cAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAg            <  1:4454203/93‑1 (MQ=255)
                                                                             cAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAg            >  2:3371008/1‑93 (MQ=255)
                                                                             cAAATAGCCGGACCACCTGATCTACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAg            <  2:4065355/93‑1 (MQ=255)
                                                                              aaaTTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGCAATCAACTATCCGTGGCTGGTCCGCCTGTTCGCTGAGGCGGGatatgaag           >  1:2928877/1‑89 (MQ=255)
                                                                               aaTTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAggg          >  2:3205917/1‑93 (MQ=255)
                                                                                aTTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGc         <  2:1017706/93‑1 (MQ=255)
                                                                                aTTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGc         <  2:2470573/93‑1 (MQ=255)
                                                                                aTTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGc         >  2:722694/1‑93 (MQ=255)
                                                                                 ttGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCt        >  1:3503155/1‑93 (MQ=255)
                                                                                 ttGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGCGGGATATCAGGGCt        >  2:3556690/1‑93 (MQ=255)
                                                                                    ccGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGa     <  1:4490488/93‑1 (MQ=255)
                                                                                     cGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGAt    <  2:1382657/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATACTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACTACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGATCG  >  NC_000913/2862791‑2862970

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