Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,372,457 A→G V260A (GTG→GCG)  yjeK ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,372,4570AG99.0% 642.1 / ‑4.5 200V260A (GTG→GCG) yjeKhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (97/101);  total (98/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 = 1.00e+00

Only 100 of 203 total aligned reads displayed.

ACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCA  >  NC_000913/4372367‑4372546
                                                                                          |                                                                                         
 cTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCAt                                                                                        >  2:597682/1‑93 (MQ=255)
  tttATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATc                                                                                       >  2:3221075/1‑93 (MQ=255)
  tttATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATc                                                                                       >  2:1930646/1‑93 (MQ=255)
   ttATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCAtaa                                                                                      >  1:2067054/1‑91 (MQ=255)
   ttATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCa                                                                                      >  1:2829462/1‑93 (MQ=255)
   ttATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCa                                                                                      <  1:2860262/93‑1 (MQ=255)
    tATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCAc                                                                                     <  1:3994551/93‑1 (MQ=255)
     aTCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACg                                                                                    >  2:3456769/1‑93 (MQ=255)
     aTCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACg                                                                                    <  2:1987164/93‑1 (MQ=255)
       cGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCTTTTGTGCGTTATCGTTCGCATCACGTa                                                                                  <  1:3525848/93‑1 (MQ=255)
         aGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAAc                                                                                <  2:2602123/93‑1 (MQ=255)
         aGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAAc                                                                                >  2:4255057/1‑93 (MQ=255)
         aGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAAc                                                                                >  2:1825215/1‑93 (MQ=255)
         aGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAAc                                                                                >  2:1327477/1‑93 (MQ=255)
          gCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACa                                                                               >  1:2164812/1‑93 (MQ=255)
          gCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACa                                                                               >  2:1661167/1‑93 (MQ=255)
           cacaTGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAg                                                                              >  1:1002939/1‑93 (MQ=255)
           cacaTGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAg                                                                              >  1:191830/1‑93 (MQ=255)
             caTGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGaa                                                                            <  1:743495/93‑1 (MQ=255)
               tGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACTTAACAGAACg                                                                          >  2:1834279/1‑93 (MQ=255)
               tGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACg                                                                          >  2:290275/1‑93 (MQ=255)
                gCAGGTAATACGGCATTACTCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGc                                                                         <  1:4617982/93‑1 (MQ=255)
                gCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGc                                                                         <  1:3941349/93‑1 (MQ=255)
                gCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGc                                                                         >  2:84645/1‑93 (MQ=255)
                  aGGTAAAACGGCATTCAGCCGGCATCGAACAACGCATTACTCAGGTTTGCCGGCGTTTGTGCGTTATCGTTCGCATCCCGTAACAGAACGctc                                                                       <  2:3322180/93‑1 (MQ=255)
                     tAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTgg                                                                    <  1:4121790/93‑1 (MQ=255)
                      aaTACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGt                                                                   <  2:211904/93‑1 (MQ=255)
                       aTACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGtt                                                                  >  1:480202/1‑93 (MQ=255)
                        tACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTc                                                                 <  2:4120005/93‑1 (MQ=255)
                         tctgCAATACGGCGGCATCGAACAACGCCTTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCa                                                                <  1:1680974/90‑1 (MQ=255)
                         aCGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCa                                                                >  1:4152549/1‑93 (MQ=255)
                           ggCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGc                                                              >  2:1064025/1‑93 (MQ=255)
                           ggCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGc                                                              >  2:4214135/1‑93 (MQ=255)
                              aTTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCaaa                                                           >  2:613007/1‑93 (MQ=255)
                              aTTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCaaa                                                           >  2:1402034/1‑93 (MQ=255)
                               ttACGCCGTCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAg                                                          <  1:789676/93‑1 (MQ=255)
                               ttACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAg                                                          >  1:1201812/1‑93 (MQ=255)
                               ttACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAg                                                          >  1:4290193/1‑93 (MQ=255)
                                tACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAga                                                         >  2:3489937/1‑92 (MQ=255)
                                tACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGt                                                         <  2:193485/93‑1 (MQ=255)
                                caCGCCGGCATCGAACAACGCATTACGCAGGGTTGCCAGCGTTTGTGCGTTATCTTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGt                                                         <  1:1398740/92‑1 (MQ=255)
                                 aCGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTa                                                        <  2:3934891/93‑1 (MQ=255)
                                 aCGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTa                                                        >  1:1076266/1‑93 (MQ=255)
                                 aCGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTa                                                        >  2:1283015/1‑93 (MQ=255)
                                  cGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTaa                                                       >  2:4165768/1‑93 (MQ=255)
                                  cGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTaa                                                       <  2:4522976/93‑1 (MQ=255)
                                   gCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAAc                                                      <  1:2122563/93‑1 (MQ=255)
                                       gCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCt                                                  <  1:3350440/93‑1 (MQ=255)
                                       gCATCGAACAACGCATTACTCAGGTGTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCt                                                  >  2:2487488/1‑93 (MQ=255)
                                        cATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTa                                                 <  1:17349/93‑1 (MQ=255)
                                          tCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTAcc                                               >  1:4238794/1‑93 (MQ=255)
                                          tCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTAcc                                               >  1:4328631/1‑93 (MQ=255)
                                          tCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTAcc                                               >  2:1474119/1‑93 (MQ=255)
                                          tCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTAcc                                               >  2:4449217/1‑93 (MQ=255)
                                          tCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTAcc                                               >  1:2022733/1‑93 (MQ=255)
                                           cGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTAccc                                              <  2:1456342/93‑1 (MQ=255)
                                            gAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCg                                             >  2:3296161/1‑93 (MQ=255)
                                            gAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCg                                             >  2:2305404/1‑93 (MQ=255)
                                                aacGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCTgcgc                                         >  2:1770903/1‑93 (MQ=255)
                                                 acGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCa                                        <  1:3276357/93‑1 (MQ=255)
                                                  cGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCaa                                       <  1:3732237/93‑1 (MQ=255)
                                                   gCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAAc                                      >  1:318433/1‑93 (MQ=255)
                                                      ttACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTa                                   <  2:521237/93‑1 (MQ=255)
                                                      ttACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTa                                   <  2:3615788/93‑1 (MQ=255)
                                                        aCTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATAACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGc                                 <  2:1539897/93‑1 (MQ=255)
                                                         cTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGcc                                >  1:4650192/1‑93 (MQ=255)
                                                          tCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCa                               >  2:2191505/1‑93 (MQ=255)
                                                            cggTTTGCCAGCGTTTGTGCGTTCTCGTTCGCCTCACGTCCCGGCGCGCTCTGGTTCGGCTAAGTACCGCCTCCCCGGCGCAGCTTAGCCATc                             <  1:4222330/92‑1 (MQ=255)
                                                            aGGTTTGCCTGCGTTTTTGCCTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATc                             <  2:2758341/93‑1 (MQ=255)
                                                            aGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATGACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATc                             >  2:2974579/1‑93 (MQ=255)
                                                            aGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATc                             <  1:3234874/93‑1 (MQ=255)
                                                            aGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATc                             <  1:229492/93‑1 (MQ=255)
                                                             ggTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCg                            >  2:1570008/1‑93 (MQ=255)
                                                              gTTTGCCAGCGTTTTTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGc                           <  1:3242267/93‑1 (MQ=255)
                                                               tttGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGcc                          <  2:281324/93‑1 (MQ=255)
                                                                 tGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTg                        <  2:4067445/93‑1 (MQ=255)
                                                                  gCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGa                       <  1:1032620/93‑1 (MQ=255)
                                                                   ccAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGAc                      >  2:3733085/1‑93 (MQ=255)
                                                                    cAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACg                     >  2:3058090/1‑93 (MQ=255)
                                                                     aGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACgg                    <  2:1889780/93‑1 (MQ=255)
                                                                     aGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACgg                    >  2:2921831/1‑93 (MQ=255)
                                                                      gCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGa                   >  2:3853575/1‑93 (MQ=255)
                                                                       cGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCTATCGCCTGACGGaa                  <  1:3083753/93‑1 (MQ=255)
                                                                        gTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAAt                 <  1:2734409/93‑1 (MQ=255)
                                                                        gTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAAt                 >  2:4339991/1‑93 (MQ=255)
                                                                         tttGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATg                >  1:1107675/1‑93 (MQ=255)
                                                                         tttGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATg                >  2:4377026/1‑93 (MQ=255)
                                                                         tttGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATg                >  1:2420399/1‑93 (MQ=255)
                                                                           tgtgCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGtt              >  2:339921/1‑93 (MQ=255)
                                                                           tgtgCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGtt              <  1:1237097/93‑1 (MQ=255)
                                                                           tgtgCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGtt              <  1:4255057/93‑1 (MQ=255)
                                                                             tgCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTc            <  2:2417788/93‑1 (MQ=255)
                                                                              gCGTTATCGTTAGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCa           >  1:3669443/1‑93 (MQ=255)
                                                                                gTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGTAATGTttcttc         >  1:4519999/1‑93 (MQ=255)
                                                                                gTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATc         >  1:3702038/1‑93 (MQ=255)
                                                                                 gtATCGTTCGCATCACGGAACAGACCGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCt        <  1:743573/92‑1 (MQ=255)
                                                                                    tCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTAcc     <  1:597682/93‑1 (MQ=255)
                                                                                    tCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTAcc     <  1:2544471/93‑1 (MQ=255)
                                                                                    tCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTAcc     <  2:4693364/93‑1 (MQ=255)
                                                                                      gTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTc   >  2:1165744/1‑93 (MQ=255)
                                                                                          |                                                                                         
ACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCA  >  NC_000913/4372367‑4372546

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