Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,191,676 C→T 31.1% H366H (CAC→CAT icd → e14 prophage; isocitrate dehydrogenase, specific for NADP+

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,191,6760CT31.1% 93.3 / 49.3 74H366H (CAC→CATicde14 prophage; isocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base C (29/20);  new base T (14/8);  total (46/28)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

GCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAAC  >  CP009273/1191531‑1191823
                                                                                                                                                 |                                                                                                                                                   
gCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTCCTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATTATGCTGCGCCAAATg                                                                                                                                                  >  1:757326/1‑149 (MQ=255)
       cggtatcgCCCCTGGTGCAAACATCGGTGACGTATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGa                                                                                                                                           <  1:759119/149‑1 (MQ=255)
        ggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATggg                                                                                                                                                >  2:918215/1‑143 (MQ=255)
         gtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCGCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTAACTCCGCTGAGATGATGCTGCACCACATGGGt                                                                                                                                               >  2:203763/1‑143 (MQ=255)
         gtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCCCTATTATTATCTCCGCTGAGATGATGCTGCGCCATATGGGt                                                                                                                                               >  2:131429/1‑143 (MQ=255)
         gtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGGACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGAcc                                                                                                                                         >  1:1223743/1‑149 (MQ=255)
         gtatcgCCCCTGGTGCAAACATCGGTGACCCGTGAGCCTTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCAACTGGGt                                                                                                                                               >  2:517035/1‑143 (MQ=255)
            tcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGTCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACCTGGGTTGGACCGac                                                                                                                                      >  1:692780/1‑148 (MQ=255)
            tcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGaa                                                                                                                                      <  1:270431/149‑1 (MQ=255)
            tcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCGCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTgg                                                                                                                                            >  2:1121326/1‑143 (MQ=255)
             cgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACAGTACTGCGCCGAAATATGCCGGTCAGTACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGa                                                                                                                                           <  2:722805/143‑1 (MQ=255)
              gCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGGTGATGCTGCGCCACAGGGGTTGGAc                                                                                                                                          >  2:481547/1‑143 (MQ=255)
              gCCCCTGGTGCAAACATCGGGGACGAATGCGCCCTGTTGGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACACAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATAGGGGGGGGAc                                                                                                                                          >  2:923476/1‑143 (MQ=255)
                cccTGGTGCAAACCTCGGTGCGGAATGCGCCCTGTTTGAAGCCACCCGCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGCAAATCCTGGCTCTATTATTCTCGCCGCTGAGATGATGCTGCCCCACATGGGTTGGACCg                                                                                                                                        <  2:841898/143‑1 (MQ=255)
                 ccTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCAAATGGGTTGGACCGa                                                                                                                                       >  2:283312/1‑143 (MQ=255)
                     gTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGAc                                                                                                                             <  1:190836/149‑1 (MQ=255)
                        cAAACAGCGGTGACGAATGCGCCCTGTTGGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGACTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCt                                                                                                                                <  2:913639/143‑1 (MQ=255)
                          aaCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAcctgat                                                                                                                        >  1:1174398/1‑144 (MQ=255)
                             aTCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTGTTATTCTCTCCGCTGAGATGATGCTGCGCCATCTGGGTTGGACTGAAGCCGCTggcct                                                                                                                           >  2:113755/1‑139 (MQ=255)
                               cGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTaa                                                                                                                         >  2:282729/1‑143 (MQ=255)
                                 gTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTaaa                                                                                                                 >  1:909223/1‑149 (MQ=255)
                                 gTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAAtt                                                                                                                       <  2:1110244/143‑1 (MQ=255)
                                 gTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAAtt                                                                                                                       <  2:1056307/143‑1 (MQ=255)
                                   gACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTCAATCCTGGCTGTATTATTCTCTCCGCTGAGCTGGTGCTGCGCCACCTGGGTTGCCCCGAAGCGGTAGACGTAAttcg                                                                                                                     >  2:702569/1‑141 (MQ=255)
                                   gACGAATGAGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAAGCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGt                                                                                                                     <  2:61890/143‑1 (MQ=38)
                                       aaggcgcCCTTTTTGAAGCCACCCACGGCACTGCGCAGAAACATGCCGGTCAGGACAAAGTAAATCCTGGCTCTACTATTCTCGCCGCGGAGATGATGCTGCTCCACATGGGTTGTACCGAAGCGGCTGACTTAATTGTTaaa                                                                                                                 <  2:66222/140‑1 (MQ=255)
                                         tGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGa                                                                                                         <  1:73849/149‑1 (MQ=255)
                                          gcgcCCTGCTTTAAGCCACCCCCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGt                                                                                                              <  2:58185/143‑1 (MQ=255)
                                                tGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGaa                                                                                                        <  2:278469/143‑1 (MQ=255)
                                                 gTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAg                                                                                                       <  2:787153/143‑1 (MQ=37)
                                                  tttGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAgg                                                                                                      >  2:1044861/1‑143 (MQ=255)
                                                    tGAAGACACCCACGGTACTGCGCCGAAATATGCCGGTTATGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgcg                                                                                                    <  2:197895/143‑1 (MQ=25)
                                                        gCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc                                                                                          >  1:463882/1‑149 (MQ=255)
                                                             ccccGGTACTGCGCCCGAATATGCCGGTCAGGACAAAGTAACTCCTGGCTCGATTATTCTCTCCGCTGAGATGATGCTGCGGCCCATCGGGTGGACCGACGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaaa                                                                                     <  1:310005/146‑1 (MQ=37)
                                                                  gTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTCTTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAAATGTTAAAGGTATGGAAGGCGCAATCAACGCGaaa                                                                                      >  2:513074/1‑143 (MQ=255)
                                                                   tACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa                                                                               >  1:664639/1‑149 (MQ=255)
                                                                           cGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGAc                                                                       >  1:1021579/1‑149 (MQ=17)
                                                                                    ccGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGt                                                              >  1:1120615/1‑149 (MQ=14)
                                                                                     cGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGAACTGAAGCGGCTGACCTGATTGTTAAAGGTGTGGAAGGCGCAATCAATGCCAAGACCGTAACTTGTGACTTCGAACgac                                                             >  1:119041/1‑147 (MQ=14)
                                                                                     cGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCg                                                                   >  2:1065359/1‑143 (MQ=255)
                                                                                        tCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGTTGCTGCGCCGCATGGGTTGGACGGAAGCGGGTGAGTTGATTGTTAAAGGTATGGAAGGCGCAATCAAAGCGAAAACCGTAGCGTCTGATTTTGAGc                                                                >  2:322706/1‑143 (MQ=21)
                                                                                        tCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTga                                                          >  1:1075035/1‑149 (MQ=255)
                                                                                         cAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgat                                                         <  1:323611/149‑1 (MQ=255)
                                                                                           ggACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTc                                                             >  2:344124/1‑143 (MQ=255)
                                                                                           ggACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAACGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTc                                                             >  2:549664/1‑143 (MQ=255)
                                                                                           ggACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGATTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTc                                                             >  2:840939/1‑143 (MQ=255)
                                                                                            gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTCGACCGAAGCGGCTGAGTGAATGGTTAAAGGTATGGAAGGCGCAAACGACGCGAAAAGCGTAACCTATGACTTCGAGCGTCt                                                            >  2:604/1‑143 (MQ=37)
                                                                                            gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATAGGGTGGTCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCCATCAACGCGAAACCCGTAACCTATGACTTCGAGCGTCTgatgga                                                      >  1:1083712/1‑149 (MQ=37)
                                                                                                aaGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAATACCGTAACCTATGACTTCGAGCGTCTgatg                                                        >  2:174325/1‑143 (MQ=255)
                                                                                                 aGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGgcg                                                 <  1:1265474/149‑1 (MQ=11)
                                                                                                    aaaTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTa                                              >  1:269484/1‑149 (MQ=11)
                                                                                                    aaaTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAACGGCTAACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGATTTCGAACGTCTGATGGAAGGCGCTa                                              >  1:797059/1‑149 (MQ=11)
                                                                                                     aaTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATGTGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGatggaagg                                                   <  2:764645/143‑1 (MQ=1)
                                                                                                          tGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATAAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAActgc                                        <  1:1183783/149‑1 (MQ=255)
                                                                                                            gCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATCTGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGctgctg                                      >  1:1148655/1‑149 (MQ=14)
                                                                                                             ctctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAAc                                           <  2:897479/143‑1 (MQ=255)
                                                                                                              tctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGACCGTCTGTTGGAAGGCGCTAAGct                                          >  2:570638/1‑143 (MQ=11)
                                                                                                              tctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGct                                          >  2:514624/1‑143 (MQ=14)
                                                                                                              gctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTTAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGaa                                    >  1:663760/2‑149 (MQ=14)
                                                                                                                tattatTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAAt                                  >  1:124315/1‑149 (MQ=255)
                                                                                                                tattatTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAAt                                  >  1:1147761/1‑149 (MQ=255)
                                                                                                                   tattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGCCTTCGAGCGTCTGATGGAAGGCGCTAAACTGCTGa                                     >  2:474746/1‑143 (MQ=255)
                                                                                                                        tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCagag                          >  1:764259/1‑149 (MQ=255)
                                                                                                                            cGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTg                      <  1:1214061/149‑1 (MQ=255)
                                                                                                                             gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTgg                     <  1:743923/149‑1 (MQ=255)
                                                                                                                             gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGGATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGCCTTCGAGCGTCTGATGGATGGCGCTAACCTGCTGAACTGTTCaga                           >  2:370352/1‑143 (MQ=255)
                                                                                                                                gagaTGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGttacgaggt                        >  2:100895/1‑136 (MQ=255)
                                                                                                                                gagaTGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGa                  >  1:1264393/1‑149 (MQ=255)
                                                                                                                                      agtctgcCCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGa                  <  2:789271/140‑1 (MQ=255)
                                                                                                                                        gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatca          <  1:282381/149‑1 (MQ=255)
                                                                                                                                         ctgcGCCACTGGGGTTGGACCGAAGTGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcgc               <  2:145973/143‑1 (MQ=255)
                                                                                                                                           gcgcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCg       >  1:1029371/1‑149 (MQ=255)
                                                                                                                                             gcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGTAAGGCGCGATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaa     <  1:961858/149‑1 (MQ=255)
                                                                                                                                             gcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaa     <  1:238035/149‑1 (MQ=255)
                                                                                                                                                acaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGCTTGGTGACGCGctcctc        >  2:255908/1‑143 (MQ=255)
                                                                                                                                                acaAGGGTTGGACCGAAGCGGCTTACTTGTTTTTTAAAGGTATTGAAGGAGCAATCAACGCGAAAACCGTAGCCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAAc  <  1:233749/149‑1 (MQ=17)
                                                                                                                                                 |                                                                                                                                                   
GCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAAC  >  CP009273/1191531‑1191823

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