Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,147 G→C 23.9% intergenic (‑81/+514) lacZ ← / ← lacI pseudogene, truncated/DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273360,1470GC23.9% 204.4 / 18.9 84intergenic (‑81/+514)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base C (17/4);  total (50/34)
Fisher's exact test for biased strand distribution p-value = 2.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATT  >  CP009273/360002‑360293
                                                                                                                                                  |                                                                                                                                                    
tCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGcc                                                                                                                                                  <  2:126983/150‑1 (MQ=255)
gatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGcc                                                                                                                                                  >  1:168224‑M1/103‑150 (MQ=255)
 atttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCg                                                                                                                                                 >  1:864272‑M1/102‑150 (MQ=255)
  tttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGt                                                                                                                                                >  2:111959‑M1/101‑150 (MQ=255)
      ggCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAGTTCCTGGCAGTTTATGGCGGGCATCCTGCCCGCCACCCTCCGGGCCGTTGCt                                                                                                                                            <  2:863516/150‑1 (MQ=255)
              aggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                   >  1:169555‑M1/89‑150 (MQ=255)
              aCCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGCCCTGCCCCCCCCCCCCCCCGCCGCTGCTTCGCaaaa                                                                                                                                   >  2:720895/1‑148 (MQ=255)
               ccGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGt                                                                                                                                  <  2:290759/150‑1 (MQ=255)
               ccGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCCCCCGGCCCTTTCTTCGCAAcct                                                                                                                                  >  2:426135/1‑148 (MQ=38)
                  gagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                                                                                               >  1:689177‑M1/85‑150 (MQ=255)
                        tcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGGGTCCCGCCCCCCCCCCCCCCGGCCGGTGCTTCTCAAAGTTCAAAtccc                                                                                                                         <  2:187899‑M1/72‑2 (MQ=255)
                        aGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                         <  1:656195/150‑1 (MQ=255)
                        aGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCGGGCCGTTGCTTCGCCACGTTCAAATca                                                                                                                          >  1:826266/1‑148 (MQ=255)
                        aGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGCGGGCCCGCCCGCCCCCCCCCCCGGCGCTTCTTCGCAACGCTCAAAACCg                                                                                                                         >  2:670652/1‑150 (MQ=255)
                         caccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                         >  2:512486‑M1/78‑149 (MQ=255)
                          cGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                       <  2:441631/150‑1 (MQ=255)
                           ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                      <  1:602273/150‑1 (MQ=255)
                               caaacatcagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcgg                                                                                                                  >  2:687465‑M1/72‑150 (MQ=255)
                                aaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                                                                 >  2:87613‑M1/71‑150 (MQ=255)
                                 ccATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                                                                 <  2:141533/149‑1 (MQ=255)
                                 aacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                               >  2:650936‑M1/70‑150 (MQ=255)
                                  cATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                                              <  2:859350/150‑1 (MQ=255)
                                   aTCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                             <  1:802050/150‑1 (MQ=255)
                                    aacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                            >  1:324786‑M1/67‑150 (MQ=255)
                                     acagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                           >  2:870323‑M1/66‑150 (MQ=255)
                                      cagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                          >  2:765315‑M1/65‑150 (MQ=255)
                                       aGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                         <  1:400947/150‑1 (MQ=255)
                                            caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                    <  1:428166/150‑1 (MQ=255)
                                            caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTCCCAGCCACCCCCCCGGGCGGTTCTTCGCAAAGTTCAAAACCGCTCCCGGCGGGGTTGTTCtt                                                                                                    >  1:670058/1‑149 (MQ=25)
                                              ccaTCCAGTGCGGGACCTCCGATGATACGCTGTCAGACATGAGAATTCCAATCCCCAATTCCTGGCGGTTTATGACGGGCGGCCTGCCCGCCACCCTCCGGACCGTTGCTTCGCAACGTTCAGATCCGCTCCCAGCGGATTTGTCCTACt                                                                                                  <  2:809301/150‑1 (MQ=18)
                                                 tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                               <  1:981376/150‑1 (MQ=255)
                                                  ccAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCCGTGCCTTGCAAAGTTCAAATCCGCTCCCGGCGGGTGTGTCCTACTctg                                                                                               >  2:864272/1‑147 (MQ=21)
                                                   aggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGga                                                                                             >  1:588050‑M1/52‑150 (MQ=255)
                                                     gTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCACCCCCCCCGCCGTTGCTTTGCCAACTTTCAAACCCCCCCCCGCGGGTTTGTCCTCCTCacgg                                                                                             >  2:177472/1‑145 (MQ=255)
                                                      cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCCCCCGGGCGGTGCTTCGCAAAGTTGAAAATCGCTCCCCGCGGATTTGTTCTAATCAGGggag                                                                                          <  2:956617‑M1/102‑1 (MQ=255)
                                                        cAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                         <  1:909252/149‑1 (MQ=255)
                                                        cAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                                        <  2:7081/150‑1 (MQ=255)
                                                          ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCCGGCCGGTGCTTCGCAACCTTCAAAACCGCCCCCCGCGGGTTTTTCCTTCTCCGGAGAGCggt                                                                                      >  1:292077/1‑148 (MQ=21)
                                                           tctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                     >  2:556551‑M1/44‑150 (MQ=255)
                                                             gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCACCCCCCCGGCCGCCGCTTCGCAAAGTTCAAATCCCCCCCCCGCGGGGTTTTTCCTCTCCGGAGAGCGtttcc                                                                                   >  1:870323/1‑147 (MQ=18)
                                                              cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                  <  1:426135/150‑1 (MQ=255)
                                                                 gactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCTCCCGGCCGTTGCTTCGCAACATTCAAATCCCCTCCCCGCGGGTTTTTCCTACTCAGGAGAgagttccccca                                                                                <  1:925558‑M1/112‑6 (MQ=255)
                                                                 cGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCCCCGGGCCGGTGCTTCGCAACGTTTCAATCCCCCCCCCGCGGGTTTTTTCTACTCAGGAGAGCGTTTAccca                                                                                >  1:628686/1‑147 (MQ=18)
                                                                  gatgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCCGCCCGGCCCCCCCCCCCCCGCGCGTTTCTTCTCCACCTTTACATCCCCCCCCCGCCGGCTTGTTCTTCTCCGGGGAGGGGTCaaccgc                                                                               >  2:324786/1‑144 (MQ=255)
                                                                  actgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGaca                                                                              >  2:934028‑M1/37‑150 (MQ=255)
                                                                    tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCCCCCGGCCGGTGCTTTGCCAAGTTCAAAAACGCCCCCGGCGGGGTTTTCCTACTCCGGAGAGAGGTCACCCacaa                                                                             >  2:471044/1‑149 (MQ=18)
                                                                         aGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAaacaac                                                                        <  2:785151/149‑1 (MQ=255)
                                                                               ttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                 >  1:471044‑M1/24‑150 (MQ=255)
                                                                                tttatttgatgcctggggaattAATTCCTGGCAGTTTATGGTGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                <  2:405342‑M1/128‑1 (MQ=255)
                                                                                 ttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                >  1:67004‑M1/22‑149 (MQ=255)
                                                                                 ttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                >  1:689656‑M1/22‑149 (MQ=255)
                                                                                  aCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                              <  2:187015/150‑1 (MQ=255)
                                                                                   atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                              >  1:187899‑M1/20‑149 (MQ=255)
                                                                                    aTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                            <  1:177472/150‑1 (MQ=255)
                                                                                     tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaaa                                                           <  1:835220/150‑1 (MQ=255)
                                                                                      tgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGGGGCCTTCCCCCCACCCCCCCGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGGTTTGTCCTACTCAGGAGAGCGTTCACCGACAAAAAACAAATAAAACGaa                                                            <  1:778234‑M1/132‑1 (MQ=255)
                                                                                       gatgcctggggaattAATTCTTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                          >  2:452127‑M1/16‑149 (MQ=255)
                                                                                        atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                         >  1:938119‑M1/15‑149 (MQ=255)
                                                                                         aaTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCCCTGCTTCTCCACCCTCCAATCCCCCCCCCGCGGGTTTTTCCCTCTCACGAGAGCGTTCACCCCCCAACaaaagaaaaaaagaaaaacc                                                       >  2:442867/1‑146 (MQ=9)
                                                                                            ctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAg                                                    >  2:319768‑M1/11‑150 (MQ=255)
                                                                                             ccGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAg                                                    <  1:735813/149‑1 (MQ=40)
                                                                                               gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                 >  1:136869‑M1/8‑150 (MQ=255)
                                                                                               gATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                 <  1:877977/150‑1 (MQ=38)
                                                                                                ggaattaaTTCCTGGCAGTTTATGGCGGGCGTCCTTCCCGCCCCCCCCCCTGCCCGTGTTTTTCAACCATTAAATCCGCCCCCCGCGGGGTTGTTCTCCTTATGAGAGCGTTTCCCCCCAAACAACAGAAAAAAAAAAAAGACCAATCtt                                                >  2:593756/7‑150 (MQ=2)
                                                                                                 tCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCttt                                               <  1:459376/150‑1 (MQ=37)
                                                                                                  ccccAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCCCCCCCCGGGCCGGTGCTTCGCAAAGTTCAAATCCGCTCCCCGCGGGTTTGTTCTACTCAGGAGAGGGGTTACCCACAAAAAAAAAAAAAAAAGAAAAGGCCAGTCCTTc                                              >  1:138599/1‑150 (MQ=2)
                                                                                                  aattaaTTCCTGGCAGTTTATGGCGGGCGCCCTTCCCCCCACCCTCCCGGGCCTTTCCTCGCAACGCTCCAATCCCCTCCCCGCGCATTTTTCCTCCTCCCGAGAGCGTTCCCCCACAAACCACTGTTTAAACGAAAAGCCCAGCCTTTc                                              >  1:565461/5‑150 (MQ=2)
                                                                                                   cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                              <  2:428065/149‑1 (MQ=35)
                                                                                                   attAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                             >  1:931196‑M1/4‑150 (MQ=255)
                                                                                                      aaTTCCTGGCAGTTTATGGCGGGTGCCCCGCCCCCCACCCCCCCGGGCGGTGTTTCTCAAAGTTTAAAACCGCTCCCCGCGGATTTATCCTTCTCTCGAGGGAGTGGACAGACAAAAAACAGATAAAAAGATAGGCTCAGTCatttgcgt                                          >  2:69124/1‑145 (MQ=2)
                                                                                                      aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                          <  2:292576/150‑1 (MQ=32)
                                                                                                       aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCGCAAAATTTAAAATCCCTCCCCGCGGGGTTTTTCTTCTCAAGGGGGGGTTCTCCCCCAAACAAAAAATAAAAAAAAAAGGCCACTCtttttact                                          >  2:605836/1‑144 (MQ=2)
                                                                                                        ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCCTTGCTTCGCAACAATTAAATCCCCTCCCCGCGGGGTTTTTCTACTCCAGGGAGGGGTTCCCCACAAAAAAAAAAAAAAAAAAAAAGGCCCGTCTTTCGACTgc                                        >  2:665759/1‑149 (MQ=2)
                                                                                                        ttCCTGGCAGTTTATGGCGGGCGGCCTTCCCCCCCCCCTCCCGGGCGTTGCTTCGCAACGTTCAAATCCCCTCCCCGCGGGGTTTTCCTTCTCAGGAGAGCGGTTACCCACAAACAAAAGATAAAAAGAAAAGCCCAGTCtttttact                                          >  1:650936/1‑143 (MQ=2)
                                                                                                         tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCGGGCCGGTGCTTCGCAACGTTCAAATCCCCTCCCCGCGGGTTTTTCCTTCTCAGGAGAGCGTTCACCCACAAACAACAGATAAAACGAAAAGCCCAGGCTTTCGACTGAg                                       >  1:41480/1‑150 (MQ=11)
                                                                                                         tCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCCGCGGGTTTTTCCTACTCAGGAGAGCGGTCACCCACAAAAAACAGATAAAAAAAAAAGCCCCGTCTTTCTACTGAg                                       >  1:448743/1‑150 (MQ=11)
                                                                                                                                    gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg              >  2:479684/1‑149 (MQ=11)
                                                                                                                                     cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg             >  2:28080/1‑149 (MQ=31)
                                                                                                                                       cgacACCCTGCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAGCGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat          <  2:223351/147‑1 (MQ=16)
                                                                                                                                       cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat          >  2:251238/1‑150 (MQ=33)
                                                                                                                                        gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt         >  1:357461/1‑150 (MQ=34)
                                                                                                                                           aCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattggg      <  2:972896/150‑1 (MQ=37)
                                                                                                                                             ccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaa    <  1:238075/150‑1 (MQ=38)
                                                                                                                                                 cGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaatt  <  1:498112/148‑1 (MQ=40)
                                                                                                                                                  |                                                                                                                                                    
TCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATT  >  CP009273/360002‑360293

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