Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,960,102 G→T 31.4% R109L (CGC→CTC)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,960,1020GT31.4% 81.5 / 44.1 69R109L (CGC→CTC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (30/17);  new base T (14/8);  total (44/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACT  >  CP009273/3959957‑3960250
                                                                                                                                                 |                                                                                                                                                    
ggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAACCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTa                                                                                                                                                   <  1:667561/149‑1 (MQ=255)
   tttGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAttt                                                                                                                                                >  1:1089068/1‑149 (MQ=255)
   tttGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAACGTGAACGCTAttt                                                                                                                                                >  1:216843/1‑149 (MQ=255)
           cctccGTTCCGCAGACAGCCCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTCGCCCAATCCGCCGTTTCAACCTCCGCCCTGGTGACACCATCCCCGGTCACTGTCGCCCGCGGTAAGAAGGTGGTCGCTATTTTg                                                                                                                                              <  2:659332/143‑1 (MQ=255)
            ctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                                                                                             >  2:704154/1‑143 (MQ=255)
            ctcacttCCGCAGACAGCCCCCACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCCACCTCCGCACTGGTGATACCACCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgct                                                                                                                                       <  1:132774/144‑1 (MQ=255)
             tccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGACACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGACCGCTATTTTgcg                                                                                                                                            >  2:572410/1‑143 (MQ=255)
             tccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGACACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgcg                                                                                                                                            >  2:1044679/1‑143 (MQ=255)
               cGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgct                                                                                                                                          >  2:914182/1‑143 (MQ=255)
               cGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGGGAACGCTGTTTTGCgct                                                                                                                                          >  2:913985/1‑143 (MQ=255)
                    cGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTCCCACCTCCGCACTGGTGATACCATCTCTGGTGAGATGCGCCCGCCGAAAGAAGGTGCACGCTCTTTGGGGctgctgg                                                                                                                                     >  2:376265/1‑142 (MQ=255)
                      cAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAAc                                                                                                                             >  1:111363/1‑149 (MQ=255)
                      cAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAAc                                                                                                                             >  1:415320/1‑149 (MQ=255)
                         aCAGCTCCCTCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGaa                                                                                                                          <  1:25284/149‑1 (MQ=255)
                             cTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTCGGCTGCTGAAAGTTAACg                                                                                                                            >  2:542226/1‑143 (MQ=255)
                                cTACCGCGCCGGTCCTGTTGATATCTACGTTTCCTCTAGCCAAATCCGCCGTTTCAACCTCCGCAATGGGGATTTCATCTTTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAg                                                                                                                         <  2:403128/143‑1 (MQ=255)
                                 tACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGAACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAAcgt                                                                                                                  >  1:1010461/1‑147 (MQ=255)
                                    cccGCCGGTCCTGATGACCTCTACGTCTCCCCTAGCCAACTCCGCCGTCTCCACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATGTTGCGCTGCTGAAAGTTAACGAAGTTaa                                                                                                                     <  2:1062293/141‑1 (MQ=255)
                                    cTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCACATCCGCCGTTTCCACCTCACCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCCCGTCTGCGCTGGCGAAAGGTACCgcggggaa                                                                                                                     >  2:216600/1‑136 (MQ=255)
                                    cTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAGCGTAGTTaa                                                                                                                     <  2:296536/143‑1 (MQ=255)
                                    cTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTaa                                                                                                                     >  2:1080811/1‑143 (MQ=255)
                                    cTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTaa                                                                                                                     <  2:363457/143‑1 (MQ=255)
                                     tCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                                                                                              >  1:1282987/1‑149 (MQ=255)
                                         cGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTAATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAACGTTAACGAAGTTAACTTCGACAAAc                                                                                                          >  1:410149/1‑149 (MQ=255)
                                                      aTCTACTGTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCg                                                                                             <  1:235390/149‑1 (MQ=255)
                                                      aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCg                                                                                             <  1:1291535/149‑1 (MQ=255)
                                                        cTACGTTTCCCCTAGCCAAACCCGCCGTTTCAACCTCCGCACTGGTGATACCATCCCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGca                                                                                           <  1:1153447/149‑1 (MQ=255)
                                                            gTTGCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCg                                                                                             <  2:497611/143‑1 (MQ=255)
                                                              ttCCCCGAGCCAAAGCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGca                                                                                           <  2:206460/143‑1 (MQ=255)
                                                                 cccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaaca                                                                                        >  2:497671/1‑143 (MQ=255)
                                                                 cccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGAGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAGACCTGACAACGCCCGcacca                                                                                        >  2:285269/1‑143 (MQ=255)
                                                                  ccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGACGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCct                                                                                 >  1:963324/1‑149 (MQ=255)
                                                                   cTAGCCAAATCCGCCGTTTCAACCTCCGCTCTGGTGATACCATCTCTGGTAAGAATCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAATCCTGaaaaaaaaaaaaaaaaa                                                                                      <  2:972145/143‑14 (MQ=255)
                                                                   cTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc                                                                                >  1:494476/1‑149 (MQ=255)
                                                                     aGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCAGGCAAGATCCGCGCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt                                                                                    >  2:857844/1‑143 (MQ=255)
                                                                        cAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCTGACAGAAGGTGAATGATATTTTGCGCTGCTGACAGTTAACGAAGTTAAGTTCGACAGACCTGAAAGCGCCCGCAACaaaagctt                                                                                 >  2:62914/1‑139 (MQ=255)
                                                                             ccgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACt                                                                      <  1:606680/149‑1 (MQ=255)
                                                                               gccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATACGCCCGCCGAAAGAAGGTGAACTCTATATTGCGCTGCTGAAAGTTAACGAAGGTAACTTCGACAAACATGAAAAACCCCACGCGAAAATCCTCTTTgat                                                                          >  2:285937/1‑142 (MQ=255)
                                                                                ctgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTaa                                                                   <  1:906591/147‑1 (MQ=255)
                                                                                          cTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTGTAGAACTTAACCCGGCTGcc                                                         >  1:749970/1‑148 (MQ=255)
                                                                                          cTCCGCACTGGTGATACCATCTCCGGTAAGATTCGCCCGCCGCAAGAAGGGGAACTCTATTTTGCGCTGCTGCAAGTTAACGAAGTGAACTTCGACAAAGCAGAAACCGCCCGTAACAAAATACTCTTTGAGAACTTAAcccc                                                               >  2:777912/1‑143 (MQ=255)
                                                                                          cTCCGCACTGGTGATACCATCTCCGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCa                                                         >  1:686264/1‑149 (MQ=255)
                                                                                            ccGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGGAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCGACAAAATCCTCTTTGCGAACTTAACCCCGc                                                             >  2:158865/1‑143 (MQ=255)
                                                                                               cACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGc                                                          >  2:16402/1‑143 (MQ=255)
                                                                                                   ggTGAGCCCAGCTCTGGTAAGAGTCGCCCGCCGAAAGTAGGTGAACGCTATTTTGCGCTTCTGAAAGTTACCGAAGTTAACTTCGACAAACCTGAAACCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGTACGc                                                      <  2:273474/143‑1 (MQ=255)
                                                                                                    gTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGGTAACGAAGTTAAGTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTCAGAACTTAAACCCGCTGCACGCa                                                     >  2:115575/1‑143 (MQ=255)
                                                                                                     tGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACtctc                                              <  1:1215924/149‑1 (MQ=255)
                                                                                                     ggATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACtctc                                              <  1:140875/148‑1 (MQ=255)
                                                                                                        tACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTc                                           >  1:1025790/1‑149 (MQ=255)
                                                                                                        tACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGGTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTc                                           >  1:895824/1‑149 (MQ=255)
                                                                                                           cATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGc                                        <  1:550088/149‑1 (MQ=255)
                                                                                                            aTCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAACGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTGTGAGAACTGAACCCCGCTGCACGCCAACTCTCg                                             >  2:1038906/1‑143 (MQ=255)
                                                                                                            aTCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGGTAACGAAGTTAACTTCGACAAACCTGACGTTTTCTGCAACAAAATCCTCTTTGAGAACTTAACCCCCCTGCACGCAAACTCTCGTCTGCg                                       >  1:1274959/1‑149 (MQ=255)
                                                                                                                   gTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTCTGGaa                                >  1:238245/1‑149 (MQ=255)
                                                                                                                   gTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGaa                                >  1:879532/1‑149 (MQ=255)
                                                                                                                    tAAGATTCGCCCGCCGAGAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGACGTTAACTTCGACAAACCTGAAAACGCCCGCAACCAAATCCTCTTTGAGAACGTCACCGCGCTGCCAGCAAACTCTCGTCTGCgga                                     >  2:1061666/1‑141 (MQ=255)
                                                                                                                    tAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTa                                     >  2:893248/1‑143 (MQ=255)
                                                                                                                    tAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGACTGCGTATGGAAc                               >  1:930473/1‑149 (MQ=255)
                                                                                                                    tAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACGTAACCCCGCTGCACGCAAACTCTCGTCTGCGTa                                     >  2:572988/1‑143 (MQ=255)
                                                                                                                      aGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAGACGCCCGCAACAAAATCCTCTTTGAGAACTTCACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGACCGt                             >  1:392525/1‑149 (MQ=255)
                                                                                                                        aTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACCAAATCCTCTTTGAGAACTTCACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGa                                 >  2:26422/1‑143 (MQ=255)
                                                                                                                        aTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGAATGGa                                 >  2:1066782/1‑143 (MQ=255)
                                                                                                                            ccccgccGGAAGAAGGTGAACTCGATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGt                             <  2:214528/142‑1 (MQ=255)
                                                                                                                                gccgAAAGAAGGTGAACGCTATTTTACGCTGCTAAAAGCTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGtt                   <  1:241379/149‑1 (MQ=255)
                                                                                                                                 ccgAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTc                  <  1:813426/149‑1 (MQ=255)
                                                                                                                                           gTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTGCTACTGAAGAt        >  1:814520/1‑149 (MQ=255)
                                                                                                                                             gAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGAttt      <  1:634074/149‑1 (MQ=255)
                                                                                                                                              aaCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGaa           <  2:523957/143‑1 (MQ=255)
                                                                                                                                              aaCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTa     >  1:834624/1‑149 (MQ=255)
                                                                                                                                                 gCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAAca  >  1:431799/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACT  >  CP009273/3959957‑3960250

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: ‑