Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,146,764 C→A 8.5% intergenic (‑27/+92) 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,7640CA 8.5% 362.1 / 16.8 127intergenic (‑27/+92)yghW/yghXDUF2623 family protein/pseudogene, dienlactone hydrolase family
Reads supporting (aligned to +/- strand):  ref base C (49/67);  new base A (6/5);  total (55/72)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

Only 100 of 128 total aligned reads displayed.

GTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGCGG  >  NC_000913/3146615‑3146911
                                                                                                                                                     |                                                                                                                                                   
gtgtGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTc                                                                                                                                                     <  1:798759/150‑1 (MQ=255)
 tgtgAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCa                                                                                                                                                    <  1:628013/150‑1 (MQ=255)
  gtgAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCat                                                                                                                                                   <  1:920285/150‑1 (MQ=255)
   tgAGTAGCCTAATACAAAACCGCGTTTATAATCGGCTCAGAATTTTGTCACATTAACCTCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAATATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCata                                                                                                                                                  <  2:735085/150‑1 (MQ=255)
   tgAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCata                                                                                                                                                  <  1:1101872/150‑1 (MQ=255)
    gAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTAACCTCATTCGCTATATATTGTCatat                                                                                                                                                 >  2:444177/1‑150 (MQ=255)
       tAGCCTAATACAAAACGGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCAtatata                                                                                                                                              <  1:751276/150‑1 (MQ=255)
        aGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAg                                                                                                                                             >  2:1132630/1‑150 (MQ=255)
        aGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGGCATATATAg                                                                                                                                             <  1:1535303/150‑1 (MQ=255)
         gCCTAATTCAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGc                                                                                                                                            <  1:192572/150‑1 (MQ=255)
         gCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGc                                                                                                                                            <  1:472474/150‑1 (MQ=255)
         gCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGc                                                                                                                                            <  1:1369142/150‑1 (MQ=255)
         gCCTAATACAAAACCGCGCTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGc                                                                                                                                            <  1:255619/150‑1 (MQ=255)
          ccTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGTCATTCGTTGCTTTTAATCCCGCCATTAAGACTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCg                                                                                                                                           <  2:1246727/150‑1 (MQ=255)
          ccTAATACAAAAACGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCTCACTGTCGGCTTGCGTTGCTTTTAATCACGCCATTAAGCCTTTACCAAAAAGGTTATTCATCGTTCACCTCATTCGCTATAAATGACAATATATAGCg                                                                                                                                           <  2:1041789/150‑1 (MQ=255)
           cTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATGACAATATATAGCGa                                                                                                                                          <  2:1220003/150‑1 (MQ=255)
            tAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGat                                                                                                                                         >  1:633421/1‑150 (MQ=255)
             aaTACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGata                                                                                                                                        <  1:1355348/150‑1 (MQ=255)
              aTACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGatat                                                                                                                                       <  1:1499654/150‑1 (MQ=255)
              aTACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGatat                                                                                                                                       >  2:808844/1‑150 (MQ=255)
               tACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATAtt                                                                                                                                      >  2:1152115/1‑150 (MQ=255)
                 cAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATAtttt                                                                                                                                    >  2:1027893/1‑150 (MQ=255)
                  aaaaCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTa                                                                                                                                   >  1:1534287/1‑150 (MQ=255)
                   aaaCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAg                                                                                                                                  >  1:823939/1‑150 (MQ=255)
                    aaCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAgg                                                                                                                                 <  2:160772/150‑1 (MQ=255)
                    aaCCGCGTTTATAATCGGCACAGAATTTTGTCAAATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAgg                                                                                                                                 <  2:23271/150‑1 (MQ=255)
                          gTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGa                                                                                                                           <  2:629846/150‑1 (MQ=255)
                             tataATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACg                                                                                                                        <  1:1200101/150‑1 (MQ=255)
                               taATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCTCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGAc                                                                                                                      >  1:313757/1‑150 (MQ=255)
                               taATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGATATTCATCGTTCGCCTCATACGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGAc                                                                                                                      >  1:1156353/1‑150 (MQ=255)
                                  tCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACttt                                                                                                                   <  2:921270/150‑1 (MQ=255)
                                    ggCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTg                                                                                                                 <  2:1685197/150‑1 (MQ=255)
                                    ggCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTg                                                                                                                 >  2:799257/1‑150 (MQ=255)
                                     gCACAGAATTTTGTCACATTAACCTCACTTTAGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTgc                                                                                                                <  1:1215790/150‑1 (MQ=255)
                                      cacaGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTgcg                                                                                                               >  1:374337/1‑150 (MQ=255)
                                      cacaGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTgcg                                                                                                               >  1:135333/1‑150 (MQ=255)
                                          gAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTAAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGAGATTTTAGGCAAAGAACGACTTTTGCGCACg                                                                                                           >  1:261878/1‑150 (MQ=255)
                                           aaTTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGc                                                                                                          >  2:928452/1‑150 (MQ=255)
                                                  tCACATTAACCGCACTGTCGGAATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTAttttt                                                                                                   <  2:135901/150‑1 (MQ=255)
                                                   cacaTTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTAtttttt                                                                                                  <  1:796887/150‑1 (MQ=255)
                                                        tAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGc                                                                                             >  2:435365/1‑150 (MQ=255)
                                                          aCCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCtt                                                                                           >  1:555575/1‑150 (MQ=255)
                                                            cGCACTGTCGGCATGCGATGCTTTTAATCCCGCCATTAAGCCTTTACAAAAATGGTTATTCATCGTTCGCCTCATTCGCTATATATTGTCATATATAGCGAGATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCAAGCTTGt                                                                                         >  2:1194830/1‑150 (MQ=255)
                                                              cACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGttt                                                                                       <  2:675531/150‑1 (MQ=255)
                                                              cACTGTCGGCATGAGTTGCTTTTAATCACGCCTTTAAGCCTTTACCAATATGGTTATTCATCGTTCACCTCATTCGCTATATATTTTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGttt                                                                                       <  1:106095/150‑1 (MQ=255)
                                                                  gTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGTTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGa                                                                                   <  2:1315760/150‑1 (MQ=255)
                                                                  gTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGa                                                                                   >  2:934409/1‑150 (MQ=255)
                                                                    cGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGTCTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACg                                                                                 <  1:444177/150‑1 (MQ=255)
                                                                     ggCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGt                                                                                <  2:381404/150‑1 (MQ=255)
                                                                      gCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTa                                                                               >  2:439584/1‑150 (MQ=255)
                                                                          tcGTTGCTTTTAATCACGCCATTAAGCCTTTACGAAAATGTTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATcc                                                                           <  1:73320/149‑1 (MQ=255)
                                                                            gTTTCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATAATGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCg                                                                         <  1:461953/150‑1 (MQ=255)
                                                                             ttGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGt                                                                        <  2:374337/150‑1 (MQ=255)
                                                                             ttGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGt                                                                        <  1:808844/150‑1 (MQ=255)
                                                                              tGCTTTTAATCCAGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTa                                                                       <  1:1152115/150‑1 (MQ=255)
                                                                               gCTTTTAATCCCGCCATTACGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGACATATATAGCGATATTTTAGGCTAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTaa                                                                      >  2:858695/1‑150 (MQ=255)
                                                                               gCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTaa                                                                      >  1:1727293/1‑150 (MQ=255)
                                                                                cTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAag                                                                     <  2:342206/150‑1 (MQ=255)
                                                                                 ttttAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAaga                                                                    <  1:945787/150‑1 (MQ=255)
                                                                                 ttttAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAaga                                                                    <  1:676183/150‑1 (MQ=255)
                                                                                  tttAATCCCGCCATTAATCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAagag                                                                   <  1:1543496/150‑1 (MQ=255)
                                                                                   ttAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAgaga                                                                  <  1:411689/150‑1 (MQ=255)
                                                                                     aaTCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATa                                                                >  1:669124/1‑150 (MQ=255)
                                                                                       tCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATaaa                                                              <  1:395093/150‑1 (MQ=255)
                                                                                         ccGCCATTAAGCCTTTACCATAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAAcc                                                            <  2:1053825/150‑1 (MQ=255)
                                                                                             cATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTTCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATc                                                        >  1:105747/1‑150 (MQ=255)
                                                                                             cATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATc                                                        >  2:649412/1‑150 (MQ=255)
                                                                                                     cTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAATAGATAAACCTATCAGgagaga                                                <  2:897829/150‑1 (MQ=255)
                                                                                                         aCCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCTCTCTATTTTTTCATGCTTGTTTAAGACGTTATCCCGTAAGAGATAAACCTATCAGGTGAGATAtt                                            <  2:1186522/150‑1 (MQ=255)
                                                                                                         aCCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCAAGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATAtt                                            <  1:1706572/150‑1 (MQ=255)
                                                                                                          ccAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATAttt                                           >  1:693592/1‑150 (MQ=255)
                                                                                                            aaaaTGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTAt                                         <  1:19923/150‑1 (MQ=255)
                                                                                                            aaaaTGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTAt                                         >  2:44184/1‑150 (MQ=255)
                                                                                                            aaaaTGGTTATTCATCATTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTAt                                         >  2:520716/1‑150 (MQ=255)
                                                                                                             aaaTGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATc                                        <  1:1644637/150‑1 (MQ=255)
                                                                                                             aaaTGGTTATTCATCGTTCACCTCATGCGCTTTATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATc                                        <  1:919055/150‑1 (MQ=255)
                                                                                                               aTGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGa                                      >  2:736988/1‑150 (MQ=255)
                                                                                                               aTGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGa                                      >  2:50017/1‑150 (MQ=255)
                                                                                                                 ggTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAAc                                    <  1:1048991/150‑1 (MQ=255)
                                                                                                                  gTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAAcc                                   >  2:1027321/1‑150 (MQ=255)
                                                                                                                     aTTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCAtt                                <  1:799257/150‑1 (MQ=255)
                                                                                                                     aTTCATCGTTAACCTCATTCGCTATATATTGTAATATATAGCGATATTTTAGTAAAAGAACGACTTTTTCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCAtt                                <  1:81562/150‑1 (MQ=255)
                                                                                                                       tCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTc                              <  1:928452/150‑1 (MQ=255)
                                                                                                                               cACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCt                      >  2:509888/1‑150 (MQ=255)
                                                                                                                                   tCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTACTTTGc                  >  1:271051/1‑150 (MQ=255)
                                                                                                                                     aTTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCa                >  2:270233/1‑150 (MQ=255)
                                                                                                                                        cGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAgga             >  1:224803/1‑149 (MQ=255)
                                                                                                                                        cGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGc             <  2:669124/150‑1 (MQ=255)
                                                                                                                                         gCTATATATGACAATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAgg                                                     >  1:1697017/1‑109 (MQ=255)
                                                                                                                                         gCTATATATGACAATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAgg                                                     <  2:1697017/109‑1 (MQ=255)
                                                                                                                                         gCTATATATGACAATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGCACCATTTCAGTGTCCTTTGCCAGGCa            >  2:736855/1‑150 (MQ=255)
                                                                                                                                         gCTATATATGACAATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCa            >  2:1647096/1‑150 (MQ=255)
                                                                                                                                          cTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCaa           <  2:1465773/150‑1 (MQ=255)
                                                                                                                                          cTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCaa           <  1:375032/150‑1 (MQ=255)
                                                                                                                                               tatTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATcg      >  2:1126785/1‑150 (MQ=255)
                                                                                                                                                atTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTc                              <  1:48973/125‑1 (MQ=255)
                                                                                                                                                 ttGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATcggc    <  1:430960/150‑1 (MQ=255)
                                                                                                                                                  tGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGc                  >  1:299363/1‑135 (MQ=255)
                                                                                                                                                  tGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGc                  <  2:299363/135‑1 (MQ=255)
                                                                                                                                                   gTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATcggcgg  <  2:261878/150‑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 < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: