Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,501,888 +C coding (618/2496 nt) b2383 ← putative PTS system enzyme IIA component, enzyme I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,501,8841.C99.4% 553.8 / ‑3.9 174coding (622/2496 nt)b2383putative PTS system enzyme IIA component, enzyme I
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (70/103);  total (71/103)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 196 total aligned reads displayed.

AATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGT  >  NC_000913/2501792‑2501971
                                                                                             |                                                                                        
aaTGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAg                                                                                           <  1:3220178/93‑1 (MQ=255)
 aTGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGc                                                                                         >  1:3872354/1‑93 (MQ=255)
  tGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGcc                                                                                        >  1:1394153/1‑93 (MQ=255)
  tGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGcc                                                                                        >  1:447774/1‑93 (MQ=255)
     cGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGccccc                                                                                    <  1:897767/93‑2 (MQ=255)
       cTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGccccctg                                                                                  >  2:2211307/1‑90 (MQ=255)
       cTAACACGTTCGCGCAGATAATCGCTGACAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGccccctg                                                                                  >  1:1945222/1‑90 (MQ=255)
        tAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGccccctga                                                                                 <  2:2273883/93‑5 (MQ=255)
          acacGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATg                                                                                <  1:3733250/93‑1 (MQ=255)
          acacGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATg                                                                                >  2:4456427/1‑93 (MQ=255)
          acacGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATg                                                                                >  2:2055800/1‑93 (MQ=255)
           cacGTCCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGc                                                                               <  2:462727/93‑1 (MQ=255)
            acGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCt                                                                              >  2:552246/1‑93 (MQ=255)
            acGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCt                                                                              >  1:2289756/1‑93 (MQ=255)
                tCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTc                                                                          <  1:5834491/93‑1 (MQ=255)
                   cgcAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGt                                                                       <  1:4222224/93‑1 (MQ=255)
                   cgcAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGt                                                                       <  1:1744716/93‑1 (MQ=255)
                     cAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGtca                                                                     <  2:2405586/93‑1 (MQ=255)
                     cAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGtca                                                                     <  2:1643824/93‑1 (MQ=255)
                     cAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGtca                                                                     >  2:5343368/1‑93 (MQ=255)
                      aGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGtcat                                                                    >  2:1872819/1‑93 (MQ=255)
                       gATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGtcatc                                                                   >  2:2622902/1‑93 (MQ=255)
                       gATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGtcatc                                                                   <  2:1959710/93‑1 (MQ=255)
                        aTAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGtcatca                                                                  >  1:4837389/1‑93 (MQ=255)
                         tAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAg                                                                 <  2:3758683/93‑1 (MQ=255)
                         tAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAg                                                                 >  2:824283/1‑93 (MQ=255)
                          aaTCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAgg                                                                <  1:4973428/93‑1 (MQ=255)
                            tCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCg                                                              >  2:4402326/1‑93 (MQ=255)
                            tCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCg                                                              >  1:2975219/1‑93 (MQ=255)
                              gCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGAc                                                            >  1:3534406/1‑93 (MQ=255)
                               cTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACg                                                           >  1:3662269/1‑93 (MQ=255)
                               cTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACg                                                           <  1:6136229/93‑1 (MQ=255)
                               cTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACg                                                           <  2:3967383/93‑1 (MQ=255)
                                tGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACgg                                                          <  1:3416202/93‑1 (MQ=255)
                                tGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACgg                                                          >  1:1404000/1‑93 (MQ=255)
                                  gcagaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGat                                                        <  2:4705507/93‑1 (MQ=255)
                                  gcagaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGat                                                        <  1:1472458/93‑1 (MQ=255)
                                   cagaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGata                                                       >  1:972660/1‑93 (MQ=255)
                                     gaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATAtt                                                     >  2:1143235/1‑93 (MQ=255)
                                     gaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATAtt                                                     <  1:6081151/93‑1 (MQ=255)
                                       ggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGc                                                   <  2:1815225/93‑1 (MQ=255)
                                       ggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGc                                                   <  1:5523647/93‑1 (MQ=255)
                                        gcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGcc                                                  <  2:831672/93‑1 (MQ=255)
                                         cagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCt                                                 >  2:3451661/1‑93 (MQ=255)
                                          agaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTg                                                >  2:2836130/1‑93 (MQ=255)
                                           gaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGc                                               <  2:3562673/93‑1 (MQ=255)
                                           gaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGc                                               <  1:2888544/93‑1 (MQ=255)
                                           gaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGc                                               <  1:2941386/93‑1 (MQ=255)
                                            aaaGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCa                                              <  1:5208931/93‑1 (MQ=255)
                                                   ggCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATtcat                                       >  1:2857735/1‑93 (MQ=255)
                                                    gcgcAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATtcatc                                      >  1:5697526/1‑93 (MQ=255)
                                                       cAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCCGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCatcatc                                   <  1:3613670/93‑1 (MQ=255)
                                                        aaaCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATcc                                  >  2:5011162/1‑93 (MQ=255)
                                                        aaaCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATcc                                  <  2:4508783/93‑1 (MQ=255)
                                                        aaaCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATcc                                  >  1:4720073/1‑93 (MQ=255)
                                                         aaCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCt                                 <  1:3077064/93‑1 (MQ=255)
                                                         aaCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCt                                 >  2:4246262/1‑93 (MQ=255)
                                                         aaCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCt                                 <  1:4012791/93‑1 (MQ=255)
                                                          aCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGACTGCAAATTCATCATCCTg                                >  2:5753280/1‑93 (MQ=255)
                                                            cTGCTCCATATTGCTGATGATGGCCGCCCCCAGCCCACTGATGCTGTTCTGTCATCAGGCGAAGGATATTGCCTGCACATTCATCATGCTGaa                              >  2:2432216/1‑93 (MQ=255)
                                                            cTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGaa                              <  1:3174382/93‑1 (MQ=255)
                                                            cTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGaa                              <  1:2415636/93‑1 (MQ=255)
                                                                tCCATATTGCTGATGATCGCCGCCCCCCGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATGTTGCCTGCAAATTCATCATCCTGAATCAg                          >  2:4186509/1‑93 (MQ=255)
                                                                   atatTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGa                       >  1:3314442/1‑93 (MQ=255)
                                                                   atatTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGa                       >  2:3847163/1‑93 (MQ=255)
                                                                   atatTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGa                       <  2:605668/93‑1 (MQ=255)
                                                                   atatTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGa                       >  1:5313739/1‑93 (MQ=255)
                                                                   atacTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGa                       <  2:2991477/93‑1 (MQ=255)
                                                                    tatTGCTGATGTTCGCCTCCTCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCGGCAAATTCATCATCCTGAATCAGCGAc                      <  2:3610301/93‑1 (MQ=255)
                                                                    tatTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGAc                      >  2:5987423/1‑93 (MQ=255)
                                                                    tatTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGAc                      >  2:1518346/1‑93 (MQ=255)
                                                                       tGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACaaa                   <  1:940817/93‑1 (MQ=255)
                                                                        gCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAAt                  <  1:3624869/93‑1 (MQ=255)
                                                                        gCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAAt                  <  2:5144893/93‑1 (MQ=255)
                                                                         cTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATg                 <  2:5652942/93‑1 (MQ=255)
                                                                         cTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATg                 <  2:2184483/93‑1 (MQ=255)
                                                                          tgatgaTCGCCTCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATgg                <  1:5177871/93‑1 (MQ=255)
                                                                          tgatgaTCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATgg                <  1:4849660/93‑1 (MQ=255)
                                                                          tgatgaTCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATgg                <  2:4788772/93‑1 (MQ=255)
                                                                           gatgatCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATtgg               >  1:4608810/1‑93 (MQ=255)
                                                                           gatgatCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATggg               >  2:1839315/1‑93 (MQ=255)
                                                                           gatgatCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATggg               <  1:4160055/93‑1 (MQ=255)
                                                                           gatgatCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATggg               <  2:3261185/93‑1 (MQ=255)
                                                                           gatgatCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATggg               <  1:469043/93‑1 (MQ=255)
                                                                           gatgatCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATggg               >  1:755776/1‑93 (MQ=255)
                                                                           gatgatCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATggg               >  1:5975839/1‑93 (MQ=255)
                                                                            atgatCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGgc              <  1:655396/93‑1 (MQ=255)
                                                                              gatCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGgcgc            <  2:2918543/93‑1 (MQ=255)
                                                                               atCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCt           <  2:5214129/93‑1 (MQ=255)
                                                                                 cgccgcCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGa         <  2:5933003/93‑1 (MQ=255)
                                                                                  gccgccCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAg        <  2:2727094/93‑1 (MQ=255)
                                                                                  gccgccCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAg        <  1:1641464/93‑1 (MQ=255)
                                                                                    cgccCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGa      <  2:4929876/93‑1 (MQ=255)
                                                                                     gccCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGAt     <  2:4894872/93‑1 (MQ=255)
                                                                                       ccccAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAg   >  1:5422696/1‑93 (MQ=255)
                                                                                       ccccAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAg   >  1:4542290/1‑93 (MQ=255)
                                                                                       ccccAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAg   >  1:5098733/1‑93 (MQ=255)
                                                                                        cccAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGt  >  1:151901/1‑93 (MQ=255)
                                                                                        cccAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGt  >  1:6259278/1‑93 (MQ=255)
                                                                                        cccAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCACCCGGAATCAGCGACAAATGGGCGCTGAGGATAGt  <  2:3321035/93‑1 (MQ=255)
                                                                                             |                                                                                        
AATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGT  >  NC_000913/2501792‑2501971

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