Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,376,953 G→C 8.9% 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 8.9% 249.1 / 12.0 79C41S (TGT→TCT) ‡ycjSputative oxidoreductase YcjS
Reads supporting (aligned to +/- strand):  ref base G (35/36);  new base C (3/4);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

TGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGA  >  NC_000913/1376805‑1377102
                                                                                                                                                    |                                                                                                                                                     
tGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGctgtctgt                                                                                                                                                      <  1:1360087/150‑1 (MQ=255)
 ggTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCtgtg                                                                                                                                                     >  2:1079871/1‑150 (MQ=255)
 gcttGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGCTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCtgtg                                                                                                                                                     >  2:184506/3‑150 (MQ=255)
     gCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAg                                                                                                                                                 <  1:648825/150‑1 (MQ=255)
      cGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAgc                                                                                                                                                <  1:63763/150‑1 (MQ=255)
         aCCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCACAGACAgccac                                                                                                                                             <  2:543657/150‑3 (MQ=255)
           cTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCACAGACAgccacca                                                                                                                                           <  1:992981/150‑2 (MQ=255)
             gCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCttt                                                                                                                                         >  2:982793/1‑150 (MQ=255)
              cTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTc                                                                                                                                        <  1:606454/150‑1 (MQ=255)
                aaGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTccc                                                                                                                                      >  1:1253036/1‑150 (MQ=255)
                aaGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGCGGCTGTCACAGACAgccaccagttcc                                                                                                                                      <  1:669464/150‑7 (MQ=255)
                 agagGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCca                                                                                                                                     >  1:576332/1‑150 (MQ=255)
                     gTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcg                                                                                                                                 <  2:284949/150‑1 (MQ=255)
                     gTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcg                                                                                                                                 <  1:946249/150‑1 (MQ=255)
                       cTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcgca                                                                                                                               >  1:211119/1‑150 (MQ=255)
                        tttGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcgcag                                                                                                                              >  1:319631/1‑150 (MQ=255)
                           gtgAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcgcaggcg                                                                                                                           >  2:1128472/1‑150 (MQ=255)
                           gtgAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCcaggcgcaggcg                                                                                                                           >  2:326131/1‑150 (MQ=255)
                            tgAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCACCGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCaggcgcaggcgc                                                                                                                          >  1:386109/1‑150 (MQ=255)
                              aaaaGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTg                                                                                                                        <  2:1010876/150‑1 (MQ=255)
                               aaaGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTgg                                                                                                                       >  1:552213/1‑150 (MQ=255)
                               aaaGTGCAATAACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTgg                                                                                                                       >  2:688185/1‑150 (MQ=255)
                                aaGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGc                                                                                                                      <  1:223347/150‑1 (MQ=255)
                                aaGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGc                                                                                                                      <  1:982793/150‑1 (MQ=255)
                                aaGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGc                                                                                                                      <  1:1221286/150‑1 (MQ=255)
                                 aGTGCATTGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGGCAGGGGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCATCGATCTGGAACTGGTGGCTGTCTGTTACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCa                                                                                                                     <  1:754789/150‑1 (MQ=255)
                                     cAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGaaa                                                                                                                 <  1:52358/150‑1 (MQ=255)
                                         gACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATa                                                                                                             <  1:831963/150‑1 (MQ=255)
                                          aCAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATAc                                                                                                            >  1:340683/1‑150 (MQ=255)
                                          aCAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATAc                                                                                                            <  1:892723/150‑1 (MQ=255)
                                             aGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACggg                                                                                                         >  1:1316207/1‑150 (MQ=255)
                                             aGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACggg                                                                                                         <  2:725302/150‑1 (MQ=255)
                                                tctcCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAAt                                                                                                      >  2:1373578/1‑150 (MQ=255)
                                                  tcCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGc                                                                                                    >  1:1302766/1‑150 (MQ=255)
                                                     gCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATc                                                                                                 <  1:743166/150‑1 (MQ=255)
                                                      cTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATcc                                                                                                >  2:1037632/1‑150 (MQ=255)
                                                       tGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCg                                                                                               <  2:1376098/150‑1 (MQ=255)
                                                       tGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGGTCATGCTTCGTACTACTGCACCCGCAACGATCTGTAACTGGTGGCTGTCTTTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCg                                                                                               <  2:1376102/150‑1 (MQ=255)
                                                         cGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATg                                                                                                     <  2:363272/142‑1 (MQ=255)
                                                         cGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATg                                                                                                     >  1:363272/1‑142 (MQ=255)
                                                         cGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATTCATCCgtg                                                                                             <  1:1128472/150‑1 (MQ=255)
                                                         cGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCgtg                                                                                             <  1:326131/150‑1 (MQ=255)
                                                           ggTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGtgtg                                                                                           <  2:1316207/150‑1 (MQ=255)
                                                               gcgATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGAc                                                                                       <  2:1110628/150‑1 (MQ=255)
                                                               gcgATAATAGGCGCAGGCAAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGAACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGAc                                                                                       >  2:53884/1‑150 (MQ=255)
                                                                      taGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGc                                                                                >  2:1204882/1‑150 (MQ=255)
                                                                                  aGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCAtgctgc                                                                    <  2:340683/150‑1 (MQ=255)
                                                                                        cGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCATGTGTTTTTTATTGAATAAAAATACTGGTATTTATACGTGTGGGACGATCCTTAGGCTATTcttctgtggg                                                              <  2:576332/150‑5 (MQ=255)
                                                                                        cGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAAt                                                                                                      <  1:951732/110‑1 (MQ=255)
                                                                                        cGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAAt                                                                                                      >  2:951732/1‑110 (MQ=255)
                                                                                        cGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCgg                                                              >  2:869922/1‑150 (MQ=255)
                                                                                          gATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTg                                                            <  1:184506/150‑1 (MQ=255)
                                                                                          gATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTg                                                            <  1:688185/150‑1 (MQ=255)
                                                                                             aaGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGaaa                                                         <  1:1037632/150‑1 (MQ=255)
                                                                                              agattCATGCTTCGTACTACTGCAACAGCAACGATCTTGAACTCGTGGCTGTCTGTAACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAAc                                                        >  2:780800/4‑150 (MQ=255)
                                                                                                   cATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGAt                                                   >  1:716088/1‑150 (MQ=255)
                                                                                                   cATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGAt                                                   <  1:139013/150‑1 (MQ=255)
                                                                                                   cATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGAt                                                   >  1:1111186/1‑150 (MQ=255)
                                                                                                      gCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGAtgtg                                                <  1:186112/150‑1 (MQ=255)
                                                                                                      gCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGAtgtg                                                >  2:338716/1‑150 (MQ=255)
                                                                                                      gCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGAtgtg                                                <  1:1204882/150‑1 (MQ=255)
                                                                                                       cTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTgg                                               >  1:526492/1‑150 (MQ=255)
                                                                                                       cTTCGAACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTgg                                               >  1:1238569/1‑150 (MQ=255)
                                                                                                           gTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAg                                           >  1:559463/1‑150 (MQ=255)
                                                                                                            tactacTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGc                                          <  2:211119/150‑1 (MQ=255)
                                                                                                                   gCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCt                                   >  1:1237501/1‑150 (MQ=255)
                                                                                                                      cccGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCAc                                >  1:942099/1‑150 (MQ=255)
                                                                                                                      cccGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCAc                                >  1:237762/1‑150 (MQ=255)
                                                                                                                       ccGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCAcc                               <  2:526492/150‑1 (MQ=255)
                                                                                                                        cGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCt                              >  2:765505/1‑150 (MQ=255)
                                                                                                                         gCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTa                             >  1:601262/1‑150 (MQ=255)
                                                                                                                          cAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCACAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGTGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTaa                            >  1:1066620/1‑150 (MQ=255)
                                                                                                                            aCGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAAcc                          <  2:529785/150‑1 (MQ=255)
                                                                                                                                 cTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGtttt                     <  1:53884/150‑1 (MQ=255)
                                                                                                                                  tGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCACAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGTGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTc                    <  2:1066620/150‑1 (MQ=255)
                                                                                                                                   ggAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCa                   >  1:1018612/1‑150 (MQ=255)
                                                                                                                                    gAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCAt                  <  1:937421/150‑1 (MQ=255)
                                                                                                                                        aggcggCTGTCACAGACAGCCACCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCt                              <  2:71741/130‑1 (MQ=255)
                                                                                                                                        aggcggCTGTCACAGACAGCCACCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCt                              >  1:71741/5‑134 (MQ=255)
                                                                                                                                         ggtggCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGa             >  2:1377410/1‑150 (MQ=255)
                                                                                                                                          gcggCTGTCACAGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGTTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGaa            >  1:1297867/3‑150 (MQ=255)
                                                                                                                                          gcggCTGTCACAGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGaa            >  1:431615/3‑150 (MQ=255)
                                                                                                                                                    gtgACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGa  >  1:104161/1‑150 (MQ=255)
                                                                                                                                                    gtgACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGa  >  1:99039/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 < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: