Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,146,763 T→C 8.4% intergenic (‑26/+93) yghW ← / ← yghX DUF2623 family protein/pseudogene, dienlactone hydrolase family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,146,7630TC 8.4% 228.2 / 9.3 83intergenic (‑26/+93)yghW/yghXDUF2623 family protein/pseudogene, dienlactone hydrolase family
Reads supporting (aligned to +/- strand):  ref base T (33/43);  new base C (2/5);  total (35/48)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-01

GTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGG  >  NC_000913/3146615‑3146911
                                                                                                                                                    |                                                                                                                                                    
gtgtGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGt                                                                                                                                                      >  2:445472/1‑149 (MQ=255)
gtgtGAGTAGCCTAATACAAAACCGCGATTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGt                                                                                                                                                      >  1:119732/1‑149 (MQ=255)
 tgtgAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTc                                                                                                                                                     <  2:307563/149‑1 (MQ=255)
  gtgAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCa                                                                                                                                                    <  1:185829/149‑1 (MQ=255)
      ttAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATTCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCatata                                                                                                                                                <  2:884637/148‑1 (MQ=255)
       tAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCatatat                                                                                                                                               >  2:605777/1‑149 (MQ=255)
         gCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAg                                                                                                                                             <  1:188464/149‑1 (MQ=255)
          ccTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATGACAATATATAGc                                                                                                                                            <  1:382189/149‑1 (MQ=255)
          ccTAATACAAAAACGCGTGTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATGACAATATATAGc                                                                                                                                            <  1:871961/149‑1 (MQ=255)
            tAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCATATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATGACAATATATAGCGa                                                                                                                                          <  1:1024893/149‑1 (MQ=255)
            tAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCATATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATGACAATATATAGCGa                                                                                                                                          <  1:1029579/149‑1 (MQ=255)
              aTACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGata                                                                                                                                        >  1:426224/1‑149 (MQ=255)
                aCAAAACCGCGATTATAATCGGCACAGAATTTTGTCACATTAACCCCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTATGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATAtt                                                                                                                                      <  2:1149270/149‑1 (MQ=255)
                   aaaCCTCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCTTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTa                                                                                                                                   <  2:464689/149‑1 (MQ=255)
                   aaaCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTa                                                                                                                                   <  1:962497/149‑1 (MQ=255)
                   aaaCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATGACAATATATAGCGaatgaggt                                                                                                                                   <  1:152606/149‑8 (MQ=255)
                   aaaCCGCGTTTATAATCGGCACAGAATTTTGTAACATTAACCGCACTTTAGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTa                                                                                                                                   <  2:720486/149‑1 (MQ=255)
                     aCCGCGGTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGGCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAgg                                                                                                                                 <  2:556871/149‑1 (MQ=255)
                      ccGCGTTTATAATCGGCACAGAATTTTGGCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGc                                                                                                                                <  2:438139/149‑1 (MQ=255)
                       cgcgTTTATAATCGGCACAGAATTTTGTCACAATAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCa                                                                                                                               <  2:1047757/149‑1 (MQ=255)
                           tttATAATCGGAACAGAATTTTGTCACATTAACCGCACTTTCGCCATGCGTTGCTTTTAATCCCGCAATTAAGACTTTAAAAAAATGGTTATTCATCGTTCACCTCATTCTCTATATATTGTCATATATATCGATATTATAGGCaaaaa                                                                                                                           >  2:883486/1‑147 (MQ=255)
                             tataATCTGCACATAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAAc                                                                                                                         <  1:288064/149‑1 (MQ=255)
                               taATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGa                                                                                                                       <  1:868856/149‑1 (MQ=255)
                                aaTCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGAc                                                                                                                      >  2:709564/1‑149 (MQ=255)
                                   cGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACttt                                                                                                                   <  1:294107/149‑1 (MQ=255)
                                   cGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACttt                                                                                                                   >  1:736576/1‑149 (MQ=255)
                                    ggCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACtttt                                                                                                                  <  1:393646/149‑1 (MQ=255)
                                     gCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTg                                                                                                                 >  2:36145/1‑149 (MQ=255)
                                     gCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTg                                                                                                                 >  1:379620/1‑149 (MQ=255)
                                      cacaGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTgc                                                                                                                >  2:433039/1‑149 (MQ=255)
                                         aGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCa                                                                                                             >  2:656076/1‑149 (MQ=255)
                                         aGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCa                                                                                                             >  2:526770/1‑149 (MQ=255)
                                             ttttGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATAAAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCt                                                                                                         >  2:592965/1‑149 (MQ=255)
                                                     caTTAACCGCACTGTCGGCATTCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTc                                                                                                 <  1:386334/149‑1 (MQ=255)
                                                        tAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCCTTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATg                                                                                              <  2:968173/149‑1 (MQ=255)
                                                         aaCCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGc                                                                                             <  2:151597/149‑1 (MQ=255)
                                                         aaCCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGc                                                                                             <  2:778775/149‑1 (MQ=255)
                                                          aCCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATATCGATATTTTAGGCAAAGAACGACTTTTTCGCACGCTATTTTTTCATGCt                                                                                            >  2:995560/1‑149 (MQ=255)
                                                           ccGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGCTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGAGATATTTTAGGCTAAGAAAGACTTTTGCGCACGCTTTTTTTTCATGCtt                                                                                           >  1:274429/1‑149 (MQ=255)
                                                               aCTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGttt                                                                                       <  2:118353/149‑1 (MQ=255)
                                                                   tCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGa                                                                                   <  1:283019/149‑1 (MQ=255)
                                                                   tCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTAACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGAGCACGCTATTTTTTCATGCTTGTTTAAGa                                                                                   >  1:703772/1‑149 (MQ=255)
                                                                     ggCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACg                                                                                 <  1:325414/149‑1 (MQ=255)
                                                                        aTGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTaa                                                                              >  1:1163874/1‑149 (MQ=255)
                                                                         tGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAAt                                                                             >  2:206790/1‑149 (MQ=255)
                                                                          gCGATGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATc                                                                            <  2:218690/149‑1 (MQ=255)
                                                                           cGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATcc                                                                           <  1:656076/149‑1 (MQ=255)
                                                                            gTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATccc                                                                          >  1:1023210/1‑149 (MQ=255)
                                                                              tGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGt                                                                        >  1:528480/1‑149 (MQ=255)
                                                                                cTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGCACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTaa                                                                      >  2:615284/1‑149 (MQ=255)
                                                                                       tCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATaa                                                               >  1:624046/1‑149 (MQ=255)
                                                                                         ccGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAAc                                                             <  1:484180/149‑1 (MQ=255)
                                                                                         ccGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAAc                                                             <  2:525111/149‑1 (MQ=255)
                                                                                            ccATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTa                                                          <  1:709564/149‑1 (MQ=255)
                                                                                                  aGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGga                                                    <  1:445472/149‑1 (MQ=255)
                                                                                                  aGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGga                                                    <  2:929130/149‑1 (MQ=255)
                                                                                                  aGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGga                                                    >  1:602152/1‑149 (MQ=255)
                                                                                                            aaaaTGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTa                                          >  1:1041382/1‑149 (MQ=255)
                                                                                                             aaaTGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTAt                                         >  2:466595/1‑149 (MQ=255)
                                                                                                             aaaTGGTTATTCATAGTTCAACTAATTCGCTATATATTGTAATATATAGCGATATTTTATGCAAAGAACGACTTTTGCGCACGCTATTTTTTCCTGCTTGTTTAAGACGTAATCCCGTAAGATATAAACCTATCAGGAGAGATATTTAt                                         <  1:592436/149‑1 (MQ=255)
                                                                                                                tGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGa                                      <  1:526770/149‑1 (MQ=255)
                                                                                                                  gTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAAc                                    >  1:774595/1‑149 (MQ=255)
                                                                                                                      ttCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATACACCTATCAGGAGAGATATTTATCGAACCAtt                                >  1:785118/1‑149 (MQ=255)
                                                                                                                       tCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCACGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAAACAttt                               >  1:20014/1‑149 (MQ=255)
                                                                                                                         ttCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCATAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCa                             <  2:734818/148‑1 (MQ=255)
                                                                                                                         aTCGTTCACCTCATTCTGTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCAAGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGCTAAACCTATCAGGAGAGATATTTATCGAACCATTTCa                             <  2:217022/149‑1 (MQ=255)
                                                                                                                         aTCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCa                             <  2:1142590/149‑1 (MQ=255)
                                                                                                                         aTCGTTCACCTCATTCGCTATATATTGTAATATATAGCGATATTATAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCa                             <  1:214574/149‑1 (MQ=255)
                                                                                                                          tCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAg                            <  2:426224/149‑1 (MQ=255)
                                                                                                                            gTTCACCTCATTCGATATATATTGTCATATATAGAGATATTGATGGCACAGAACGACTTTTGCGCACGCTATATTTTCATGCTTGTTTAAGCCGTAATCCCGTACGAGATAAACCTATCAGGCGAGATATTTATCGAACCATTTCAgtg                          <  1:989584/149‑1 (MQ=255)
                                                                                                                             ttaaCCTCATTCGCTATATATTGTCATATATAGCGAGATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGAAATAGAAAAACCTATCAGGAGAGATATTTATCGAACCATTTCAgtgt                         >  2:1083670/4‑149 (MQ=255)
                                                                                                                             ttCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAgtgt                         <  2:379620/149‑1 (MQ=255)
                                                                                                                             ttCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTAGTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAgtgt                         >  2:510139/1‑149 (MQ=255)
                                                                                                                               cACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTcc                       <  1:433039/149‑1 (MQ=255)
                                                                                                                                aCCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCt                      >  1:524301/1‑149 (MQ=255)
                                                                                                                                 ccgcATTCGCTTTATATTGTCATATATAGCGATATTTTTGTAAAAGAACTACTTTTGCGCACGCTATTTATTCTTGCTTGTTTAAGACGTATTCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGGGTCGtt                     <  2:568843/146‑1 (MQ=255)
                                                                                                                                  ctctTTCGCTATATTTTGTCATATATAGCGATATTTTAGGCAAAGTACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCttt                    <  2:119732/149‑1 (MQ=255)
                                                                                                                                  cTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCttt                    >  1:838363/1‑149 (MQ=255)
                                                                                                                                     aTTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGcc                 >  1:526158/1‑149 (MQ=255)
                                                                                                                                         gCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGc             >  1:664981/1‑149 (MQ=255)
                                                                                                                                         gCTATATATGACAATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGc             >  1:458776/1‑149 (MQ=255)
                                                                                                                                         gCTATATATGACAATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGc             >  1:157269/1‑149 (MQ=255)
                                                                                                                                            atatatTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTAAGTGTCCTTTGCCAGGCAAg          >  2:240598/1‑149 (MQ=255)
                                                                                                                                                atTGTCATATATAGCGATATTTTAGGAAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCACGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATcg      <  2:74614/149‑1 (MQ=255)
                                                                                                                                                   gTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCCCGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATcggcg   <  2:129927/149‑1 (MQ=255)
                                                                                                                                                    tCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATcggcgg  >  1:497504/1‑149 (MQ=255)
                                                                                                                                                    tCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATcggcgg  <  2:607104/149‑1 (MQ=255)
                                                                                                                                                    tCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATcggcgg  <  2:115378/149‑1 (MQ=255)
                                                                                                                                                    tCATATATAGCGATATTTTAGGCAAAGAACGACTATTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCATGCAAGAAcggcgg  >  1:900097/1‑149 (MQ=255)
                                                                                                                                                    tCATATATAGCGATAATTTAGGCAAAGAACGACTTTTGGGCACGCTATTTTATCAAGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATcggcgg  <  2:246126/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGG  >  NC_000913/3146615‑3146911

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: