Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,181 G→C 13.0% intergenic (‑115/+480) 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,1810GC13.0% 350.4 / 8.4 120intergenic (‑115/+480)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (49/53);  new base C (15/2);  total (64/56)
Fisher's exact test for biased strand distribution p-value = 2.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 126 total aligned reads displayed.

TACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTAT  >  CP009273/360033‑360330
                                                                                                                                                     |                                                                                                                                                      
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)
    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)
       caaacaactgagaaatttttaattcgctctgaaactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGAtt                                                                                                                                                >  1:893514‑M1/142‑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)
                   ccAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCCGTGCCTTGCAAAGTTCAAATCCGCTCCCGGCGGGTGTGTCCTACTctg                                                                                                                                    >  2:864272/1‑147 (MQ=21)
                    aggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGga                                                                                                                                  >  1:588050‑M1/52‑150 (MQ=255)
                       cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCCCCCGGGCGGTGCTTCGCAAAGTTGAAAATCGCTCCCCGCGGATTTGTTCTAATCAGGggag                                                                                                                               <  2:956617‑M1/102‑1 (MQ=255)
                         cAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                                              <  1:909252/149‑1 (MQ=255)
                         cAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                                                                             <  2:7081/150‑1 (MQ=255)
                          taattcgctctgaaactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                                              >  2:468445‑M1/123‑149 (MQ=255)
                          taattcgctctgaaactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                                                                             <  2:671582‑M1/28‑1 (MQ=255)
                           aattcgctctgaaactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGATAGCg                                                                                                                             >  2:131404‑M1/122‑149 (MQ=255)
                            tctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                                          >  2:556551‑M1/44‑150 (MQ=255)
                            attcgctctgaaactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                                            >  1:279664‑M1/121‑149 (MQ=255)
                              gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCACCCCCCCGGCCGCCGCTTCGCAAAGTTCAAATCCCCCCCCCGCGGGGTTTTTCCTCTCCGGAGAGCGtttcc                                                                                                                        >  1:870323/1‑147 (MQ=18)
                               cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                                                       <  1:426135/150‑1 (MQ=255)
                                gctctgaaactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                                                       >  1:1566‑M1/117‑150 (MQ=255)
                                  gactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCTCCCGGCCGTTGCTTCGCAACATTCAAATCCCCTCCCCGCGGGTTTTTCCTACTCAGGAGAgagttccccca                                                                                                                     <  1:925558‑M1/112‑6 (MQ=255)
                                  cGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCCCCGGGCCGGTGCTTCGCAACGTTTCAATCCCCCCCCCGCGGGTTTTTTCTACTCAGGAGAGCGTTTAccca                                                                                                                     >  1:628686/1‑147 (MQ=18)
                                   gatgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCCGCCCGGCCCCCCCCCCCCCGCGCGTTTCTTCTCCACCTTTACATCCCCCCCCCGCCGGCTTGTTCTTCTCCGGGGAGGGGTCaaccgc                                                                                                                    >  2:324786/1‑144 (MQ=255)
                                     tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCCCCCGGCCGGTGCTTTGCCAAGTTCAAAAACGCCCCCGGCGGGGTTTTCCTACTCCGGAGAGAGGTCACCCacaa                                                                                                                  >  2:471044/1‑149 (MQ=18)
                                       aactgatggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                                                 <  2:714510‑M1/40‑1 (MQ=255)
                                          aGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAaacaac                                                                                                             <  2:785151/149‑1 (MQ=255)
                                             tggcgtagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGa                                                                                                          >  1:243157‑M1/104‑150 (MQ=255)
                                               gcgtagaagcgaagaaataattccctttataccatgtcccgtattaccgcctttgagtgagctgatgccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTGCTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                                                         >  2:582140‑M1/102‑149 (MQ=255)
                                                ttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                                      >  1:471044‑M1/24‑150 (MQ=255)
                                                 tttatttgatgcctggggaattAATTCCTGGCAGTTTATGGTGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                                                     <  2:405342‑M1/128‑1 (MQ=255)
                                                  ttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                                                     >  1:689656‑M1/22‑149 (MQ=255)
                                                  ttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                                                     >  1:67004‑M1/22‑149 (MQ=255)
                                                  tagaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                                                     <  1:268764‑M1/52‑1 (MQ=255)
                                                   aCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                                                   <  2:187015/150‑1 (MQ=255)
                                                    gaagcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                                                    >  1:129821‑M1/97‑149 (MQ=255)
                                                    atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                                                   >  1:187899‑M1/20‑149 (MQ=255)
                                                     aTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                                                 <  1:177472/150‑1 (MQ=255)
                                                      tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaaa                                                                                                <  1:835220/150‑1 (MQ=255)
                                                       gcgaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaaa                                                                                                >  2:571701‑M1/94‑150 (MQ=255)
                                                        gatgcctggggaattAATTCTTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                                                               >  2:452127‑M1/16‑149 (MQ=255)
                                                         gaagaaataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                                                              <  1:961824‑M1/59‑1 (MQ=255)
                                                         atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                                                              >  1:938119‑M1/15‑149 (MQ=255)
                                                              ccGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAg                                                                                         <  1:735813/149‑1 (MQ=40)
                                                               ataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                                                        >  2:73176‑M1/86‑150 (MQ=255)
                                                               ataataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                                                        <  1:748658‑M1/65‑1 (MQ=255)
                                                                gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                                      >  1:136869‑M1/8‑150 (MQ=255)
                                                                 ggaattaaTTCCTGGCAGTTTATGGCGGGCGTCCTTCCCGCCCCCCCCCCTGCCCGTGTTTTTCAACCATTAAATCCGCCCCCCGCGGGGTTGTTCTCCTTATGAGAGCGTTTCCCCCCAAACAACAGAAAAAAAAAAAAGACCAATCtt                                                                                     >  2:593756/7‑150 (MQ=2)
                                                                 aataccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                                      <  2:361489‑M1/67‑1 (MQ=255)
                                                                   ccccAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCCCCCCCCGGGCCGGTGCTTCGCAAAGTTCAAATCCGCTCCCCGCGGGTTTGTTCTACTCAGGAGAGGGGTTACCCACAAAAAAAAAAAAAAAAGAAAAGGCCAGTCCTTc                                                                                   >  1:138599/1‑150 (MQ=2)
                                                                    cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                                   <  2:428065/149‑1 (MQ=35)
                                                                     ccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                                   <  1:724585‑M1/70‑1 (MQ=255)
                                                                     ccctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                                   <  1:454591‑M1/70‑1 (MQ=255)
                                                                       ctttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                                                >  2:68045‑M1/78‑150 (MQ=255)
                                                                       aaTTCCTGGCAGTTTATGGCGGGTGCCCCGCCCCCCACCCCCCCGGGCGGTGTTTCTCAAAGTTTAAAACCGCTCCCCGCGGATTTATCCTTCTCTCGAGGGAGTGGACAGACAAAAAACAGATAAAAAGATAGGCTCAGTCatttgcgt                                                                               >  2:69124/1‑145 (MQ=2)
                                                                       aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                               <  2:292576/150‑1 (MQ=32)
                                                                        tttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                                                 <  1:958683‑M1/72‑1 (MQ=255)
                                                                        tttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                               <  2:287960‑M1/74‑1 (MQ=255)
                                                                         ttataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                                                >  2:958683‑M1/76‑148 (MQ=255)
                                                                         ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCCTTGCTTCGCAACAATTAAATCCCCTCCCCGCGGGGTTTTTCTACTCCAGGGAGGGGTTCCCCACAAAAAAAAAAAAAAAAAAAAAGGCCCGTCTTTCGACTgc                                                                             >  2:665759/1‑149 (MQ=2)
                                                                         ttCCTGGCAGTTTATGGCGGGCGGCCTTCCCCCCCCCCTCCCGGGCGTTGCTTCGCAACGTTCAAATCCCCTCCCCGCGGGGTTTTCCTTCTCAGGAGAGCGGTTACCCACAAACAAAAGATAAAAAGAAAAGCCCAGTCtttttact                                                                               >  1:650936/1‑143 (MQ=2)
                                                                          tataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                               >  2:517214‑M1/75‑148 (MQ=255)
                                                                          tataccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                               >  2:429230‑M1/75‑148 (MQ=255)
                                                                          tCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCCGCGGGTTTTTCCTACTCAGGAGAGCGGTCACCCACAAAAAACAGATAAAAAAAAAAGCCCCGTCTTTCTACTGAg                                                                            >  1:448743/1‑150 (MQ=11)
                                                                            taccatgtcccgtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAg                                                                            >  1:706328‑M1/73‑149 (MQ=255)
                                                                                       gtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTa                                                                <  2:177985‑M1/89‑1 (MQ=255)
                                                                                       gtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTa                                                                <  1:26062‑M1/89‑1 (MQ=255)
                                                                                       gtattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTa                                                                <  2:318476‑M1/89‑1 (MQ=255)
                                                                                        tattaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTAt                                                               <  1:914750‑M1/90‑1 (MQ=255)
                                                                                         attaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTAtt                                                              <  2:395965‑M1/91‑1 (MQ=255)
                                                                                          ttaccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTAttt                                                             >  1:807083‑M1/59‑150 (MQ=255)
                                                                                           taccgcctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTAttt                                                             <  2:671326‑M1/92‑1 (MQ=255)
                                                                                                 ctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCt                                                      <  2:587876‑M1/99‑1 (MQ=255)
                                                                                                 ctttgagtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCt                                                      >  2:663889‑M1/52‑150 (MQ=255)
                                                                                                      cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                                                  >  2:28080/1‑149 (MQ=31)
                                                                                                       gtgagctgataccgctcgccgcagccgaacgccctaggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTgatgcctctagat                                                >  2:42083‑M1/46‑144 (MQ=255)
                                                                                                        cgacACCCTGCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAGCGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                               <  2:223351/147‑1 (MQ=16)
                                                                                                         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)
                                                                                                                      ccGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCt                                 >  2:430079/1‑150 (MQ=255)
                                                                                                                        gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                                                  >  1:5382/1‑131 (MQ=31)
                                                                                                                         ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                                                 <  2:5382/131‑1 (MQ=32)
                                                                                                                          tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTg                               >  2:448819/1‑148 (MQ=255)
                                                                                                                            cTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGa                           <  2:925558/150‑1 (MQ=255)
                                                                                                                               cGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGTGGAATTGGGAATTAATTCTTGAAGAc                          <  1:28080/148‑1 (MQ=255)
                                                                                                                                gCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGaa                       <  1:750191/150‑1 (MQ=255)
                                                                                                                                 cAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGaaa                      <  1:956617/150‑1 (MQ=255)
                                                                                                                                 cAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGaaa                      <  1:405342/150‑1 (MQ=255)
                                                                                                                                      ttCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGcc                 <  2:820523/150‑1 (MQ=255)
                                                                                                                                          aaTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGt             >  2:836828/1‑150 (MQ=255)
                                                                                                                                           ggtctagggCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                                   <  2:168864‑M1/70‑1 (MQ=255)
                                                                                                                                              cGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATa         >  2:212414/1‑150 (MQ=255)
                                                                                                                                                     ggcggATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTAt  <  1:492577/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTAT  >  CP009273/360033‑360330

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