Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,146,763 T→C 5.5% intergenic (‑26/+93) yghW ← / ← yghX DUF2623 domain‑containing protein YghW/putative hydrolase fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,146,7630TC 5.5% 277.5 / 4.8 90intergenic (‑26/+93)yghW/yghXDUF2623 domain‑containing protein YghW/putative hydrolase fragment
Reads supporting (aligned to +/- strand):  ref base T (41/44);  new base C (2/3);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

GTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCA  >  NC_000913/3146615‑3146901
                                                                                                                                                    |                                                                                                                                          
gtgtGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTATCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTc                                                                                                                                           <  2:380154/150‑1 (MQ=255)
 tgtgAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCa                                                                                                                                          <  1:394264/150‑1 (MQ=255)
 tgtgAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCa                                                                                                                                          <  1:775418/150‑1 (MQ=255)
  gtgAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCat                                                                                                                                         <  1:1313129/150‑1 (MQ=255)
        aGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAg                                                                                                                                   >  1:609037/1‑150 (MQ=255)
         gCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGc                                                                                                                                  <  2:1339994/150‑1 (MQ=255)
         gCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGc                                                                                                                                  >  2:611337/1‑150 (MQ=255)
           cTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGa                                                                                                                                >  1:1011090/1‑150 (MQ=255)
            tAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATGACAATATATAGCGaa                                                                                                                               <  1:475361/150‑2 (MQ=255)
               tACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATGACAATATATAGCGaatga                                                                                                                            <  1:335517/150‑5 (MQ=255)
                aCAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATAttt                                                                                                                           >  2:1255709/1‑150 (MQ=255)
                aCAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATGACAATATATAGCGaatgag                                                                                                                           <  2:975720/150‑6 (MQ=255)
                 cAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATAtttt                                                                                                                          >  2:706847/1‑150 (MQ=255)
                 cAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATAtttt                                                                                                                          >  1:958758/1‑150 (MQ=255)
                  aaaaCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTa                                                                                                                         <  2:521881/150‑1 (MQ=255)
                   aaaCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAg                                                                                                                        <  1:1390275/150‑1 (MQ=255)
                   aaaCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAg                                                                                                                        <  1:188206/150‑1 (MQ=255)
                     aCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGc                                                                                                                      <  1:656373/150‑1 (MQ=255)
                     aCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGc                                                                                                                      >  2:1064750/1‑150 (MQ=255)
                     aCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGc                                                                                                                      >  1:1306977/1‑150 (MQ=255)
                      ccGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCa                                                                                                                     <  2:1011090/150‑1 (MQ=255)
                        gcgTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCaaa                                                                                                                   >  1:975743/1‑150 (MQ=255)
                            ttATAATCGTCAGAGAATTTTGTCAGATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAAc                                                                                                               <  2:1314064/150‑1 (MQ=255)
                               taATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCaa                                                                                                                    <  2:558582/142‑1 (MQ=255)
                               taATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCaa                                                                                                                    >  1:558582/1‑142 (MQ=255)
                               taATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGAc                                                                                                            >  2:449888/1‑150 (MQ=255)
                                 aTCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACtt                                                                                                          <  1:565885/150‑1 (MQ=255)
                                      cacaGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTgcg                                                                                                     <  2:808806/150‑1 (MQ=255)
                                      cacaGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTgcg                                                                                                     <  1:1321723/150‑1 (MQ=255)
                                       acaGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTgcgc                                                                                                    <  2:430302/150‑1 (MQ=255)
                                       acaGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTgcgc                                                                                                    >  2:270266/1‑150 (MQ=255)
                                        caGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCa                                                                                                   >  1:1420707/1‑150 (MQ=255)
                                          gAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACg                                                                                                 >  1:805921/1‑150 (MQ=255)
                                            aTTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCt                                                                                               >  1:1006994/1‑150 (MQ=255)
                                            aTTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCt                                                                                               >  1:491262/1‑150 (MQ=255)
                                            aTTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCt                                                                                               >  1:632837/1‑150 (MQ=255)
                                                   cacaTTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTAtttttt                                                                                        >  2:389094/1‑150 (MQ=255)
                                                     caTTAACCGCACTGTCGGCATGCGTTGCTATAAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATAGTTCACATCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCa                                                                                      >  2:552437/1‑150 (MQ=255)
                                                         aaCCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCt                                                                                  >  1:108670/1‑150 (MQ=255)
                                                          aCCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCtt                                                                                 >  2:973531/1‑150 (MQ=255)
                                                              cACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGttt                                                                             <  2:680858/150‑1 (MQ=255)
                                                                 tGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAg                                                                          <  1:293185/150‑1 (MQ=255)
                                                                           cGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCat                                                                                                                                         <  1:95452/77‑1 (MQ=255)
                                                                           cGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCat                                                                                                                                         >  2:95452/1‑77 (MQ=255)
                                                                                       tCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATaaa                                                    <  2:632837/150‑1 (MQ=255)
                                                                                        cccGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAAc                                                   <  2:492147/150‑1 (MQ=255)
                                                                                        cccGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAAc                                                   <  1:1064750/150‑1 (MQ=255)
                                                                                        cccGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAAc                                                   <  2:1306977/150‑1 (MQ=255)
                                                                                         ccGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAAcc                                                  <  2:805921/150‑1 (MQ=255)
                                                                                         ccGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAAcc                                                  >  1:452556/1‑150 (MQ=255)
                                                                                              aTTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCa                                             >  1:1220453/1‑150 (MQ=255)
                                                                                                tAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAgg                                           <  2:140431/150‑1 (MQ=255)
                                                                                                  aGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGgag                                         <  1:449888/150‑1 (MQ=255)
                                                                                                  aGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGgag                                         <  2:1220453/150‑1 (MQ=255)
                                                                                                  aGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGgag                                         <  1:611337/150‑1 (MQ=255)
                                                                                                     cTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGgagaga                                      >  2:534338/1‑150 (MQ=255)
                                                                                                     cTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGgagaga                                      >  2:620759/1‑150 (MQ=255)
                                                                                                       ttACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGata                                    <  2:975743/150‑1 (MQ=255)
                                                                                                       ttACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGata                                    >  2:1073825/1‑150 (MQ=255)
                                                                                                         aCCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACATATAAGGAGATATAtt                                  >  1:1382898/1‑150 (MQ=255)
                                                                                                          ccAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATAttt                                 >  2:1416272/1‑150 (MQ=255)
                                                                                                           cAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTa                                <  2:179936/150‑1 (MQ=255)
                                                                                                           cAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTa                                >  2:272984/1‑150 (MQ=255)
                                                                                                            aaaaTGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTAt                               >  1:521875/1‑150 (MQ=255)
                                                                                                            aaaaTGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTAt                               <  2:1006994/150‑1 (MQ=255)
                                                                                                             aaaTGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATc                              <  2:108670/150‑1 (MQ=255)
                                                                                                              aaTGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCg                             <  1:389094/150‑1 (MQ=255)
                                                                                                              aaTGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCg                             <  2:1202347/150‑1 (MQ=255)
                                                                                                                tGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGaa                           >  1:514898/1‑150 (MQ=255)
                                                                                                                tGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGaa                           <  2:522421/150‑1 (MQ=255)
                                                                                                                 ggTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAAc                          <  2:452556/150‑1 (MQ=255)
                                                                                                                  gTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAAcc                         >  1:968102/1‑150 (MQ=255)
                                                                                                                     aTTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCAtt                      <  1:1255709/150‑1 (MQ=255)
                                                                                                                       tCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAAc                          <  1:390663/144‑1 (MQ=255)
                                                                                                                       tCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAAc                          >  2:390663/1‑144 (MQ=255)
                                                                                                                        cATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCa                   >  1:1236769/1‑150 (MQ=255)
                                                                                                                        cATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCa                   >  1:803805/1‑150 (MQ=255)
                                                                                                                         aTCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAg                  <  2:609037/150‑1 (MQ=255)
                                                                                                                           ccttCACCTCATTCGCAATATATTGTCATATATAGCGATATTTAAGGCAAAGAACCACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAgtg                <  2:233571/148‑1 (MQ=255)
                                                                                                                           cGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAgtg                <  2:958758/150‑1 (MQ=255)
                                                                                                                             ttCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTc              >  1:864504/1‑150 (MQ=255)
                                                                                                                                  cTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTg         >  2:1197188/1‑150 (MQ=255)
                                                                                                                                     aTTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCa      >  1:1267754/1‑150 (MQ=255)
                                                                                                                                     aTTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCa      <  2:1267754/150‑1 (MQ=255)
                                                                                                                                      ttCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAg     >  2:1201044/1‑150 (MQ=255)
                                                                                                                                      atCGCTATATATGACAATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAg     >  2:559977/2‑150 (MQ=255)
                                                                                                                                         gCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCa  <  2:491262/150‑1 (MQ=255)
                                                                                                                                         gCTATATATGACAATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTGTAAGACGTAATCCCGTAATAGATAAACCTATCAGGAGAGATATTTATCGATCCATTTCAGTGTCCTTTGCAATGCa  >  2:475498/1‑150 (MQ=255)
                                                                                                                                               tatTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTc                    <  2:737859/126‑1 (MQ=255)
                                                                                                                                               tatTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTc                    >  1:737859/1‑126 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
GTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCA  >  NC_000913/3146615‑3146901

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