Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 858,665 A→C 6.7% S133A (TCA→GCA)  ybiU ← DUF1479 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913858,6650AC 6.7% 264.1 / 6.7 94S133A (TCA→GCA) ybiUDUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base A (48/38);  new base C (2/5);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

GACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGT  >  NC_000913/858516‑858813
                                                                                                                                                     |                                                                                                                                                    
gACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGa                                                                                                                                                      <  2:1511551/150‑1 (MQ=255)
    cTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAgcg                                                                                                                                                  >  1:723351/1‑150 (MQ=255)
    cTCACATCCGGGTTAAACCATTGCTATCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGCCTGCGCATTCGCCATTTCTTCACTCTGGCGCGCCTGCATTTGCGCCTGCGACCAGTAGATCTGGTAAATCTCGGTACGTGAAgcg                                                                                                                                                  >  1:979567/1‑150 (MQ=255)
        caTCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGCCGATTGACAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGa                                                                                                                                              <  2:1250425/150‑1 (MQ=255)
         aTCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCATGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAg                                                                                                                                             >  1:967915/1‑150 (MQ=255)
          tCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGc                                                                                                                                            >  1:814400/1‑150 (MQ=255)
           ccGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCAATCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCATGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCTCTGAGCg                                                                                                                                           >  2:1274704/1‑150 (MQ=255)
             gtgTTAAACCTTTTCTTTCGATCAGTTTCAACTGTCCAGAGACGATTGAGAAACGTCTGCGCATTCGCCATTTCTTGACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGATGCGCTGAGCGTc                                                                                                                                         <  1:59301/148‑1 (MQ=255)
                 tAAACCATTGCTTTCCATCACTTTCAAATTTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGCTCAGCGCTTCACGTCCCga                                                                                                                                     <  2:1893209/150‑1 (MQ=255)
                 tAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACTACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGCTCAGCGCTTCACGTCCCga                                                                                                                                     <  2:1390416/150‑1 (MQ=255)
                 tAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCga                                                                                                                                     >  1:1669889/1‑150 (MQ=255)
                 tAAACCATTGCTTTCCATCACTTTAAAATGTCCAGAGTCGATTGTGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTATATGTCGGGACGCTCAGCGCTTCACGTCCCga                                                                                                                                     <  1:999106/150‑1 (MQ=255)
                    aCCATTGCTTTTCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGCTCAGCGCTTCACGTCCCgagat                                                                                                                                  <  1:614196/150‑4 (MQ=255)
                      cATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGaaa                                                                                                                                >  1:719917/1‑150 (MQ=255)
                            tttCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCg                                                                                                                          >  1:1507698/1‑150 (MQ=255)
                              tCCATCACATTCAAATGTCCACAGACGATTGAGAACCGACTCCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGcc                                                                                                                        <  2:215286/150‑1 (MQ=255)
                               ccATCACTTTCAAATGTCCACCGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTCAAGCGCTGAGCGTCCCGAAGAAATTATCGCCg                                                                                                                       <  1:368710/150‑1 (MQ=255)
                               ccATCACTTTCAAATGTCCACAGACGATTGAGAACCGACTGCGCATTCGCCATTTATTCACTCTGGCGGGCCTGCATTTGAGCATGCGAACAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCACGAAGAAATTATCGCCg                                                                                                                       >  1:1817500/1‑150 (MQ=255)
                               ccATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGAGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCCAAGAAATTATCGCCg                                                                                                                       >  1:1825730/1‑150 (MQ=255)
                                 aTCACTTTCAATTGCCCACCGCCGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGCCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCggg                                                                                                                     <  2:460102/150‑1 (MQ=255)
                                  tCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCgggg                                                                                                                    >  1:809914/1‑150 (MQ=255)
                                   cACTTTCAAATGTCCACCGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCGGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCGCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGc                                                                                                                   <  1:890496/150‑1 (MQ=255)
                                   cACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGc                                                                                                                   >  2:660168/1‑150 (MQ=255)
                                    cgtttCAATTGTCCACTGACTATTGATAAACGACTCCGCATTCGCCATTTCTTCACTCTTGCGGGCCTGTATTTGCGCCTGCGACCATTAGAGGGGGTAAATCGCGGGACGTGAAGCGCTTTGCGTCCCGAAGAAATTATCGCCGGGGcc                                                                                                                  <  2:870131/148‑1 (MQ=255)
                                     cTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCt                                                                                                                 >  2:1904664/1‑150 (MQ=255)
                                     cTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCt                                                                                                                 >  1:1035524/1‑150 (MQ=255)
                                       ttCAAATGTCCACAGACGATTGAGACACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCttt                                                                                                               >  1:1668600/1‑150 (MQ=255)
                                       ttCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCttt                                                                                                               <  2:751917/150‑1 (MQ=255)
                                           aaTGTCTACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAg                                                                                                           <  2:1071121/150‑1 (MQ=255)
                                           aaTGTCGACAGACGATTGAGAAACGACTGCGCATTCTACATTTCTTGACTCTGTCGGGCCTGAATTTGCGCCTGCGACCAGTAGTTGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAg                                                                                                           <  2:227420/150‑1 (MQ=255)
                                             tGTCCACAGACGATTGAGAAACGACTGCGCATTCGCC‑TTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGAcc                                                                                                        >  1:942959/1‑150 (MQ=255)
                                             tGTCCACAGACGAATGAGAACCGACTGCGCAAACGCCATTTCATCACTCTGGCGCGCCTGCAAATGCGACTGCGACCAGAAGATGAGGTAAATCTCGGGACGTGAAGCGCTGAGCGACCCGAAGAAATTAAGGGCGGGGCCTTTGTAGAc                                                                                                         >  2:343705/1‑150 (MQ=255)
                                                   caGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTc                                                                                                   >  2:1985399/1‑150 (MQ=255)
                                                        gATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGACGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGc                                                                                              <  2:646865/150‑1 (MQ=255)
                                                         aTTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCTGAATATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCg                                                                                             <  1:219404/150‑1 (MQ=255)
                                                          ttGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCgg                                                                                            <  2:23337/150‑1 (MQ=255)
                                                            gagaAACGACTGCGCATTCGCCATTTCTTCCCTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGGGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGtt                                                                                          <  2:588223/150‑1 (MQ=255)
                                                              taaaCGACTGCGCATTCGCCATTTAATCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTTAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGc                                                                                        <  2:863550/149‑1 (MQ=255)
                                                              gaAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGc                                                                                        >  2:65424/1‑150 (MQ=255)
                                                              gaAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGc                                                                                        >  2:1111083/1‑150 (MQ=255)
                                                                aaCGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCgg                                                                                      >  2:1173822/1‑150 (MQ=255)
                                                                 ccGACTGCGAATTCGCCATTTCTTCACTCTGGCGGGCCGGCATTTGCGCCTGCGACCAGTAGATGGTGTAAATCTCGGGACGTTAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGt                                                                                     <  2:1753943/149‑1 (MQ=255)
                                                                    aCTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCa                                                                                  >  2:1402465/1‑150 (MQ=255)
                                                                    aCTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCa                                                                                  <  1:1186554/150‑1 (MQ=255)
                                                                     cTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCGGCGACCAGTAGATGGGGTAAGTCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAg                                                                                 <  2:1358419/150‑1 (MQ=255)
                                                                      tGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGa                                                                                >  1:1634291/1‑150 (MQ=255)
                                                                      tGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGa                                                                                >  1:1375338/1‑150 (MQ=255)
                                                                       gcgcATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGACGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGAt                                                                               <  2:663448/150‑1 (MQ=255)
                                                                            tgccccttttCTTCCCTCTGCCGGGCCTGCATTTGCGCCTGCGACCAGTTGAT‑GGTAAACTCGCGGGACGTGACGCGCTGAGCGTCCGGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCa                                                                         <  1:1544138/143‑1 (MQ=255)
                                                                            agcGCCATTTCTTCACTCTGGCGGGCCTGCATTGGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGCCGGGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTcc                                                                          <  1:620916/148‑1 (MQ=255)
                                                                              cGCCATTTCTTCACGCGGGCGGGCCTGCATTGGCGCCTGCGACCAGTAGATGGGGTAACTCTCGGGACGGGAAGCGCTGAGCGGCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAg                                                                        <  1:1538236/150‑1 (MQ=255)
                                                                               gCCATTTCTTCACTGTGGCGGCCCTGCATATGCGCCTGCGACCAGTAGATGGGGTAAAGCTCGGGCCGGGCAGCGCGGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGc                                                                       <  1:1232440/150‑1 (MQ=255)
                                                                               gCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGc                                                                       >  2:662262/1‑150 (MQ=255)
                                                                               gCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGc                                                                       >  2:272538/1‑150 (MQ=255)
                                                                                   tttCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCg                                                                   >  1:1561911/1‑150 (MQ=255)
                                                                                     tcttcACTCTGGCGGGCCTGCTTTTGCGCCTGCTACCAGTAGATGGGGTAAATCTCGGGTCGTGACGCGCTGAGCGTCCCGAAGAAAATATCGCCGGGGCCTGTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGAc                                                                 <  2:1031186/150‑1 (MQ=255)
                                                                                      cttcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGTCGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCCGCATCGACt                                                                >  2:209950/1‑150 (MQ=255)
                                                                                      cttcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACt                                                                >  2:535261/1‑150 (MQ=255)
                                                                                      cttcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACt                                                                >  2:549133/1‑150 (MQ=255)
                                                                                      attcACGCTGGCGGGCCTGCATTTGCGCCTGCGACATGTAGATGTGTTAAATCTCTGGACGGGATGCGGTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACt                                                                <  1:634098/149‑1 (MQ=255)
                                                                                         cACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGAt                                                             >  2:1881422/1‑150 (MQ=255)
                                                                                          ccgctGGCGGGCCTGCAGTAGCGCCTGCGACCAGTAGATGGGGTAATTCTCGGGACGGGAAGCGCTTAGCGCCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATc                                                            <  1:933380/147‑1 (MQ=255)
                                                                                          aCTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATc                                                            >  2:1412033/1‑150 (MQ=255)
                                                                                            gctGGCGGTCCGGCATTGGCGCCTGGCACCTGTTGTTGGGGTAAATCTCGGGACGTGTACCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCGCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATccc                                                          <  1:1521915/149‑1 (MQ=255)
                                                                                               ggCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGc                                                       >  2:1127709/1‑150 (MQ=255)
                                                                                               ggCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGc                                                       >  2:1399541/1‑150 (MQ=255)
                                                                                                    gCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAgc                                                  <  1:576513/150‑1 (MQ=255)
                                                                                                     ccTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAgcg                                                 <  2:1507698/150‑1 (MQ=255)
                                                                                                         cATTTGCGCCGGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTg                                             <  2:356098/150‑1 (MQ=255)
                                                                                                                 ccTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGACGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCggg                                     <  1:717177/150‑1 (MQ=255)
                                                                                                                 ccTGCGACCAGTAGATGGGGTAAATCTCGGGACGGTAAGCGCGGAGCGGCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCggg                                     <  1:947989/150‑1 (MQ=255)
                                                                                                                      gACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGaaa                                >  2:1595905/1‑150 (MQ=255)
                                                                                                                      gACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGaaa                                >  1:1407997/1‑150 (MQ=255)
                                                                                                                         cAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATgg                             >  2:1644631/1‑150 (MQ=255)
                                                                                                                           gTAGAGGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCTAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCCGCCTAGCGCTTGTTCGCGGGGGAAATGGcg                           >  2:113221/1‑149 (MQ=255)
                                                                                                                               acgcggTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGTAATTATCGCCGGGGCCTTTGTAGAGCTCGTCACAGCGTGTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCtttt                       <  1:99352/146‑1 (MQ=255)
                                                                                                                                gggggTAAATCTCGGGCCGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTa                      <  1:1015561/149‑1 (MQ=255)
                                                                                                                                  gggTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCATTTTAta                    >  2:1719098/1‑149 (MQ=255)
                                                                                                                                    ttAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCAc                  <  2:814400/149‑1 (MQ=255)
                                                                                                                                    gTAAATCTCGGGACTTGAAGCGCTGAGCGTCCAGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATTGTCCAGCATCGACTGATCCCAGCCTAGCGGTTGTTCGCGGGGGAAATGGCCTTTTATCAc                  <  2:1386889/150‑1 (MQ=255)
                                                                                                                                     tAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGGTGCGGGCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCAcc                 <  1:1346609/150‑1 (MQ=255)
                                                                                                                                      aaaTCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCg                >  2:1447994/1‑150 (MQ=255)
                                                                                                                                        tcttcGGGACGCTCAGCGCTTCACGTCCAGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGATCCGCTCCAGATACTCCAGCATAGACTGATCCCAGCCTTGCGCTTGTTCGCGGGCGAAATGGCCTATTTTCCCCGca              >  1:279600/4‑150 (MQ=255)
                                                                                                                                        aTCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGca              >  1:809887/1‑150 (MQ=255)
                                                                                                                                        aTCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGGCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATCGCGTGTTATCACCGca              >  2:1885611/1‑150 (MQ=255)
                                                                                                                                         tctcGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGcac             >  2:1917785/1‑150 (MQ=255)
                                                                                                                                          ctcGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGcaca            >  2:1965888/1‑150 (MQ=255)
                                                                                                                                          cgcGGGACGTGAAGCGCTGAGCGGCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGcaca            <  1:1982451/148‑1 (MQ=255)
                                                                                                                                             gggACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAAcc         >  2:663167/1‑150 (MQ=255)
                                                                                                                                              ggACGTGCAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCGTGTTCGCGGGGGAAATGGCCTTTTATCACCGCGCACCcg        >  2:1222849/1‑150 (MQ=255)
                                                                                                                                                aCGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACcgcg      <  2:809914/150‑1 (MQ=255)
                                                                                                                                                 cGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGa     <  1:1402465/150‑1 (MQ=255)
                                                                                                                                                   tGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACg   >  1:243587/1‑150 (MQ=255)
                                                                                                                                                    gAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGt  >  1:1167894/1‑150 (MQ=255)
                                                                                                                                                    gAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGt  >  2:1868732/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGT  >  NC_000913/858516‑858813

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: ‑