Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,376,953 G→C 5.5% C41S (TGT→TCT) ‡ ycjS → putative oxidoreductase YcjS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,9530GC 5.5% 490.5 / 12.9 146C41S (TGT→TCT) ‡ycjSputative oxidoreductase YcjS
Reads supporting (aligned to +/- strand):  ref base G (76/60);  new base C (3/5);  total (79/67)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 161 total aligned reads displayed.

TGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGA  >  NC_000913/1376805‑1377102
                                                                                                                                                    |                                                                                                                                                     
tGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGctgtctgt                                                                                                                                                      >  2:1404350/1‑150 (MQ=255)
tGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGctgtctgt                                                                                                                                                      >  1:1793788/1‑150 (MQ=255)
tGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGctgtctgt                                                                                                                                                      <  1:1606066/150‑1 (MQ=255)
      cGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAgc                                                                                                                                                >  1:698864/1‑150 (MQ=255)
          ccTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCACAGACAgccacc                                                                                                                                            <  1:659981/150‑1 (MQ=255)
          ccTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGCGGCTGTCACAGACAgccacc                                                                                                                                            <  1:1929447/150‑1 (MQ=255)
              cTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTc                                                                                                                                        <  1:913499/150‑1 (MQ=255)
               tAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTTCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTcc                                                                                                                                       <  2:1403138/150‑1 (MQ=255)
               tAAGAGGTCTTTCTGAAAAGTGCAATGACAAGCTCTCCGCAGCGGGTCTCGATAATAGGCGCAGGCCAGGTGGCGGAAAAGGTACATGCTTCGTACTACTGCACCCGCAACCATCTGGAACTGGTGGCTGTCTGTGACACCCGCCTTTcc                                                                                                                                       >  2:1425545/1‑150 (MQ=255)
                 agagGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCca                                                                                                                                     <  2:864673/150‑1 (MQ=255)
                 agagGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCca                                                                                                                                     <  1:1247247/150‑1 (MQ=255)
                  gagGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCACGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTcaaag                                                                                                                                    >  1:1727933/1‑146 (MQ=255)
                    ggTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggc                                                                                                                                  >  1:748228/1‑150 (MQ=255)
                     gTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcg                                                                                                                                 >  2:975936/1‑150 (MQ=255)
                     gTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcg                                                                                                                                 <  2:2342824/150‑1 (MQ=255)
                     gTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcg                                                                                                                                 >  2:2445234/1‑150 (MQ=255)
                       cTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcgca                                                                                                                               >  2:259058/1‑150 (MQ=255)
                        tttGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcgcag                                                                                                                              >  2:1095283/1‑150 (MQ=255)
                            tgAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCaggcgcaggcgc                                                                                                                          >  2:1384016/1‑150 (MQ=255)
                             gAAAAGTGCATTTACAAGCTCTCCGCTGCGTGTCTCGATAATAGGCGCAGGCGAGGTGGGGGATAAGGTTCATGATGCGTAATACTGCTCACTCAACTATCTGGAACTGTTGGATTTCTTTTACATCCGCGTTTCCCAGGCGCAGGCGCt                                                                                                                         <  2:2315448/150‑1 (MQ=255)
                             gAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCt                                                                                                                         >  2:361622/1‑150 (MQ=255)
                                aaGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGc                                                                                                                      <  1:67594/150‑1 (MQ=255)
                                aaGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGc                                                                                                                      <  2:1839693/150‑1 (MQ=255)
                                aaGTGCAATGACAAGCTCTCCGCTGCGGGGCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGc                                                                                                                      <  2:2141807/150‑1 (MQ=255)
                                  gTGCAATGACAAGCTCTCCTCTTCGGGTAGCGATAATAGGCGCAGGCCAGGTGGCGGATTAGGGTCATGCTTCGTAATACTTCCCCCTCAATGAGCTGGAAGTGGTGGCTCTCTGTTTCACCCGCCTTTCCCAGGCGCAGGCGCTGGCAg                                                                                                                    <  2:1589899/150‑1 (MQ=255)
                                    gCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGCGGCTGTCACAGACAgccacc                                                                                                                                            >  2:1140617/1‑124 (MQ=255)
                                    gCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGCGGCTGTCACAGACAgccacc                                                                                                                                            <  1:1140617/124‑1 (MQ=255)
                                        tGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAAt                                                                                                              >  1:821279/1‑150 (MQ=255)
                                         gACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACGGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATa                                                                                                             >  2:79785/1‑150 (MQ=255)
                                           cAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACg                                                                                                           >  2:1249940/1‑150 (MQ=255)
                                           cAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACg                                                                                                           >  1:681238/1‑150 (MQ=255)
                                             aGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACggg                                                                                                         >  1:386536/1‑150 (MQ=255)
                                             aGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACggg                                                                                                         <  1:361622/150‑1 (MQ=255)
                                               ctctCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGaa                                                                                                       >  2:1461393/1‑150 (MQ=255)
                                                   ccGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAATTGGTGGCTGTGTGTGAGAGCCGACTTTGCGAGTCTCAGGCGCTGTCAGATAAATTCGGGTATgaa                                                                                                   <  2:507966/150‑3 (MQ=255)
                                                      cTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATcc                                                                                                <  1:1000160/150‑1 (MQ=255)
                                                       tGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCg                                                                                               <  2:436417/150‑1 (MQ=255)
                                                           ggTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGtgtg                                                                                           <  1:220289/150‑1 (MQ=255)
                                                            gTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTgg                                                                                          >  2:1144538/1‑150 (MQ=255)
                                                              cgcgGTAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGa                                                                                        >  2:1610195/1‑150 (MQ=255)
                                                              cgcgATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGa                                                                                        >  2:103094/1‑150 (MQ=255)
                                                               gcgATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGAc                                                                                       <  1:2342554/150‑1 (MQ=255)
                                                               gcgATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGAc                                                                                       <  1:1760788/150‑1 (MQ=255)
                                                                cgATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACg                                                                                      >  2:1934988/1‑150 (MQ=255)
                                                                cgATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACg                                                                                      >  1:118877/1‑150 (MQ=255)
                                                                   taataGGCGCAGGCCAGGTGGCGGATAAGGTTGATGCTTCGTACTACTTCACCCGCAAGTATCTGGAACTGGTGGCTGTCTGTGACAGCCGCGTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATc                                                                                   <  1:1509166/150‑1 (MQ=255)
                                                                      taGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGc                                                                                >  2:1467098/1‑150 (MQ=255)
                                                                        ggCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAg                                                                              >  1:2045187/1‑150 (MQ=255)
                                                                        ggCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAg                                                                              >  1:1630872/1‑150 (MQ=255)
                                                                          cgcAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGc                                                                            >  1:1023420/1‑150 (MQ=255)
                                                                          cgcAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGc                                                                            <  1:2448313/150‑1 (MQ=255)
                                                                               gCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCAtgc                                                                       <  2:813275/150‑1 (MQ=255)
                                                                                ccAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCAtgct                                                                      >  2:2023303/1‑150 (MQ=255)
                                                                                 cAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCAtgctg                                                                     >  2:2089909/1‑150 (MQ=255)
                                                                                  aGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCAtgctgc                                                                    >  2:2081501/1‑150 (MQ=255)
                                                                                    gTGGCGGATAAGGTTCATGCTTCGTAGTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGctgctg                                                                  <  1:612584/150‑1 (MQ=255)
                                                                                        cGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCgg                                                              >  1:2228323/1‑150 (MQ=255)
                                                                                          gATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTg                                                            <  2:1805765/150‑1 (MQ=255)
                                                                                          gATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTg                                                            <  1:2117777/150‑1 (MQ=255)
                                                                                           aTAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGa                                                           >  2:1668089/1‑150 (MQ=255)
                                                                                               ggTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAAcc                                                       >  1:1143691/1‑150 (MQ=255)
                                                                                               ggTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAAcc                                                       >  2:152450/1‑150 (MQ=255)
                                                                                                 ttCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTg                                                     <  2:2346632/150‑1 (MQ=255)
                                                                                                  tCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGa                                                    <  1:1285486/150‑1 (MQ=255)
                                                                                                  tCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGa                                                    <  1:461586/150‑1 (MQ=255)
                                                                                                    aTGATTCGTACTACTGCACCCGCAAAGATCTGGAACTTGTGTCTTTCTGTGACAGCATCCTTTCCCATGCTCAGGCGCTGGCATAAAAATAAGGGAATTCATCCGTTTGGGACGATCCGCAGGCCATGCTGCTGGCGTTGATACCTGAtt                                                  <  2:2263035/150‑2 (MQ=255)
                                                                                                        ttCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGt                                              <  2:40331/150‑1 (MQ=255)
                                                                                                         tCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGtt                                             >  1:882281/1‑150 (MQ=255)
                                                                                                           gTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAg                                           >  2:1865732/1‑150 (MQ=255)
                                                                                                           gTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAg                                           <  2:1386744/150‑1 (MQ=255)
                                                                                                              ctactGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGt                                        >  2:779847/1‑150 (MQ=255)
                                                                                                                actGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCCGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCt                                      <  1:204841/150‑1 (MQ=255)
                                                                                                                actGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCt                                      >  2:1557229/1‑150 (MQ=255)
                                                                                                                  tGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGc                                    <  1:861029/150‑1 (MQ=255)
                                                                                                                  tGCACCCGCAACGAGCTGGAACTGGTGGCTGGCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGc                                    <  2:604834/150‑1 (MQ=255)
                                                                                                                   gCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCt                                   >  2:1931005/1‑150 (MQ=255)
                                                                                                                   gCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCt                                   >  2:719277/1‑150 (MQ=255)
                                                                                                                   gCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCt                                   >  2:2469317/1‑150 (MQ=255)
                                                                                                                     aCCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCa                                 >  1:234357/1‑150 (MQ=255)
                                                                                                                        cGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCt                              >  1:561949/1‑150 (MQ=255)
                                                                                                                          cAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCAcc                               >  1:486588/1‑147 (MQ=255)
                                                                                                                           aaCGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGTACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAAc                           <  2:750382/150‑1 (MQ=255)
                                                                                                                            aCGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCAGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAAcc                          >  1:1385524/1‑150 (MQ=255)
                                                                                                                             cGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCg                         <  2:722430/150‑1 (MQ=255)
                                                                                                                               aTCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGtt                       >  1:1376330/1‑150 (MQ=255)
                                                                                                                                  tGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCCGGCCATGCTGCTGGCGGGGAACCCTGATGTGGTTAGCGTCTGCTCCCCTACCCGTTTTc                    >  1:1512233/1‑150 (MQ=255)
                                                                                                                                       cTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTAc               <  1:824029/150‑1 (MQ=255)
                                                                                                                                         ggtggCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGa             <  2:2468935/150‑1 (MQ=255)
                                                                                                                                         ggcggCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGa             >  1:709806/4‑150 (MQ=255)
                                                                                                                                         ggcggCTGTCACAGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCAt                  <  2:1781336/142‑1 (MQ=255)
                                                                                                                                           tggCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAAc           <  2:821279/150‑1 (MQ=255)
                                                                                                                                             gCTGTCTGTGACAGCCGCCTTTCCCAGGAGCAGGCGCTGGCAGAAAAATACGTGAATGCATCCCTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACaa         >  2:922574/1‑149 (MQ=255)
                                                                                                                                               tgtctgtGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATAc       >  2:2062950/1‑150 (MQ=255)
                                                                                                                                                gtctgtGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATAcc      >  1:1237462/1‑150 (MQ=255)
                                                                                                                                                 tctgtGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAGCATAccc     >  1:648293/1‑150 (MQ=255)
                                                                                                                                                 tctgtGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATAccc     <  2:386536/150‑1 (MQ=255)
                                                                                                                                                 tctgtGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATAccc     <  1:1610195/150‑1 (MQ=255)
                                                                                                                                                 tctgtGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATAccc     <  2:893435/150‑1 (MQ=255)
                                                                                                                                                  ctgtGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAGAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCt    >  2:678413/1‑150 (MQ=255)
                                                                                                                                                   tgtgACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCCCCTAACCGTTTTCATTACGAACATACCCTg   >  1:188446/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGA  >  NC_000913/1376805‑1377102

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