Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,688,199 A→T 5.3% N325K (AAT→AAA bcsC ← cellulose biosynthesis protein BcsC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,688,1990AT 5.3% 245.4 / 2.9 76N325K (AAT→AAAbcsCcellulose biosynthesis protein BcsC
Reads supporting (aligned to +/- strand):  ref base A (36/36);  new base T (2/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

GGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCGGGAATGGCTTTACCCGCCATACCAGAGTCCACCGCCGCTA  >  NC_000913/3688051‑3688344
                                                                                                                                                    |                                                                                                                                                 
ggAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGAt                                                                                                                                                  >  1:702819/1‑150 (MQ=255)
ggAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGAt                                                                                                                                                  <  1:232784/150‑1 (MQ=255)
   aCAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCg                                                                                                                                               <  1:884436/150‑1 (MQ=255)
         gTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTg                                                                                                                                         >  1:883067/1‑150 (MQ=255)
                  ggTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGgcgc                                                                                                                                <  2:821259/150‑1 (MQ=255)
                     cAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGAt                                                                                                                             >  1:510859/1‑150 (MQ=255)
                       ggATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCACCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCg                                                                                                                           >  2:906424/1‑150 (MQ=255)
                           tattGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCtt                                                                                                                       <  2:636333/150‑1 (MQ=255)
                           tattGGCTTTCAGCGCAGCATCGAGCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCtt                                                                                                                       <  2:564196/150‑1 (MQ=255)
                            attGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCttt                                                                                                                      >  1:675554/1‑150 (MQ=255)
                             ttGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTTCTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTc                                                                                                                     <  2:1039927/150‑1 (MQ=255)
                              tGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCAGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTc                                                                                                                     <  1:924908/149‑1 (MQ=255)
                              tGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCAGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTc                                                                                                                     >  2:924908/1‑149 (MQ=255)
                               ggCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTg                                                                                                                   >  2:552477/1‑150 (MQ=255)
                                gCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTga                                                                                                                  <  2:106552/150‑1 (MQ=255)
                                  tttCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTgaga                                                                                                                >  2:1133142/1‑150 (MQ=255)
                                       gcgcAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACg                                                                                                           <  2:702819/150‑1 (MQ=255)
                                       gcgcAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACg                                                                                                           >  1:17144/1‑150 (MQ=255)
                                       gcgcAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACg                                                                                                           >  1:11846/1‑150 (MQ=255)
                                        cgcAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGc                                                                                                          >  2:846784/1‑150 (MQ=255)
                                         gcagcaTCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGcc                                                                                                         <  2:114938/150‑1 (MQ=255)
                                            gcaTCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCAATGCAGTGgcgaatct                                                                                                                            <  2:826439/128‑6 (MQ=255)
                                            gcaTCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCAATGCAGTGgcgaatct                                                                                                                            >  1:826439/1‑123 (MQ=255)
                                             caTCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGc                                                                                                     <  1:465174/150‑1 (MQ=255)
                                             caTCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGc                                                                                                     <  1:544310/150‑1 (MQ=255)
                                             caTAGCCCTGCTGGATCGCCAGCAAGTACCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTCCCCGTCCAGAGCCAGGGCTGTTTACAGATTCGCCACTGAATTGGCGCGATCGCCTATCTGAGACTACGCCTGGc                                                                                                     >  2:704401/1‑150 (MQ=255)
                                                   ccTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGc                                                                                                     >  1:401193/1‑144 (MQ=255)
                                                   ccTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGc                                                                                                     <  2:401193/144‑1 (MQ=255)
                                                     tgctgGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGcc                                                                                             >  1:21056/1‑150 (MQ=255)
                                                           aTCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGTTTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCaga                                                                                       >  1:1163422/1‑150 (MQ=255)
                                                           aTCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGTTTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCaga                                                                                       >  1:1168359/1‑150 (MQ=255)
                                                           aTCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCaga                                                                                       <  2:684761/150‑1 (MQ=255)
                                                           aTCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCaga                                                                                       <  2:337626/150‑1 (MQ=255)
                                                              gccagccaGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGccc                                                                                            >  1:1302771/1‑142 (MQ=255)
                                                              gccagccaGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGccc                                                                                            <  2:1302771/142‑1 (MQ=255)
                                                              gccagccaGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCt                                                                                    >  1:385792/1‑150 (MQ=255)
                                                               ccagccagTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATTCGCCTGGCCCAGCGCCCCCAGAGCtt                                                                                   <  2:1122837/150‑1 (MQ=255)
                                                                  gccagTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCAc                                                                                <  2:909408/150‑1 (MQ=255)
                                                                     agTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTgg                                                                             >  2:369612/1‑149 (MQ=255)
                                                                     agTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGt                                                                             >  2:949431/1‑150 (MQ=255)
                                                                       tAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCt                                                                           <  2:1181637/150‑1 (MQ=255)
                                                                         gCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCttt                                                                         >  2:986879/1‑150 (MQ=255)
                                                                            gTTTACTTTCAGCAGACTGTTGCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCgg                                                                      <  1:1190119/150‑1 (MQ=255)
                                                                                 cTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCTCCAGAGCTTCACTGTCTTTCGGGTTCg                                                                 >  2:511379/1‑150 (MQ=255)
                                                                                 cTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCg                                                                 <  2:17144/150‑1 (MQ=255)
                                                                                 cTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCg                                                                 <  2:11846/150‑1 (MQ=255)
                                                                                 atttCAGCAGTCTGTTCGATTTGTCGTTGTTGCTGCTGTGAGGGTCTAGTGGGAGGGCTTTTTCCTTATTCGCTACTGGTTTGTCTCGTTGGACTTTATGAGATTATTCATGTCGTAGCGCACCGATAGCTTCACTTTCTTTCGGtttag                                                                 <  2:1355854/149‑3 (MQ=255)
                                                                                  tttCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGTTTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGc                                                                <  2:1168359/150‑1 (MQ=255)
                                                                                  tttCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGTTTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGc                                                                <  2:1163422/150‑1 (MQ=255)
                                                                                   ttCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGcc                                                               <  2:21056/150‑1 (MQ=255)
                                                                                           aCTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGc                                                       <  1:36651/150‑1 (MQ=255)
                                                                                              gTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCtg                                                    >  2:640471/1‑150 (MQ=255)
                                                                                               ttCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCtgt                                                   <  1:1050734/150‑1 (MQ=255)
                                                                                                 ccATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCtgttg                                                 >  2:257591/1‑150 (MQ=255)
                                                                                                   aTTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTa                                               <  2:675554/150‑1 (MQ=255)
                                                                                                   aTTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTa                                               <  2:475860/150‑1 (MQ=255)
                                                                                                         cGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCgg                                         <  1:986879/150‑1 (MQ=255)
                                                                                                         cGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCgg                                         >  2:47327/1‑150 (MQ=255)
                                                                                                          gttgttGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCggg                                        <  2:510859/150‑1 (MQ=255)
                                                                                                            tgttgCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCGGGaa                                      >  1:897148/1‑150 (MQ=255)
                                                                                                              ttgCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCGGGAATg                                    <  1:50444/150‑1 (MQ=255)
                                                                                                               tgctgcTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCGGGAATgg                                   >  1:723316/1‑150 (MQ=255)
                                                                                                               tgctgcTGTGCGGGTCCAGTGCGAGGGATTTTTCCAGATTCGCCACTGCAATGGCGCGATCCCCATTCTGACAATACGCCTGGCCCAGCGCCCCCAGAGCTTCCCTGTCTTTCGGGTTCGCCGGCACAGCCTGATGTAGTTCGGGAATgg                                   >  2:1232098/1‑150 (MQ=255)
                                                                                                                gctgctGTGTGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCg                                                                 >  1:212070/1‑119 (MQ=255)
                                                                                                                gctgctGTGTGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCg                                                                 <  2:212070/119‑1 (MQ=255)
                                                                                                                 ctgctgTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCGGGAATGGCt                                 >  2:134471/1‑150 (MQ=255)
                                                                                                                     tgtgCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCGGGAATGGCTTTAc                             >  2:233417/1‑150 (MQ=255)
                                                                                                                        gCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCGGGAATGGCTTTACCCg                          >  1:813591/1‑150 (MQ=255)
                                                                                                                        gCGGGTCCAGTGAGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCGGGAATGGCTTTACCCg                          >  2:909032/1‑150 (MQ=255)
                                                                                                                              ccAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCGGGAATGGCTTTACCCGcc                        >  1:74641/1‑146 (MQ=255)
                                                                                                                              ccAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCGGGAATGGCTTTACCCGcc                        <  2:74641/146‑1 (MQ=255)
                                                                                                                               cAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCGGGAATGGCTTTACCCGCCATAcc                   <  2:883067/150‑1 (MQ=255)
                                                                                                                                         gCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGACTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCGGGAATGGCTTTACCCGCCATACCAGAGTCCAcc         >  2:980182/1‑150 (MQ=255)
                                                                                                                                          cTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCGGGAATGGCTTTACCCGCCATACCAGAGTCCAccg        <  1:233417/150‑1 (MQ=255)
                                                                                                                                           tttttCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCGGGAATGGCTTTACCCGCCATACCAGAGTCCAccgc       >  1:192545/1‑150 (MQ=255)
                                                                                                                                                ccAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCGGGAATGGCTTTACCCGCCATACCAGAGTCCACCGCCGCTa  >  2:1023887/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGAACAGGCGTTCTGCCCGGTCAGGATTATTGGCTTTCAGCGCAGCATCGCCCTGCTGGATCGCCAGCCAGTAGCGGTTTACTTTCAGCAGACTGTTCCATTTGTCGTTGTTGCTGCTGTGCGGGTCCAGTGCGAGGGCTTTTTCCAGATTCGCCACTGCATTGGCGCGATCGCCTTTCTGAGAATACGCCTGGCCCAGCGCCCCCAGAGCTTCACTGTCTTTCGGGTTCGCCCGCACCGCCTGTTGTAGTTCGGGAATGGCTTTACCCGCCATACCAGAGTCCACCGCCGCTA  >  NC_000913/3688051‑3688344

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: