Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,146,762 G→A 8.3% intergenic (‑25/+94) 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,7620GA 8.3% 243.3 / 9.6 84intergenic (‑25/+94)yghW/yghXDUF2623 family protein/pseudogene, dienlactone hydrolase family
Reads supporting (aligned to +/- strand):  ref base G (35/42);  new base A (2/5);  total (37/47)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCG  >  NC_000913/3146615‑3146910
                                                                                                                                                   |                                                                                                                                                    
gtgtGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGt                                                                                                                                                     >  2:445472/1‑149 (MQ=255)
gtgtGAGTAGCCTAATACAAAACCGCGATTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGt                                                                                                                                                     >  1:119732/1‑149 (MQ=255)
gtggGAGTAGCCTAATACAAAACCGCTGTTATCATCGGCTCAGAATTTTGTCACATTATCCGCACGGTCGGCATGCGTTGCTTTTAATCCCGCCTATACGCCTTTACAATAATGGTTATTCATCGTTCACCTCATTCGCTCTATATTgg                                                                                                                                                     <  2:774093/149‑2 (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                                                                                                                >  1:379620/1‑149 (MQ=255)
                                     gCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTg                                                                                                                >  2:36145/1‑149 (MQ=255)
                                      cacaGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTgc                                                                                                               >  2:433039/1‑149 (MQ=255)
                                         aGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCa                                                                                                            >  2:526770/1‑149 (MQ=255)
                                         aGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCa                                                                                                            >  2:656076/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:778775/149‑1 (MQ=255)
                                                         aaCCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGc                                                                                            <  2:151597/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                                                   <  2:929130/149‑1 (MQ=255)
                                                                                                  aGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGga                                                   >  1:602152/1‑149 (MQ=255)
                                                                                                  aGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGga                                                   <  1:445472/149‑1 (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)
                                                                                                                                                   |                                                                                                                                                    
GTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCG  >  NC_000913/3146615‑3146910

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: