Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,980,744 A→T 17.1% intergenic (‑109/+20) ygeA ← / ← araE Asp/Glu_racemase family protein/arabinose transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,980,7440AT17.1% 256.8 / 43.3 105intergenic (‑109/+20)ygeA/araEAsp/Glu_racemase family protein/arabinose transporter
Reads supporting (aligned to +/- strand):  ref base A (47/40);  new base T (10/8);  total (57/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.88e-01

Only 100 of 112 total aligned reads displayed.

CTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCG  >  NC_000913/2980645‑2980843
                                                                                                    |                                                                                                   
cTGATGTGTGTTACCGCATCAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACa                                                                                                     <  2:2131437/100‑1 (MQ=255)
cTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACa                                                                                                     <  1:2321052/100‑1 (MQ=255)
cTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACa                                                                                                     <  2:2508217/100‑1 (MQ=255)
cTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACa                                                                                                     >  2:1722253/1‑100 (MQ=255)
  gATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGc                                                                                                   <  2:573812/100‑1 (MQ=255)
  gATGTGTGTTACCCCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAg                                                                                                    <  1:1764162/99‑1 (MQ=255)
    tgtgtgTTACCGCAACAAGGCTGAATAAGGAACGAGACAAGCGGCTCAACGGCCAAGTGCTCAATCTCTATTAAAGAAAAAAGGGCCGGATGTACAGcac                                                                                                 <  1:180338/100‑1 (MQ=255)
    tgtgtgTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGcac                                                                                                 >  1:6970/1‑100 (MQ=255)
     gtgtgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGcaca                                                                                                >  2:509469/1‑100 (MQ=255)
       gtgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATc                                                                                              >  1:1206656/1‑100 (MQ=255)
       gtgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATc                                                                                              >  2:1886840/1‑100 (MQ=255)
           tACCGCAACAAGGCTGACTACGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGCAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                          <  2:1943816/100‑1 (MQ=255)
           tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                          <  1:259313/100‑1 (MQ=255)
           tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGc                                                                                          >  2:1377683/1‑100 (MQ=255)
            aaagaaaaaaGGCTGAATAAGGAACGAGACAAACACATCAACGGCCAAATGCCCAAAAAAAAAAAAAAAAAAAAAGGGCCGGATGTACAGCACATCCGGc                                                                                          <  1:511583/92‑1 (MQ=255)
            aCCGCAACAAGGCTGATTAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                          <  1:839874/99‑1 (MQ=38)
            aCCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGcc                                                                                         <  1:2421197/100‑1 (MQ=39)
            aCCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGcc                                                                                         <  2:936469/100‑1 (MQ=39)
             ccGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGccc                                                                                        <  2:1730504/100‑1 (MQ=39)
             ccGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCAGGATGTACAGCACATCCGGcc                                                                                         >  2:2352696/1‑99 (MQ=255)
              cGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCg                                                                                       >  1:1924779/1‑100 (MQ=255)
              cGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCg                                                                                       >  2:1343949/1‑100 (MQ=255)
               gCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGt                                                                                      >  2:1360813/1‑100 (MQ=255)
               gCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGt                                                                                      >  1:913802/1‑100 (MQ=255)
                 aacaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGCCCCGTGa                                                                                    <  2:1362703/100‑1 (MQ=255)
                   caaGGCTGAATAAGGAACGAGACAAACGCCGCAACGGCCAGGTGCCCAATCTCTATTAACGAATAAAGGGCCGGATGTGCTGTACATCCGGccctttttt                                                                                  <  1:454209/100‑7 (MQ=25)
                     aGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATc                                                                                >  2:2511345/1‑100 (MQ=255)
                     aGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATc                                                                                >  1:1310510/1‑100 (MQ=255)
                     aGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGAGCCGGATGTACAGCACATCCCGcccttttttcg                                                                                <  2:1707342/100‑9 (MQ=255)
                       gCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAg                                                                              <  2:192646/100‑1 (MQ=255)
                         tGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGAc                                                                            >  1:1640187/1‑100 (MQ=255)
                          gAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACg                                                                           >  2:216652/1‑100 (MQ=255)
                          gAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGACCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAg                                                                              >  1:1740822/1‑97 (MQ=255)
                            aTAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCTAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGcc                                                                         <  2:1866122/100‑1 (MQ=255)
                            aTAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGcc                                                                         >  2:1855680/1‑100 (MQ=255)
                                ggAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGata                                                                     >  1:1267340/1‑100 (MQ=255)
                                    cGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTtc                                                                 >  1:2433725/1‑100 (MQ=255)
                                     gggaCAAACGCCTCAACGGCCAAGTGCCCAATCTCTATGACCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTtc                                                                 <  1:1371631/97‑1 (MQ=255)
                                       gaCAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCa                                                              <  2:2037655/100‑1 (MQ=255)
                                       gaCAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCa                                                              <  2:1015993/100‑1 (MQ=255)
                                         cAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAAc                                                            <  1:2284964/100‑1 (MQ=255)
                                         cAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAAc                                                            >  1:2489687/1‑100 (MQ=255)
                                            aCGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTtc                                                         >  1:2261821/1‑100 (MQ=255)
                                             cGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTtct                                                        <  2:2131479/100‑1 (MQ=255)
                                              gCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTtctc                                                       >  1:1788605/1‑100 (MQ=255)
                                              gCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTtctc                                                       >  2:765993/1‑100 (MQ=255)
                                               ccTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCg                                                      >  2:2559517/1‑100 (MQ=255)
                                                cTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGc                                                     >  1:258200/1‑100 (MQ=255)
                                                  cAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCt                                                   >  2:704061/1‑100 (MQ=255)
                                                  cAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCt                                                   >  1:2074801/1‑100 (MQ=255)
                                                     cGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGcc                                                <  2:1026651/100‑1 (MQ=255)
                                                     cGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGcc                                                <  2:1380205/100‑1 (MQ=255)
                                                      ggCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCa                                               <  2:775982/100‑1 (MQ=255)
                                                      ggCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCa                                               >  2:1427970/1‑100 (MQ=255)
                                                       gCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCAt                                              <  1:1377675/100‑1 (MQ=255)
                                                       gCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCAt                                              <  1:1886830/100‑1 (MQ=255)
                                                              gCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTg                                       >  2:320434/1‑100 (MQ=255)
                                                              gCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTg                                       <  1:697334/100‑1 (MQ=255)
                                                               cccAATCTCTATTAACGAAAAAATGGCCGTATGTACAGCACATCCGGCCCGTGAAATCAGACTCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGc                                      <  2:482146/100‑1 (MQ=255)
                                                               cccAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGc                                      <  2:1299526/100‑1 (MQ=255)
                                                                ccAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCg                                     <  1:1738329/100‑1 (MQ=255)
                                                                  aaTCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGtt                                   >  2:332626/1‑100 (MQ=255)
                                                                  aaTCTCTATTAACGAAAAAAGAGCCGGATGTACAGCACATACGGGCCGTGAGACCAGCCGGAGATCTTCCTCAACCGCTAGCAAGCCATCAGt                                          >  1:2024679/1‑93 (MQ=255)
                                                                    tctcTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCg                                 <  2:311923/100‑1 (MQ=255)
                                                                     ctctATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGa                                <  2:614505/100‑1 (MQ=255)
                                                                      tctATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGat                               >  1:994977/1‑100 (MQ=255)
                                                                            aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTc                         <  1:842007/100‑1 (MQ=255)
                                                                            aaCGAAAAAAGGGCCGGATGTACAGCACATCCGACCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGt                           >  1:304795/1‑98 (MQ=255)
                                                                             aCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTcc                        <  2:265213/100‑1 (MQ=255)
                                                                             aCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTcc                        >  2:2238412/1‑100 (MQ=255)
                                                                              cGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCa                       >  2:1076705/1‑100 (MQ=255)
                                                                               gAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCat                      >  1:516617/1‑99 (MQ=255)
                                                                               gAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAg                      >  1:197542/1‑100 (MQ=255)
                                                                               gAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCCCGCCTGCCATCAGTTTGCGGTCGATATGTTCCa                       >  1:1157874/1‑99 (MQ=255)
                                                                                 aaaaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGc                     <  1:353502/100‑1 (MQ=255)
                                                                                  aaaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCg                    <  1:1268462/100‑1 (MQ=255)
                                                                                   aaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGt                   >  2:743452/1‑100 (MQ=255)
                                                                                    aaaGGGCCGGATGTGCAGAACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTg                  <  2:411372/100‑1 (MQ=255)
                                                                                    aaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTg                  <  2:50415/100‑1 (MQ=255)
                                                                                    aaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTg                  <  2:398628/100‑1 (MQ=255)
                                                                                      cgggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGAc                >  2:399421/2‑100 (MQ=39)
                                                                                      atggCCGGATGTACATCACATCCGGCACGTGATATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTg                  <  1:2094406/96‑1 (MQ=255)
                                                                                       gggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACa               >  1:1211727/1‑100 (MQ=39)
                                                                                        ggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCTGATATTTCTCCACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAt              >  1:1705235/1‑100 (MQ=38)
                                                                                        ggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAt              >  2:1870669/1‑100 (MQ=39)
                                                                                        ggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAt              >  2:85025/1‑100 (MQ=39)
                                                                                        ggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAt              >  1:2084170/1‑100 (MQ=39)
                                                                                        ggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAt              >  1:744264/1‑100 (MQ=255)
                                                                                        ggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAt              >  2:2463355/1‑100 (MQ=255)
                                                                                         gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtt             >  2:1654182/1‑100 (MQ=255)
                                                                                         gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtt             >  2:2246383/1‑100 (MQ=255)
                                                                                         gCCGGATGTACAGCACATACCGCGCGGGAAATCTATAGCCGCTACATCTCAACTTAACGTCCGCCATCAGGTTACGTTCACAATGCTCCAGCGTcaca               >  1:1904512/1‑98 (MQ=255)
                                                                                           cGGATGTACAGCACATCCGGCCCGTGGAATCAGACGCCGATATTTCGCAACTTCTCGCCTGCCATCAGTTTGCGTTCTATATGTTCCAGCGTGACAtttt           <  1:150883/100‑1 (MQ=255)
                                                                                           cGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtttt           <  1:16215/100‑1 (MQ=255)
                                                                                                gTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGtt      >  1:2376291/1‑100 (MQ=255)
                                                                                                 tACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCCGCGTGACATTTTTGGttt     >  1:501345/1‑100 (MQ=255)
                                                                                                  aCAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTc    >  2:1467460/1‑100 (MQ=255)
                                                                                                    aGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCg  >  2:1426285/1‑100 (MQ=255)
                                                                                                    aGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCg  >  2:74424/1‑100 (MQ=255)
                                                                                                    aGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCg  >  2:84814/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
CTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCG  >  NC_000913/2980645‑2980843

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