Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,376,952 T→A 8.9% C41S (TGT→AGT) ‡ ycjS → putative oxidoreductase YcjS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,9520TA 8.9% 235.7 / 10.9 79C41S (TGT→AGT) ‡ycjSputative oxidoreductase YcjS
Reads supporting (aligned to +/- strand):  ref base T (35/37);  new base A (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.78e-01

CCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAA  >  NC_000913/1376803‑1377092
                                                                                                                                                     |                                                                                                                                            
ccTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGctgtct                                                                                                                                              >  2:352665/1‑150 (MQ=255)
  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:892723/150‑1 (MQ=255)
                                            aCAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATAc                                                                                                  >  1:340683/1‑150 (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                                                                                           >  1:363272/1‑142 (MQ=255)
                                                           cGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATg                                                                                           <  2:363272/142‑1 (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:1111186/1‑150 (MQ=255)
                                                                                                     cATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGAt                                         >  1:716088/1‑150 (MQ=255)
                                                                                                     cATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGAt                                         <  1:139013/150‑1 (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:237762/1‑150 (MQ=255)
                                                                                                                        cccGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCAc                      >  1:942099/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)
                                                                                                                                                     |                                                                                                                                            
CCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAA  >  NC_000913/1376803‑1377092

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: