Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,955,039 C→A 5.5% H335N (CAC→AAC)  rep → DNA helicase and single‑stranded DNA‑dependent ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,955,0390CA 5.5% 1281.9 / 32.7 324H335N (CAC→AAC) repDNA helicase and single‑stranded DNA‑dependent ATPase
Reads supporting (aligned to +/- strand):  ref base C (162/144);  new base A (14/4);  total (176/148)
Fisher's exact test for biased strand distribution p-value = 5.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 330 total aligned reads displayed.

AAGCGGCGAACATCCTGATCGCCAATAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAG  >  CP009273/3954890‑3955188
                                                                                                                                                     |                                                                                                                                                     
aaGCGGCGAACATCCTGATCGCCAATAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATc                                                                                                                                                       <  2:2766263/150‑1 (MQ=255)
   cggcgACCATCCTGATCGCCAATAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACt                                                                                                                                                    <  2:3107626/150‑1 (MQ=255)
    ggcgAACATCCTGATCGCCAATAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACtt                                                                                                                                                   <  2:4556427/150‑1 (MQ=255)
      cgAACATCCTGATCGCCAATAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTCTTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCg                                                                                                                                                 >  2:3035821/1‑150 (MQ=255)
       gAACATCCTGATCGCCAATAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGt                                                                                                                                                >  1:4846737/1‑150 (MQ=255)
        aaCATCCTGATCGCCAATAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTc                                                                                                                                               <  2:4922541/150‑1 (MQ=255)
           aTCCTGATCGCCAATAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAAt                                                                                                                                            <  2:2900615/150‑1 (MQ=255)
           aTCCTGATCGCCAATAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATAACTTCGTCaac                                                                                                                                            >  1:2111538/1‑149 (MQ=255)
            tCCTGATCGCCAATAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATa                                                                                                                                           <  2:3503908/150‑1 (MQ=255)
              cTGATCGCCAATAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATaac                                                                                                                                         <  1:658/150‑2 (MQ=255)
               tGATCGCCAATAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATaaaa                                                                                                                                        >  2:1775267/1‑150 (MQ=255)
               tGATCGCCAATAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATaaaa                                                                                                                                        <  2:1008367/150‑1 (MQ=255)
                gATCGCCAATAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAAc                                                                                                                                       >  1:2628434/1‑150 (MQ=255)
                      cAATAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTa                                                                                                                                 <  1:472155/150‑1 (MQ=255)
                       aaTAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTAc                                                                                                                                <  2:2110980/150‑1 (MQ=255)
                        aTAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACa                                                                                                                               >  1:2326308/1‑150 (MQ=255)
                           aCCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAg                                                                                                                            <  2:1655161/150‑1 (MQ=255)
                                  cGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGc                                                                                                                     >  1:1158940/1‑150 (MQ=255)
                                   gTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATAACTTCGTCAATAAAACGCAGTACAAAGATTACGcc                                                                                                                    >  1:4813104/1‑150 (MQ=255)
                                     cTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCAt                                                                                                                  <  2:2727014/150‑1 (MQ=255)
                                      tttGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCAtt                                                                                                                 <  1:970278/150‑1 (MQ=255)
                                       ttGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTc                                                                                                                <  2:4378969/150‑1 (MQ=255)
                                       ttGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTc                                                                                                                <  2:652923/150‑1 (MQ=255)
                                       ttGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTc                                                                                                                <  2:4069768/150‑1 (MQ=255)
                                         gAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCtt                                                                                                              <  1:2902965/150‑1 (MQ=255)
                                          aaaaGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCttt                                                                                                             <  2:944762/150‑1 (MQ=255)
                                           aaaGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTa                                                                                                            >  2:2569155/1‑150 (MQ=255)
                                           aaaGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTa                                                                                                            >  2:149605/1‑150 (MQ=255)
                                            aaGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTAt                                                                                                           <  1:3772771/150‑1 (MQ=255)
                                             aGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATc                                                                                                          >  1:1814774/1‑150 (MQ=255)
                                               cGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgc                                                                                                        >  2:1033070/1‑150 (MQ=255)
                                                gTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgcg                                                                                                       >  1:3576993/1‑150 (MQ=255)
                                                gTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgcg                                                                                                       >  2:869066/1‑150 (MQ=255)
                                                gTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgcg                                                                                                       >  2:2138963/1‑150 (MQ=255)
                                                gTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgcg                                                                                                       <  1:3042601/150‑1 (MQ=255)
                                                  cTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGt                                                                                                     <  2:77491/150‑1 (MQ=255)
                                                  cTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGt                                                                                                     <  2:1170766/150‑1 (MQ=255)
                                                  cTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGt                                                                                                     >  2:1852037/1‑150 (MQ=255)
                                                   tGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCCTTCTTTATCGCGGTa                                                                                                    >  2:3550156/1‑150 (MQ=255)
                                                     ttCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAAc                                                                                                  <  2:2484487/150‑1 (MQ=255)
                                                       ctcCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCa                                                                                                <  1:3946792/150‑1 (MQ=255)
                                                        tcCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCAt                                                                                               <  2:2295495/150‑1 (MQ=255)
                                                           gAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATAACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAg                                                                                            >  2:699428/1‑150 (MQ=255)
                                                               tGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTcgc                                                                                        <  1:1960577/150‑1 (MQ=255)
                                                               tGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATAACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTcgc                                                                                        >  2:3178177/1‑150 (MQ=255)
                                                                 ggTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCgg                                                                                      >  2:3745710/1‑150 (MQ=255)
                                                                   ttATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGgt                                                                                    <  2:1087933/150‑1 (MQ=255)
                                                                     aTGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGgtgt                                                                                  <  2:567922/150‑1 (MQ=255)
                                                                      tGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGtt                                                                                 >  1:3299355/1‑150 (MQ=255)
                                                                           cgGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGaaa                                                                            >  1:4232404/1‑150 (MQ=255)
                                                                           cgGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATAACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGaaa                                                                            >  1:856464/1‑150 (MQ=255)
                                                                             gAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGaaaaa                                                                          <  2:3702700/150‑1 (MQ=255)
                                                                              aGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAAt                                                                         >  1:3203999/1‑150 (MQ=255)
                                                                                  aaaaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATAACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCt                                                                     >  1:2756961/1‑150 (MQ=255)
                                                                                   aaaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTg                                                                    >  1:1019294/1‑150 (MQ=255)
                                                                                    aaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGa                                                                   >  2:3328457/1‑150 (MQ=255)
                                                                                       tATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGc                                                                <  2:4136400/150‑1 (MQ=255)
                                                                                        aTTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCa                                                               >  1:1055802/1‑150 (MQ=255)
                                                                                          tAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATAACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCaaa                                                             >  2:3528299/1‑150 (MQ=255)
                                                                                           aaGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCaaaa                                                            >  1:3806720/1‑150 (MQ=255)
                                                                                             gcgcGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAAcc                                                          <  2:2422143/150‑1 (MQ=255)
                                                                                              cgcgAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCg                                                         <  1:249722/150‑1 (MQ=255)
                                                                                                 gAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCAt                                                      <  1:4992841/150‑1 (MQ=255)
                                                                                                   aTAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGAAAAACCGCATcc                                                    >  2:5055700/1‑150 (MQ=255)
                                                                                                    tAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATccc                                                   <  2:4599702/150‑1 (MQ=255)
                                                                                                     aaCGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATAACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCg                                                  >  1:1666067/1‑150 (MQ=255)
                                                                                                      aCGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGt                                                 <  1:4162294/150‑1 (MQ=255)
                                                                                                        gaagaaCATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTAc                                               <  1:1209197/150‑1 (MQ=255)
                                                                                                        gaagaaCATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTAc                                               >  1:1693810/1‑150 (MQ=255)
                                                                                                        gaagaaCATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTAc                                               >  1:425630/1‑150 (MQ=255)
                                                                                                            aaCATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACaaaa                                           >  2:1142549/1‑150 (MQ=255)
                                                                                                                 gAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCt                                      <  2:3768589/150‑1 (MQ=255)
                                                                                                                    gCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCtggt                                   <  2:152653/150‑1 (MQ=255)
                                                                                                                     cTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCtggtg                                  <  1:2959379/150‑1 (MQ=255)
                                                                                                                      tGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCtggtgg                                 >  1:4346462/1‑150 (MQ=255)
                                                                                                                        aGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTa                               >  2:1266509/1‑150 (MQ=255)
                                                                                                                          cgcgTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTAcg                             <  1:3790242/150‑1 (MQ=255)
                                                                                                                          cgcgTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTAcg                             >  2:1970165/1‑150 (MQ=255)
                                                                                                                           gcgTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTAcgt                            <  1:3704065/150‑1 (MQ=255)
                                                                                                                              ttACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTAcgtcgt                         <  2:1814774/150‑1 (MQ=255)
                                                                                                                                aCTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGttt                       >  1:222372/1‑150 (MQ=255)
                                                                                                                                 cTGGCGAGCTGATCGCCCATAACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGtttt                      <  1:1112384/150‑1 (MQ=255)
                                                                                                                                   ggCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTtc                    <  2:4468477/150‑1 (MQ=255)
                                                                                                                                       aGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGAACAAAATATCTGGTGGTACGTCGTTTTtctctc                >  1:1394194/1‑150 (MQ=255)
                                                                                                                                        gCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCg               >  2:4795365/1‑150 (MQ=255)
                                                                                                                                         cTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGt              >  2:4021412/1‑150 (MQ=255)
                                                                                                                                          tGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTc             >  1:3421485/1‑150 (MQ=255)
                                                                                                                                            aTCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCt           >  1:3365147/1‑150 (MQ=255)
                                                                                                                                            aTCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCt           >  1:5100073/1‑150 (MQ=255)
                                                                                                                                            aTCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCt           >  2:4998585/1‑150 (MQ=255)
                                                                                                                                               gCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGaa        <  2:3812564/150‑1 (MQ=255)
                                                                                                                                                 ccATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAAt      >  1:620713/1‑150 (MQ=255)
                                                                                                                                                 ccATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAAt      >  2:2940368/1‑150 (MQ=255)
                                                                                                                                                  cATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATc     >  2:3971647/1‑150 (MQ=255)
                                                                                                                                                  cATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATc     >  2:3971899/1‑150 (MQ=255)
                                                                                                                                                   aTCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCa    >  2:3603451/1‑150 (MQ=255)
                                                                                                                                                    tCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCaa   >  1:4446763/1‑150 (MQ=255)
                                                                                                                                                    tCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCaa   >  1:1359370/1‑150 (MQ=255)
                                                                                                                                                     cACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAg  <  1:362030/150‑1 (MQ=255)
                                                                                                                                                     cACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAg  >  2:4197301/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGCGGCGAACATCCTGATCGCCAATAACCCGCACGTCTTTGAAAAGCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAG  >  CP009273/3954890‑3955188

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