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.C98.6% 690.9 / 1.2 216coding (622/2496 nt)b2383putative PTS system enzyme IIA component, enzyme I
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (112/101);  total (114/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 233 total aligned reads displayed.

AATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGT  >  NC_000913/2501792‑2501971
                                                                                             |                                                                                        
aaTGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGGGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAg                                                                                           >  1:583336/1‑93 (MQ=255)
 aTGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGc                                                                                         >  1:1232220/1‑93 (MQ=255)
 aTGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGc                                                                                         >  2:3583411/1‑93 (MQ=255)
   gTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGccc                                                                                      >  2:6267662/1‑93 (MQ=255)
     cGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGccccc                                                                                    <  1:201130/93‑2 (MQ=255)
       cTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGccccctg                                                                                  >  1:2993591/1‑90 (MQ=255)
       cTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGccccctg                                                                                  <  1:4706617/93‑4 (MQ=255)
           cacGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGc                                                                               <  2:2090445/93‑1 (MQ=255)
            acGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCt                                                                              >  1:171185/1‑93 (MQ=255)
             cGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTg                                                                             >  2:1178313/1‑93 (MQ=255)
              gTTCGCGCAGAGAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCCCCCCCAGCCCCCTGATGCTGt                                                                            <  2:3059040/93‑1 (MQ=255)
               ttCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGtt                                                                           >  1:2468164/1‑93 (MQ=255)
                tCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTc                                                                          <  2:4427810/93‑1 (MQ=255)
                   cgcAAATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCTCCCCCAGCCCCCTGATGCTGTTCTGt                                                                       <  1:603450/93‑1 (MQ=255)
                    gcAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGtc                                                                      <  2:3641011/93‑1 (MQ=255)
                       gATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGtcatc                                                                   <  1:4929847/93‑1 (MQ=255)
                       gATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGtcatc                                                                   >  1:5245073/1‑93 (MQ=255)
                       gATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGtcatc                                                                   >  2:1903372/1‑93 (MQ=255)
                       gATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGtcatc                                                                   >  2:6065012/1‑93 (MQ=255)
                         tAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAg                                                                 <  1:5033787/93‑1 (MQ=255)
                         tAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAg                                                                 >  2:738255/1‑93 (MQ=255)
                           aTCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCTCCAGCCCCCTGATGCTGTTCTGTCATCAGGc                                                               >  2:3969549/1‑93 (MQ=255)
                           aTCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGc                                                               >  1:2348023/1‑93 (MQ=255)
                            tCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCg                                                              >  1:2909661/1‑93 (MQ=255)
                            tCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCg                                                              >  2:1515947/1‑93 (MQ=255)
                             cGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGa                                                             <  2:1395145/93‑1 (MQ=255)
                              gCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGAc                                                            >  1:2611884/1‑93 (MQ=255)
                               cTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACg                                                           >  1:378134/1‑93 (MQ=255)
                               cTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACg                                                           >  2:3841620/1‑93 (MQ=255)
                                 ggcagaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGa                                                         <  1:6647423/93‑1 (MQ=255)
                                  gcagaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGat                                                        <  2:2495660/93‑1 (MQ=255)
                                  gcagaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCCGTTCTGTCATCAGGCGACGGat                                                        <  2:4646556/93‑1 (MQ=255)
                                  gcagaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCATCCTGATGCTGTTCTGTCATCAGGCGACGGat                                                        >  1:4897378/1‑93 (MQ=255)
                                    agaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGatat                                                      >  2:6178111/1‑93 (MQ=255)
                                     gaggcagaGAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATAtt                                                     >  2:2534635/1‑93 (MQ=255)
                                     gaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATAtt                                                     <  2:6388955/93‑1 (MQ=255)
                                     gaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATAtt                                                     <  1:1108314/93‑1 (MQ=255)
                                     gaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATAtt                                                     >  1:5168640/1‑93 (MQ=255)
                                      aggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTg                                                    >  2:5067389/1‑93 (MQ=255)
                                       ggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGc                                                   >  2:6198148/1‑93 (MQ=255)
                                          agaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTg                                                >  1:2759571/1‑93 (MQ=255)
                                           gaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGc                                               >  1:534155/1‑93 (MQ=255)
                                           gaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGc                                               >  1:4181161/1‑93 (MQ=255)
                                           gaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGc                                               >  1:722029/1‑93 (MQ=255)
                                           gaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGc                                               >  1:6689756/1‑93 (MQ=255)
                                           aaaaTTTTTGCGCAAACCTGCTCCATATGGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGc                                               <  1:3590472/92‑1 (MQ=255)
                                            aaaGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCa                                              >  2:1631359/1‑93 (MQ=255)
                                               gTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAAt                                           >  2:4567576/1‑93 (MQ=255)
                                                tttGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAAtt                                          >  1:3610493/1‑93 (MQ=255)
                                                tttGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAAtt                                          >  1:4360779/1‑93 (MQ=255)
                                                  tGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATtca                                        <  2:5013182/93‑1 (MQ=255)
                                                  tGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATtca                                        <  2:6620578/93‑1 (MQ=255)
                                                    gcgcAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATtcatc                                      >  1:392831/1‑93 (MQ=255)
                                                    gcgcAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATtcatc                                      >  2:1416253/1‑93 (MQ=255)
                                                    gcgcAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATtcatc                                      >  1:4058188/1‑93 (MQ=255)
                                                     cgcAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATtcatca                                     <  2:3448976/93‑1 (MQ=255)
                                                      gcAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTcatcat                                    >  2:2051620/1‑93 (MQ=255)
                                                      gcAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTcatcat                                    >  1:4815472/1‑93 (MQ=255)
                                                        aaaCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATcc                                  >  1:3852029/1‑93 (MQ=255)
                                                        aaaCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATcc                                  <  2:302076/93‑1 (MQ=255)
                                                         aaCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCt                                 >  2:95117/1‑93 (MQ=255)
                                                         aaCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCt                                 <  1:5124508/93‑1 (MQ=255)
                                                         aaCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCt                                 <  1:3206216/93‑1 (MQ=255)
                                                         aaCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCt                                 >  1:1794610/1‑93 (MQ=255)
                                                         aaCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCt                                 >  1:5941790/1‑93 (MQ=255)
                                                          aCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTg                                >  2:3340789/1‑93 (MQ=255)
                                                           ccTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGa                               >  2:1330940/1‑93 (MQ=255)
                                                            cTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGaa                              <  2:5530433/93‑1 (MQ=255)
                                                            cTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGaa                              <  2:4099245/93‑1 (MQ=255)
                                                             tGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAAt                             >  1:5411776/1‑93 (MQ=255)
                                                              gCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATc                            <  1:1343593/93‑1 (MQ=255)
                                                              gCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATc                            >  2:6409405/1‑93 (MQ=255)
                                                              gCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATc                            <  2:5321470/93‑1 (MQ=255)
                                                               cTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCa                           <  2:5810363/93‑1 (MQ=255)
                                                                  cATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCg                        >  2:4118708/1‑93 (MQ=255)
                                                                  cATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCg                        >  2:5636975/1‑93 (MQ=255)
                                                                    tatTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGAc                      >  1:2233731/1‑93 (MQ=255)
                                                                     atTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACa                     <  1:1563779/93‑1 (MQ=255)
                                                                         cTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATg                 <  2:5783696/93‑1 (MQ=255)
                                                                         cTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCACCATCCTGAATCAGCGACAAATg                 >  1:302936/1‑93 (MQ=255)
                                                                          tgatgaTCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATgg                <  2:2899901/93‑1 (MQ=255)
                                                                          tgatgaTCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATgg                <  1:6336373/93‑1 (MQ=255)
                                                                             tgatCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGgcg             >  2:6580846/1‑93 (MQ=255)
                                                                              gatCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGgcgc            >  2:6453138/1‑93 (MQ=255)
                                                                               atCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCt           <  1:1147681/93‑1 (MQ=255)
                                                                               atCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCt           <  2:2690501/93‑1 (MQ=255)
                                                                                 cgccgcCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGa         <  2:2233735/93‑1 (MQ=255)
                                                                                 cgccgcCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGa         >  2:1291879/1‑93 (MQ=255)
                                                                                 cgccgcCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGa         <  2:670314/93‑1 (MQ=255)
                                                                                  gccgccCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAg        <  1:88252/93‑1 (MQ=255)
                                                                                   ccgccCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAgg       <  2:4784445/93‑1 (MQ=255)
                                                                                   ccgccCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAgg       <  1:4629209/93‑1 (MQ=255)
                                                                                     gccCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGAt     >  1:5904476/1‑93 (MQ=255)
                                                                                     gccCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGAt     >  1:3707096/1‑93 (MQ=255)
                                                                                     ccccccAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGAt     <  2:1152406/92‑1 (MQ=255)
                                                                                      cccccAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATa    >  2:4456017/1‑93 (MQ=255)
                                                                                      cccccAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATa    <  2:3693830/93‑1 (MQ=255)
                                                                                      cccccAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATa    <  1:3995059/93‑1 (MQ=255)
                                                                                       ccccAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAg   <  2:1753160/93‑1 (MQ=255)
                                                                                       ccccAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAg   <  1:3727866/93‑1 (MQ=255)
                                                                                             |                                                                                        
AATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGT  >  NC_000913/2501792‑2501971

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