Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,991,798 Δ1 bp 5.0% intergenic (+39/‑296) ygeG → / → ygeH SycD‑like chaperone family TPR‑repeat‑containing protein/putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,991,7980C. 5.0% 273.3 / 9.5 99intergenic (+39/‑296)ygeG/ygeHSycD‑like chaperone family TPR‑repeat‑containing protein/putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (44/49);  new base . (3/2);  total (47/52)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-02

CAATGCAATGTTTCTATCACATTATTAACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTATTTCAGATTACATCAC  >  NC_000913/2991650‑2991946
                                                                                                                                                    |                                                                                                                                                      
caatgcaatgTTTCTATCACATTATTAACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCa                                                                                                                                                       >  2:542917/1‑150 (MQ=255)
 aatgcaatgTTTCTATCACATTATTAACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAg                                                                                                                                                      <  2:1172648/150‑1 (MQ=255)
  atgcaatgTTTCTATCACATTATTAACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGt                                                                                                                                                     <  1:735002/150‑1 (MQ=255)
   tgcaatgTTTCTATCACATTATTAACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGtt                                                                                                                                                    >  2:613002/1‑150 (MQ=255)
    tcaatgTTTCTATCACATTATTAACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTa                                                                                                                                                   <  1:505177/149‑1 (MQ=255)
     caatgTTTCTATCACATTATTAACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGAATATATTGAACTCTTAAACGATAATTCAGAAGAAAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTaa                                                                                                                                                  <  2:1087595/150‑1 (MQ=255)
      aatgTTTCTATCACATTATTAACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAc                                                                                                                                                 >  2:529097/1‑150 (MQ=255)
       atgTTTCTATAACATTATTAACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAc‑‑c                                                                                                                                              <  1:841358/150‑1 (MQ=255)
        tgTTTCTATCACTATATTAACAATTGTGAGGAAGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑Ca                                                                                                                                             <  2:812901/150‑1 (MQ=255)
             cTATCACATTATTAACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACtata                                                                                                                                        <  2:1019463/150‑1 (MQ=255)
              tATCACATTATTAACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATaa                                                                                                                                       <  1:1264554/150‑1 (MQ=255)
              tATCACATTATTAACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATaa                                                                                                                                       <  1:119805/150‑1 (MQ=255)
                tCACATTATTAACATTTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATTGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAAcc                                                                                                                                     <  1:893171/150‑1 (MQ=255)
                tCACATTATTAACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTTAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAAcc                                                                                                                                     <  2:637830/150‑1 (MQ=255)
                tCACATTATTAACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTG‑‑GTTAACTGCACTATAAcc                                                                                                                                     <  2:944324/150‑1 (MQ=255)
                   aattattAACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTG‑‑GTTAACTGCACTATAACCAGc                                                                                                                                  <  1:1238952/149‑1 (MQ=255)
                       attaACAATTGTGAGGATGATAGTGTTAAGAGTAAAGATCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACct                                                                                                                              >  2:1258496/1‑150 (MQ=255)
                         taACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACctct                                                                                                                            <  1:825424/150‑1 (MQ=255)
                         taACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACctct                                                                                                                            >  1:287218/1‑150 (MQ=255)
                          aaCAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGAAAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTAGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCtt                                                                                                                           >  1:423060/1‑150 (MQ=255)
                           aCAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTa                                                                                                                          >  2:950593/1‑150 (MQ=255)
                              aTTGTGAGGATGATAGTGTTTAGAGTAAAGGGCAGGCATATATTGAACTCTTAAACGATAATTCTGATGATTATGGCTAAGAATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGt                                                                                                                       <  2:738534/150‑1 (MQ=255)
                              aTTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGt                                                                                                                       >  1:292894/1‑150 (MQ=255)
                               ttGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGtt                                                                                                                      <  1:322828/150‑1 (MQ=255)
                                    aGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTTAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCtt                                                                                                                 >  2:1149314/1‑150 (MQ=255)
                                    aGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCtt                                                                                                                 <  1:747201/150‑1 (MQ=255)
                                    aGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCtt                                                                                                                 <  2:599831/150‑1 (MQ=255)
                                     ggATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTa                                                                                                                >  1:620967/1‑150 (MQ=255)
                                      gatgatAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTaa                                                                                                               <  1:876756/150‑1 (MQ=255)
                                        tgatAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAAt                                                                                                             >  2:207106/1‑150 (MQ=255)
                                          atAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCa                                                                                                           <  2:567725/150‑1 (MQ=255)
                                            aGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGc                                                                                                         >  2:133741/1‑150 (MQ=255)
                                             ttgtTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCa                                                                                                        <  1:1183590/149‑1 (MQ=255)
                                               gtTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATc                                                                                                      >  1:412843/1‑150 (MQ=255)
                                                 tAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAAATCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAG‑‑CACTATAAACAGCACCTCTTATGTTTTCTTAAATCAGCATCtt                                                                                                    >  2:457733/1‑150 (MQ=255)
                                                   agagTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTaa                                                                                                  >  1:689918/1‑150 (MQ=255)
                                                    gagTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTaaa                                                                                                 <  1:613002/150‑1 (MQ=255)
                                                      gTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTaaaaa                                                                                               >  1:395733/1‑149 (MQ=255)
                                                         aaGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGt                                                                                            <  2:1249724/150‑1 (MQ=255)
                                                           gcgcAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTa                                                                                          >  1:778799/1‑150 (MQ=255)
                                                           gcgcAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑AACTATAACCAGCACCTCTTATGTTTTATTAAATCAGCATCTTAAAACAGTTa                                                                                          >  2:1146285/1‑150 (MQ=255)
                                                             gcAGGCATATATTGAACTCTTAAACGATTATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTTGTTATAGTGCAGTTAAC‑‑CACTATAACCAGAACATCTTATGTTTTCTTATATCAGAATCTTAAAACAGTTAag                                                                                        >  1:191227/1‑150 (MQ=255)
                                                              cAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAaga                                                                                       >  1:482745/1‑150 (MQ=255)
                                                                  cATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCAt                                                                                   <  2:79666/150‑1 (MQ=255)
                                                                          gAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATaaa                                                                           >  2:20397/1‑150 (MQ=255)
                                                                             ctctTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGTTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACATCTAATGTTTTATTAAATCAGCATCTTAAAACAGTAAATAGAAAAATATAAAtaa                                                                        >  2:1292850/1‑148 (MQ=255)
                                                                             ctctTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGATAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCAACTATTATGTTATCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACAATAAATCa                                                                        >  1:973244/1‑150 (MQ=255)
                                                                              tctTAAAAGATAATTCAGAAGATAATTTATAAGCATTACAATAATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCat                                                                       <  1:21954/150‑1 (MQ=255)
                                                                                ttAATCGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCatta                                                                     <  1:1110133/150‑1 (MQ=255)
                                                                                 tAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTTGTTAAATTTGCTGTTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCattat                                                                    <  1:199329/150‑1 (MQ=255)
                                                                                 tAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAACTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCattat                                                                    >  2:692467/1‑150 (MQ=255)
                                                                                  aaaCGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAAAAGTTAAGAGCATTCTATAAATCattatt                                                                   >  1:628908/1‑150 (MQ=255)
                                                                                   aaCGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTAttt                                                                  <  2:1316787/150‑1 (MQ=255)
                                                                                       aTAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCa                                                              >  1:258801/1‑150 (MQ=255)
                                                                                       aTAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCa                                                              <  2:491754/150‑1 (MQ=255)
                                                                                          aTTCATAAGATAATGTCTAATCATTACAATCATTTTTTTAAATTTTCTTGTTATAGTTCAGTTAAA‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATg                                                           <  1:183220/150‑1 (MQ=255)
                                                                                          aTTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAATAGCATACTATAAATCATTATTTTTCAATg                                                           >  1:161503/1‑150 (MQ=255)
                                                                                            tCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCg                                                         <  2:882460/150‑1 (MQ=255)
                                                                                              agaagaTAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGtt                                                       >  2:1296924/1‑150 (MQ=255)
                                                                                               gaagaTAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATATTTCAATGCGTTa                                                      >  2:625302/1‑150 (MQ=255)
                                                                                               gaagaTAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGAATACTATAAATCATTATTTTTCAATGCGTTa                                                      >  2:489420/1‑150 (MQ=255)
                                                                                                  gaTAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTagag                                                   <  2:863847/150‑1 (MQ=255)
                                                                                                   aTAAAGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTTGAGGTAAC‑‑CACTATCACCAGGACCTCTTATGTTTTCTTTTATCTGCTTCTTAATACAGTTAAGAGCATACTATAAATCATTATTTTTCATTGCGTTAgaga                                                  <  1:666528/150‑1 (MQ=255)
                                                                                                       tGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGAGTTAGAGAaaaa                                              >  2:607451/1‑150 (MQ=255)
                                                                                                        ggCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAAt                                             <  2:1194333/150‑1 (MQ=255)
                                                                                                            aaGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTAc                                         >  2:258247/1‑150 (MQ=255)
                                                                                                              ttaTTACATTCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑AACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACtt                                       <  1:116113/148‑1 (MQ=255)
                                                                                                                       tCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAAt                              <  2:620967/150‑1 (MQ=255)
                                                                                                                       tCATTTGGTTAAATTTGATGGTTATAGTTCAGTTAAC‑‑CACTATAACCAGAAACTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAAt                              <  1:319843/150‑1 (MQ=255)
                                                                                                                         aTTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTg                            >  1:970378/1‑150 (MQ=255)
                                                                                                                         aTTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTg                            <  1:950593/150‑1 (MQ=255)
                                                                                                                           ttGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGaa                          >  2:669422/1‑150 (MQ=255)
                                                                                                                           ttGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGaa                          <  1:1199690/150‑1 (MQ=255)
                                                                                                                             ggTTAAATTTGCTTGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATa                        <  1:963331/150‑1 (MQ=255)
                                                                                                                             ggTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAATCAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATa                        >  1:543270/1‑150 (MQ=255)
                                                                                                                                tAAATTTGCTGTTTATAGTGCAGTTAAC‑‑CACTATTACCAGAACCTCTTATTTTTTCTTAAATCAGCATCTTAATATAGTTATGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAgg                     <  1:632463/150‑1 (MQ=255)
                                                                                                                                 aaaTTTGCTTGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGc                    <  2:809671/150‑1 (MQ=255)
                                                                                                                                  aaTTTGCTTGTTATAGTGCAGTTAAC‑‑CACTATAAACAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCt                   <  1:1268002/150‑1 (MQ=255)
                                                                                                                                     ttGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTAtt                <  1:657256/150‑1 (MQ=255)
                                                                                                                                     ttGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGGTTTCTTAAATCAGCATCTTAAAACAGTTAACAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTAtt                <  2:1143345/150‑1 (MQ=255)
                                                                                                                                      tGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTAttt               >  2:282353/1‑150 (MQ=255)
                                                                                                                                      tGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGATTTCTTAAATCAGCATCTTAAAACAGGTAAGAGCATACTATAAATCATTATTTTTCAATGCGCTAGAGAACAAACTACTTTATCAAAAATGAATAAGGCTAttt               <  2:41070/150‑1 (MQ=255)
                                                                                                                                      tGCTGGTTATAGTG‑‑GTTAACTGCACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTAttt               >  2:1214072/1‑150 (MQ=39)
                                                                                                                                      tGCTGGTTATAGTG‑‑GTTAACTGCACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTAttt               >  2:379854/1‑150 (MQ=39)
                                                                                                                                       gCTGGTTATAGTGCAGTTAAA‑‑CAATATAACCAGAAACTATTATGTTTTCTTAAATCATCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTATTTc              <  1:216432/150‑1 (MQ=255)
                                                                                                                                        cTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTATTTCa             <  2:1032225/150‑1 (MQ=255)
                                                                                                                                        cTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTATTTCa             <  2:1091911/150‑1 (MQ=255)
                                                                                                                                        cTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTATTTCa             >  2:475029/1‑150 (MQ=255)
                                                                                                                                        cTGGTTATAGTG‑‑GTTAACTGCACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTATTTTATCAAAATTGAATAAGGCTATTTCa             >  1:311061/1‑150 (MQ=255)
                                                                                                                                         tGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTATTTCAg            <  1:1292293/150‑1 (MQ=255)
                                                                                                                                          ggTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTATTTCAGa           >  1:917386/1‑150 (MQ=255)
                                                                                                                                              ataGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTATTTCAGATTAc       <  1:1133732/150‑1 (MQ=255)
                                                                                                                                                aGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTATTTCAGATTACAt     >  1:598523/1‑150 (MQ=255)
                                                                                                                                                aGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTATTTCAGATTACAt     >  1:1027435/1‑150 (MQ=255)
                                                                                                                                                 gTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACGATCTACTTTATCAAAATTGAATAAGGCTATTTCAGATTACAta    >  2:753780/1‑149 (MQ=255)
                                                                                                                                                 gTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTATTTCAGATTACAta    >  1:468381/1‑149 (MQ=255)
                                                                                                                                                 gTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTATTTCAGATTACATc    <  1:717829/150‑1 (MQ=255)
                                                                                                                                                 gTGCAGTTAAC‑‑CACTATAACCAGCACCTCGTATGTTTTCTTAACTCAGCATCTTAAAACTGGTAAGTGCATACTATAACTGATTATTTTTCAATGCGTTAGAGAACAATCTAGTTTATCAAAATTGAATAAGGCTATTTCAGTTTTCATc    >  2:812532/1‑150 (MQ=255)
                                                                                                                                                   gCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTATTTCAGATTACATCAc  >  1:178397/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CAATGCAATGTTTCTATCACATTATTAACAATTGTGAGGATGATAGTGTTAAGAGTAAAGCGCAGGCATATATTGAACTCTTAAACGATAATTCAGAAGATAATGGCTAAGCATTACAATCATTTGGTTAAATTTGCTGGTTATAGTGCAGTTAAC‑‑CACTATAACCAGCACCTCTTATGTTTTCTTAAATCAGCATCTTAAAACAGTTAAGAGCATACTATAAATCATTATTTTTCAATGCGTTAGAGAACAATCTACTTTATCAAAATTGAATAAGGCTATTTCAGATTACATCAC  >  NC_000913/2991650‑2991946

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: