Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 993,254 A→G 6.1% intergenic (+20/+23) pepN → / ← ssuB aminopeptidase N/aliphatic sulfonate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,2540AG 6.1% 417.3 / 10.4 131intergenic (+20/+23)pepN/ssuBaminopeptidase N/aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (64/59);  new base G (6/2);  total (70/61)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 133 total aligned reads displayed.

CGCTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTC  >  NC_000913/993107‑993400
                                                                                                                                                   |                                                                                                                                                  
cGCTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCaa                                                                                                                                                   >  1:1035492/1‑149 (MQ=255)
cGCTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCaa                                                                                                                                                   <  4:1401401/149‑1 (MQ=255)
  cTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATa                                                                                                                                                 <  2:338120/149‑1 (MQ=255)
    gATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAgc                                                                                                                                               <  4:795504/149‑1 (MQ=255)
     aTTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAgcg                                                                                                                                              <  1:632678/149‑1 (MQ=255)
      ttCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAgcgc                                                                                                                                             >  3:1278290/1‑149 (MQ=255)
         gCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGcgccgc                                                                                                                                          >  3:1045750/1‑149 (MQ=255)
              aaaCGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTc                                                                                                                                     <  3:121490/149‑1 (MQ=255)
               aaCGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCg                                                                                                                                    <  1:309560/149‑1 (MQ=255)
                aCGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCgg                                                                                                                                   <  4:660903/149‑1 (MQ=255)
                 cGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                                                                  <  2:960688/149‑1 (MQ=255)
                 cGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                                                                  <  1:1329849/149‑1 (MQ=255)
                    tACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGaa                                                                                                                               <  3:766426/149‑1 (MQ=255)
                     aCGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAAt                                                                                                                              <  1:811535/149‑1 (MQ=255)
                     aCGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAAt                                                                                                                              <  1:264856/149‑1 (MQ=255)
                       gATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAAttt                                                                                                                            >  3:1169922/1‑149 (MQ=255)
                            cAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAcccc                                                                                                                       >  2:3552/1‑149 (MQ=255)
                            cAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAccc                                                                                                                        >  1:675886/1‑148 (MQ=255)
                             aaaCGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCt                                                                                                                      >  3:157081/1‑149 (MQ=255)
                              aaCGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTg                                                                                                                     >  2:1067452/1‑149 (MQ=255)
                                cGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGtt                                                                                                                   <  2:1368441/149‑1 (MQ=255)
                                    aGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCAGCAATAGCGCCACCATTCGGGGAATTTACCCCTGTTTtct                                                                                                               >  4:1164500/1‑149 (MQ=255)
                                     ggAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTtctc                                                                                                              <  3:684497/149‑1 (MQ=255)
                                     ggAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTtctc                                                                                                              >  3:1318032/1‑149 (MQ=255)
                                        gaAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAcc                                                                                                                         >  2:695525/1‑135 (MQ=255)
                                        gaAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAcc                                                                                                                         <  1:695525/135‑1 (MQ=255)
                                           aaTGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGa                                                                                                        <  4:541530/149‑1 (MQ=255)
                                                gcgcgGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTc                                                                                                   <  3:303274/149‑1 (MQ=255)
                                                 cgcgGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCa                                                                                                  <  1:9504/149‑1 (MQ=255)
                                                   cggcgCTGGAACAGTTGAAAGGGCTGGGAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGa                                                                                                <  2:1209508/149‑1 (MQ=255)
                                                   cggcgCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGa                                                                                                <  3:749047/149‑1 (MQ=255)
                                                   cggcgCAGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATATCTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGa                                                                                                <  1:528945/149‑1 (MQ=255)
                                                          ggAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAACGCCGCCATTCGGGGAATTTACCACTGTTTTCTCAGGCGAGTTTCAGATTcacca                                                                                         >  1:434867/1‑149 (MQ=255)
                                                             aCAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGt                                                                                      >  4:1122260/1‑149 (MQ=255)
                                                              cAGTTGAAAGGGCTGGAATATCTCTCTGGCGATCTGTACTAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGtt                                                                                     <  1:337317/149‑1 (MQ=255)
                                                              cAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTTTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCTATTGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGtt                                                                                     <  4:108491/149‑1 (MQ=255)
                                                                gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGc                                                                                   >  3:1435828/1‑149 (MQ=255)
                                                                   gAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATc                                                                                >  2:460069/1‑149 (MQ=255)
                                                                     aaGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAc                                                                              >  2:281586/1‑149 (MQ=255)
                                                                     aaGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAc                                                                              <  4:177335/149‑1 (MQ=255)
                                                                      aGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAcc                                                                             <  2:632439/149‑1 (MQ=255)
                                                                       gggCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAccc                                                                            <  4:940002/149‑1 (MQ=255)
                                                                       gggCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCAccc                                                                            >  1:592392/1‑149 (MQ=255)
                                                                          cTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGtt                                                                         <  1:163503/149‑1 (MQ=255)
                                                                          cTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGtt                                                                         <  1:586475/149‑1 (MQ=255)
                                                                            ggAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGc                                                                       <  1:1144820/149‑1 (MQ=255)
                                                                               aaaTCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAac                                                                    >  2:206128/1‑149 (MQ=255)
                                                                                 aTCTCTCTGGCAATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGc                                                                                                          <  4:1159384/109‑1 (MQ=255)
                                                                                 aTCTCTCTAGCAATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGc                                                                                                          >  3:1159384/1‑109 (MQ=255)
                                                                                    tctctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTc                                                               >  3:976397/1‑149 (MQ=255)
                                                                                       ctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGc                                                            <  2:592392/149‑1 (MQ=255)
                                                                                        tGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAg                                                                                                            <  4:206237/100‑1 (MQ=255)
                                                                                        tGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAg                                                                                                            >  3:206237/1‑100 (MQ=255)
                                                                                        tGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGtt                                                                         >  4:204851/1‑135 (MQ=255)
                                                                                        tGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGtt                                                                         <  3:204851/135‑1 (MQ=255)
                                                                                          gCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTc                                                         <  3:289240/149‑1 (MQ=255)
                                                                                            gATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCa                                                       >  4:775714/1‑149 (MQ=255)
                                                                                              tCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGc                                                     >  4:508790/1‑149 (MQ=255)
                                                                                                tGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTc                                                               <  4:1144878/137‑1 (MQ=255)
                                                                                                tGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTc                                                               >  3:1144878/1‑137 (MQ=255)
                                                                                                     gagaAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTa                                              >  2:540292/1‑149 (MQ=255)
                                                                                                     gagaAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTGTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTa                                              >  4:878951/1‑149 (MQ=255)
                                                                                                          gaTAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGc                                         >  2:762615/1‑149 (MQ=255)
                                                                                                          gaTAACTAAAGCACTGGCTGGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTGACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGc                                         >  4:1258070/1‑149 (MQ=255)
                                                                                                           aTAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCa                                        <  4:435387/149‑1 (MQ=255)
                                                                                                           aTAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCa                                        <  2:39295/149‑1 (MQ=255)
                                                                                                             aaCTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCGCCGCTAATCGCAcc                                      <  4:1024779/149‑1 (MQ=255)
                                                                                                               cTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGa                                    >  3:340035/1‑149 (MQ=255)
                                                                                                                tAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGc                                                                       >  4:394016/1‑113 (MQ=255)
                                                                                                                tAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGc                                                                       <  3:394016/113‑1 (MQ=255)
                                                                                                                tAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGAc                                   >  3:641203/1‑149 (MQ=255)
                                                                                                                tAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGAc                                   >  2:163605/1‑149 (MQ=255)
                                                                                                                tAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGAc                                   >  1:40648/1‑149 (MQ=255)
                                                                                                                  aaGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCATGCGAGTTTCAGATTCACCACGTTGCTTCACCCGTTGCAACACCTCCGCTTCGAGCTCCGCTAATCACACCGAccc                                 >  3:1252874/1‑149 (MQ=255)
                                                                                                                  aaGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGAccc                                 <  1:206128/149‑1 (MQ=255)
                                                                                                                  aaGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGAccc                                 >  4:656421/1‑149 (MQ=255)
                                                                                                                      aCTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAAt                             <  1:1067452/149‑1 (MQ=255)
                                                                                                                      aCTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAAt                             <  2:761506/149‑1 (MQ=255)
                                                                                                                          gCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGAc                         <  2:1437688/149‑1 (MQ=255)
                                                                                                                          gCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGAc                         >  4:1108501/1‑149 (MQ=255)
                                                                                                                            ttGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGAcgc                       >  3:43158/1‑149 (MQ=255)
                                                                                                                            ttGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGAcgc                       >  3:1113141/1‑149 (MQ=255)
                                                                                                                            ttGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGAcgc                       >  4:57863/1‑149 (MQ=255)
                                                                                                                            ttGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGAcgc                       >  4:643011/1‑149 (MQ=255)
                                                                                                                             tGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGAc                                   <  3:757003/136‑1 (MQ=255)
                                                                                                                                   ataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACggg                <  4:808422/149‑1 (MQ=255)
                                                                                                                                   ataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACggg                >  1:1260221/1‑149 (MQ=255)
                                                                                                                                   ataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACggg                >  1:914618/1‑149 (MQ=255)
                                                                                                                                      aCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGa             >  1:422966/1‑149 (MQ=255)
                                                                                                                                       ccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGat            >  3:421558/1‑149 (MQ=255)
                                                                                                                                       ccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGat            >  3:583722/1‑149 (MQ=255)
                                                                                                                                       ccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGat            >  4:787327/1‑149 (MQ=255)
                                                                                                                                         gAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATc                                           <  1:1214841/116‑1 (MQ=255)
                                                                                                                                         gAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATc                                           >  2:1214841/1‑116 (MQ=255)
                                                                                                                                         gAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGatat          >  3:181314/1‑149 (MQ=255)
                                                                                                                                          aaTGGCGGCAATAGCGCCGCCATTCGGTGAATTTACCCCGTTTtctc                                                                                                              <  2:410906/47‑1 (MQ=255)
                                                                                                                                              gcggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACt     >  4:385686/1‑149 (MQ=255)
                                                                                                                                               cggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTg    >  1:1146485/1‑149 (MQ=255)
                                                                                                                                                ggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGt   <  4:1318032/149‑1 (MQ=255)
                                                                                                                                                 gcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTc  <  3:878951/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTC  >  NC_000913/993107‑993400

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