Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,074,839 T→C 5.9% H57H (CAT→CAC yihW → putative transcriptional regulator for sulphoquinovose utilization

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,074,8390TC 5.9% 298.1 / 7.2 102H57H (CAT→CACyihWputative transcriptional regulator for sulphoquinovose utilization
Reads supporting (aligned to +/- strand):  ref base T (57/38);  new base C (3/3);  total (60/42)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.16e-01

Only 100 of 104 total aligned reads displayed.

TAACCCGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGAC  >  NC_000913/4074692‑4074988
                                                                                                                                                   |                                                                                                                                                     
tAACCCGCGGCACGACCAACTCCTCATGCTGATCGCAGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCa                                                                                                                                                     >  2:1214217/1‑150 (MQ=255)
  aCCCGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACg                                                                                                                                                   >  1:1366060/1‑150 (MQ=255)
   cccGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACgg                                                                                                                                                  >  1:5665/1‑150 (MQ=255)
     cgcgGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCAGCCGGGATATTCGTAAATTAAGCGTGCAAGGCCTGATTACGCGCCATCTCGGTg                                                                                                                                                >  2:1609628/1‑150 (MQ=255)
       cgGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGCCGGTCCGCCGGGATATTCGTAAATTATGCGAGCAAGGCCTGATTACGCGCCATCACGGTGgc                                                                                                                                              >  1:1285636/1‑150 (MQ=255)
         gCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGgcgc                                                                                                                                            >  1:696126/1‑150 (MQ=255)
         gCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGgcgc                                                                                                                                            >  2:483271/1‑150 (MQ=255)
          cACGACCAACTCCTCATGCTTATCGGCTAGGTTGGGTATATTAATATTGATGAGCTGGCAAAGCTGCTGGATGGCTCCACGGATACGGTCCGCCGGGATATTCGTAAATTACGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcg                                                                                                                                           <  1:1092786/150‑1 (MQ=255)
           aCGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgt                                                                                                                                          <  1:1021971/150‑2 (MQ=255)
            cGACCAACTCCTCATGCTGATAGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCAGATTACGCGCCACAACGGTGGCGCggg                                                                                                                                         >  1:822116/1‑150 (MQ=255)
             gACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGt                                                                                                                                        >  2:431807/1‑150 (MQ=255)
             caCCAACTCCTCATGCTCATCCCCGAGCGTGGGTATATGAAAATACATGAGCTGCCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCCTCAATTCAGCGAGCAAGGCCTGATTACGCGCCACCGTGATGGcgcgtaa                                                                                                                                        <  2:334568/149‑4 (MQ=255)
               ccAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCg                                                                                                                                      <  2:595764/150‑1 (MQ=255)
                cAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGCCGGTCCGCCGGGATATTCGTAAATTAAGCGCGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCgg                                                                                                                                     >  2:674174/1‑150 (MQ=255)
                cAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCCTCACGGTGGCGCGGGTCgg                                                                                                                                     >  2:968646/1‑150 (MQ=255)
                cAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCgg                                                                                                                                     <  2:1686091/150‑1 (MQ=255)
                 aaCTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCggg                                                                                                                                    >  2:230644/1‑150 (MQ=255)
                 aaCTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCggg                                                                                                                                    >  1:1105172/1‑150 (MQ=255)
                         aTGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGc                                                                                                                            >  1:51993/1‑150 (MQ=255)
                          tGCTGATCGCCGAGCTTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcg                                                                                                                           >  2:77632/1‑150 (MQ=255)
                          tGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcg                                                                                                                           >  2:807848/1‑150 (MQ=255)
                           gCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcgt                                                                                                                          <  1:1215164/150‑1 (MQ=255)
                           gCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGCCAAATCTGCTGGATGCCTCCACGCAGACGGTCCGCGGGGATATTCGTAATTTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGGGGCGCGGGTCGGGCTTCCAGggt                                                                                                                          <  1:290875/150‑1 (MQ=255)
                             tGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcgtcg                                                                                                                        <  2:207352/150‑1 (MQ=255)
                              taTCGCCGTTCGTGGGTTTATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcgtcgt                                                                                                                       <  2:950919/149‑1 (MQ=255)
                              gATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGcgtcgt                                                                                                                       >  1:843049/1‑150 (MQ=255)
                              gATCGCCGAGCGTGGGTATATGAATAATGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTACGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTCATTACGCGCCATCACGCTGGCGCGCGTCGGGCTTCCAGcctcgt                                                                                                                       >  1:76335/1‑150 (MQ=255)
                               aTCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGtt                                                                                                                      >  1:786424/1‑150 (MQ=255)
                               aTCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGtt                                                                                                                      <  2:1014919/150‑1 (MQ=255)
                                   ccGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGAGGGCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATa                                                                                                                  <  1:726789/150‑1 (MQ=255)
                                    cGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGAAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATAc                                                                                                                 >  1:731205/1‑150 (MQ=255)
                                        cGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCg                                                                                                             <  2:546236/150‑1 (MQ=255)
                                          tGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCTACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGtt                                                                                                           >  2:779559/1‑150 (MQ=255)
                                          tGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGtt                                                                                                           <  2:1539535/150‑1 (MQ=255)
                                              tattTGAATATTGATGCCCTGGCAAATCTGCTGGATTTCTCAACGAAGAGTGTCCGCCGCGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGTGCTTGCAGCGTCGTTAATACGGCGTTCGag                                                                                                       <  2:213452/150‑1 (MQ=255)
                                                       atTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAg                                                                                              <  2:1037085/150‑1 (MQ=255)
                                                        ttGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAgg                                                                                             >  1:1119886/1‑150 (MQ=255)
                                                           atgaGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGttt                                                                                          >  1:566134/1‑150 (MQ=255)
                                                             gaGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCg                                                                                        <  2:786424/150‑1 (MQ=255)
                                                             gaGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTGAATACGGCGTTCGAGCAGCGTGAGGTTTCg                                                                                        >  1:1210083/1‑150 (MQ=255)
                                                               gCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCa                                                                                      <  1:700493/150‑1 (MQ=255)
                                                                 tGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCaaa                                                                                    >  1:593244/1‑150 (MQ=255)
                                                                  ggCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCATAc                                                                                   >  1:1537488/1‑150 (MQ=255)
                                                                  ggCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAAc                                                                                   <  1:1238222/150‑1 (MQ=255)
                                                                   gCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCTAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAAcc                                                                                  <  1:1720105/150‑1 (MQ=255)
                                                                          tgctgGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCCTCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGaaa                                                                           >  2:1484321/1‑150 (MQ=255)
                                                                          tgctgGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGaaa                                                                           >  1:1479367/1‑150 (MQ=255)
                                                                             tgGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGAGTACGCGCCATCACGGTGGCGCGGGGCGGGCTTCCAGCGTCGGTAATACGGCGTTCGAGCAGCGTGACGGTTCGCAAACCGAGGaaaaaa                                                                        >  1:226866/1‑150 (MQ=255)
                                                                               gATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGaaaaaaaa                                                                      >  1:690474/1‑150 (MQ=255)
                                                                                   tctcCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGa                                                                  >  1:512099/1‑150 (MQ=255)
                                                                                   tctcCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGa                                                                  >  1:1683272/1‑150 (MQ=255)
                                                                                    ctcCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGAt                                                                 >  1:480824/1‑150 (MQ=255)
                                                                                    ctcCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGAt                                                                 >  2:1296342/1‑150 (MQ=255)
                                                                                     tcCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGAtt                                                                >  2:922481/1‑150 (MQ=255)
                                                                                      ccACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTg                                                               >  1:623715/1‑150 (MQ=255)
                                                                                      ccACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCGTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTg                                                               >  2:1556641/1‑150 (MQ=255)
                                                                                      ccACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCAGCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTg                                                               <  2:668442/150‑1 (MQ=255)
                                                                                       cccGCAGACGGTCCGCCGGGATATTCGTAAATTACGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGc                                                              <  2:1418713/148‑1 (MQ=255)
                                                                                       cACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGc                                                              >  2:1072712/1‑150 (MQ=255)
                                                                                       cACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGc                                                              >  1:1575118/1‑150 (MQ=255)
                                                                                        aCGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGTTGGCGCGGGTAGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGcc                                                             <  1:936582/150‑1 (MQ=255)
                                                                                        aCGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTCAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGcc                                                             >  2:542071/1‑150 (MQ=255)
                                                                                          gCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGa                                                           >  2:662059/1‑150 (MQ=255)
                                                                                                 gTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGTGTCGGGCTTCCAGCGTCGTAAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAATCCGAGGAAAAAAAAGCGATTGCCGAAGCGGTg                                                    >  2:1379122/1‑150 (MQ=255)
                                                                                                  tCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTgg                                                   <  1:807848/150‑1 (MQ=255)
                                                                                                     gccgGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAg                                                >  1:1336721/1‑150 (MQ=255)
                                                                                                      ccgGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTACTACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAACGCGATTGCCGACGCGGTGGCAGa                                               >  2:1605157/1‑150 (MQ=255)
                                                                                                       ctggATATTGGTAAATTACGCGATCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGAc                                              <  1:1074983/148‑1 (MQ=255)
                                                                                                       cgGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGAc                                              >  2:360343/1‑150 (MQ=255)
                                                                                                       cgGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGAc                                              <  1:1386410/150‑1 (MQ=255)
                                                                                                        gggATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACt                                             >  2:465946/1‑150 (MQ=255)
                                                                                                        gggATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACt                                             <  2:1292707/150‑1 (MQ=255)
                                                                                                               tCGTAACTTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGGGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTcc                                      <  2:851298/150‑1 (MQ=255)
                                                                                                               tCGTAAATTCAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGGGAGGTTTCGCAAACCGAGGAAAAAACAGCGATTGCCGAAGCGGTGGCAGACTATATTcc                                      <  2:733242/150‑1 (MQ=255)
                                                                                                                cGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTACt                                     >  2:791518/1‑150 (MQ=255)
                                                                                                                 gTAAATTATGCGAGCAAGGCCTGATTACGCGCCCTCACGGTGGCGCGGGGCGGGCTTCCACCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTg                                    <  1:601294/150‑1 (MQ=255)
                                                                                                                 gTAAATTAAGCGAGCAAGGCCTGCTTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTg                                    <  2:1017580/150‑1 (MQ=255)
                                                                                                                    aaTTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATg                                 >  1:1123667/1‑150 (MQ=255)
                                                                                                                    aaTTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATg                                 >  2:379852/1‑150 (MQ=255)
                                                                                                                    aaTTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGGGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATg                                 <  2:837395/150‑1 (MQ=255)
                                                                                                                      ttAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGt                               <  1:779559/150‑1 (MQ=255)
                                                                                                                        aatcGATCAAGGCATTATTACGCGCCATCACGGTGGCGCGGGTCTGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTc                             <  2:1617966/147‑1 (MQ=255)
                                                                                                                         aGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTca                            <  2:566134/150‑1 (MQ=255)
                                                                                                                          gCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGAGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTcaa                           >  2:864854/1‑150 (MQ=255)
                                                                                                                           cGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCAGAAGCGGTGGCAGACTATATTCCTGATGGTTcaac                          >  1:829481/1‑150 (MQ=255)
                                                                                                                           cGAGCAAGGCCTGATTACGCGCCATCACGGGGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTcaac                          <  2:1604153/150‑1 (MQ=255)
                                                                                                                                  ggCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATAttt                   >  2:1732284/1‑150 (MQ=255)
                                                                                                                                  ggCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATAttt                   >  2:1442781/1‑150 (MQ=255)
                                                                                                                                  ggCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATAttt                   >  1:902130/1‑150 (MQ=255)
                                                                                                                                    ccTGATTACGCGCCATCACGGTGGCGCGGGTGGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCTAGGAAAAAAAAGCTATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTAt                 <  1:25757/150‑1 (MQ=255)
                                                                                                                                     cTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAAGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATc                >  2:852403/1‑150 (MQ=255)
                                                                                                                                     cTGATTACGCGCCATCACGGGGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATc                <  2:596526/150‑1 (MQ=255)
                                                                                                                                            cgggcCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTAAAAACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGTAAAAATAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTg         <  1:506552/147‑1 (MQ=255)
                                                                                                                                            cgcgCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTg         <  2:1119886/150‑1 (MQ=255)
                                                                                                                                            cgcgCCACCGTGATGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTg         >  1:895468/1‑150 (MQ=255)
                                                                                                                                            cgcgCCACCGTGATGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTg         >  1:1404643/1‑150 (MQ=255)
                                                                                                                                               gcCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACTGGGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAACCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTa      <  1:230644/150‑1 (MQ=255)
                                                                                                                                               gcCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTa      <  1:431807/150‑1 (MQ=255)
                                                                                                                                                cattCACGGTGGGGCGGGTCGGGCGTCGAGGGTCGTTAATACGGCGTTCGAGAAGCGTGATGTTTCGCAAACCGAGGAAAAAATTTCGATTGCAGTATCGGTTGCAGACTATATTCCTGATGGTTCAACAATATTTATGACCATTGGTAc     <  2:1344796/147‑1 (MQ=255)
                                                                                                                                                   tCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGAc  <  2:76335/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TAACCCGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCCAGCGTCGTTAATACGGCGTTCGAGCAGCGTGAGGTTTCGCAAACCGAGGAAAAAAAAGCGATTGCCGAAGCGGTGGCAGACTATATTCCTGATGGTTCAACAATATTTATCACCATTGGTACGAC  >  NC_000913/4074692‑4074988

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