Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,204,161 G→A 34.4% E8K (GAA→AAA)  rpsU → 30S ribosomal subunit protein S21

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,204,1610GA34.4% 51.6 / 48.9 61E8K (GAA→AAA) rpsU30S ribosomal subunit protein S21
Reads supporting (aligned to +/- strand):  ref base G (11/29);  new base A (9/12);  total (20/41)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-01

AAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGT  >  CP009273/3204026‑3204304
                                                                                                                                       |                                                                                                                                                
aaaGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGa                                                                                                                                           <  2:1019436/143‑1 (MQ=255)
      cgcgGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGGGAGAGGCACATGCCGGTAATTAAAGTACGGGAAAACGAGCCGTTCGACGt                                                                                                                               >  1:736927/1‑149 (MQ=255)
        cgGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCAGTTCGACATAg                                                                                                                             >  1:876978/1‑149 (MQ=255)
         ggcgCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGa                                                                                                                                  <  2:590944/140‑1 (MQ=255)
          gTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGAc                                                                                                                                 <  2:735483/143‑1 (MQ=255)
              gttACAAAACAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACACGCCGGCAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAg                                                                                                                             <  2:69532/142‑1 (MQ=255)
               ttACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGc                                                                                                                            <  2:621871/143‑1 (MQ=255)
                     agcagcagcAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCg                                                                                                                      >  2:146526/1‑143 (MQ=255)
                      gcagcagcaATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTc                                                                                                               >  1:557090/1‑149 (MQ=255)
                          cagcaATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGAAATTAAAGTACGAAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCt                                                                                                                 >  2:303526/1‑143 (MQ=255)
                                    aGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTc                                                                                                       >  2:1067006/1‑143 (MQ=255)
                                    aGTAAAACTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCCTCAAACGtca                                                                                                       >  2:763650/1‑141 (MQ=255)
                                     gTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGaa                                                                                               <  1:1176751/149‑1 (MQ=255)
                                       aagatttggctgtaTTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAATCCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGaaaa                                                                                             <  1:166996/136‑1 (MQ=255)
                                        aacttCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGaaaaa                                                                                            <  1:1107056/146‑1 (MQ=255)
                                         aaTTCCGCACCATTTTGAAATAAGCTGGCGTTGAGGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTCCCCCGGAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCg                                                                                                 <  2:246186/143‑1 (MQ=255)
                                          aTTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGGGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGa                                                                                                <  2:278142/143‑1 (MQ=255)
                                            tCCGCACCATTTTGAAATAAGCCGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAg                                                                                        <  1:609809/149‑1 (MQ=255)
                                             ccGCAACATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCTAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAgg                                                                                       <  1:1221402/149‑1 (MQ=255)
                                              cGCACCATTTTGAAATAAGCTGGCTTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAATTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGaaaaa                                                                                            <  2:842588/143‑1 (MQ=255)
                                              cGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGgt                                                                                      <  1:1080494/149‑1 (MQ=255)
                                                cACCATTTTGAAATACGCGGGCGTTGGTCCCAGCGGCAATCCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGc                                                                                          <  2:104819/143‑1 (MQ=255)
                                                  ccATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAg                                                                                        >  2:1126224/1‑143 (MQ=255)
                                                     ttttTAAATAAGCTGGCGTTGATGCCAGCGGCAAACCAAAATAAACAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGCGCCGTTCGACGTAGCTCTGCCTCGCTTCAAGCAAAC‑AAACAAAAAAGCAGgtg                                                                                     <  2:597222/143‑1 (MQ=255)
                                                     ttttGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTACGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGgtg                                                                                     >  2:373027/1‑143 (MQ=255)
                                                           aaTAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGCGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTgg                                                                               >  2:329147/1‑143 (MQ=255)
                                                           aaTAAGCCGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAg                                                                         <  1:231450/149‑1 (MQ=255)
                                                            aTAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAACCCAAGGTGAGAGGTACATGCAGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGc                                                                              <  2:862/143‑1 (MQ=255)
                                                                    gCGTTGATGCGAGCGGCAAACCGAATTAATCAAGGGTGAGCGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCCAAAAAGCAGGTGTTCTGGCAGAAGGtc                                                                      <  2:1053210/143‑1 (MQ=255)
                                                                     cGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgtcg                                                               <  1:153161/149‑1 (MQ=255)
                                                                      gTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTcgtcgtcgt                                                              <  1:560621/149‑1 (MQ=255)
                                                                         ggtGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgt                                                                 <  2:762356/141‑1 (MQ=255)
                                                                         gATGCCAGCGGCGAACCGAATGAAGCAAAGGTGAGAGTCGCATGCCGGGAAGTAAGGTACGTAAAAACGAACGATTCTACTAGGGTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTCCTGGCGGAAGTtggtcgc                                                                 >  2:293508/1‑142 (MQ=255)
                                                                          aTGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGAAAGTTCGTCGTCGTGAGt                                                          >  1:627807/1‑149 (MQ=255)
                                                                              cAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGa                                                            <  2:498026/143‑1 (MQ=255)
                                                                                 cGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGGGAGTTCTATGa                                                   >  1:1126171/1‑149 (MQ=255)
                                                                                      aaCCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAAc                                              <  1:605357/149‑1 (MQ=255)
                                                                                      aaCCGAATTAACCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCCACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAAc                                              <  1:203577/149‑1 (MQ=255)
                                                                                       aCCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAAcc                                             >  1:264506/1‑149 (MQ=255)
                                                                                           aaTTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaaaa                                               <  2:661303/143‑1 (MQ=255)
                                                                                            aTTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTa                                        >  1:747766/1‑149 (MQ=255)
                                                                                            aTTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGGGAGTTCTAGGAAAAACCGACTa                                        >  1:683455/1‑149 (MQ=255)
                                                                                              tAAACAAAGGTGAGAGGCACATCCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTAcc                                      <  1:405928/149‑1 (MQ=255)
                                                                                                 tCAAAGGTGAGAGGCACATGCCGGTAATTAAAGCACGTAAAACCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGaa                                   <  1:1084911/149‑1 (MQ=255)
                                                                                                 gccgaGGTGGGAGGCNCATGTCGGTAATGAAAGTAGGGACAAGCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CGGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGaa                                   <  1:1160364/145‑1 (MQ=255)
                                                                                                  cAAAGGTCAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAAc                                  <  1:663875/149‑1 (MQ=255)
                                                                                                  cAAAGGGGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAAc                                  <  1:516597/149‑1 (MQ=255)
                                                                                                                acaTGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTACGCGCGCTaaa                    >  1:499449/1‑149 (MQ=255)
                                                                                                                      cGGTAATTAAAGTAGGGGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGGGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCt              <  1:1090099/149‑1 (MQ=255)
                                                                                                                      agttAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCGCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGAAAGAAAAAAAAAAAACCGACTACCGAACGTAAGCGCGCTaaa                    <  2:412058/140‑1 (MQ=255)
                                                                                                                       ggTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAACAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGCCTACCGAACGTAAGCGCGCTAAAg                   >  2:447429/1‑143 (MQ=255)
                                                                                                                        gTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGATCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTCCCGAACGTAAGCGCGCTAAAGCTTCTGc            <  1:849540/149‑1 (MQ=255)
                                                                                                                        gTAATTAAAGTACGTAAAAACGAGCCGTTCGACCTCGGCCTGCGTTGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTACTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGc            >  1:634598/1‑149 (MQ=255)
                                                                                                                             tAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCtc             >  2:217936/1‑142 (MQ=255)
                                                                                                                              aaaGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGaa      <  1:903688/149‑1 (MQ=255)
                                                                                                                               aaGTCCGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGaaa     <  1:1207117/149‑1 (MQ=255)
                                                                                                                               aaGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCa           <  2:230520/143‑1 (MQ=255)
                                                                                                                               aaGTACGTAAAAACGAGCCGTTTGTCTTAGTTCTGCGTCGCTTCTTG‑GTTCAAAGTAAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGaaa     <  1:364496/149‑1 (MQ=255)
                                                                                                                               aaGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTAT‑TAACTAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGaaa     <  1:364565/149‑1 (MQ=255)
                                                                                                                                aGTCCGTGAAAACGAGCCGTTCGACGTAGCTTTGCGTCGCTTCAAGCGTTC‑CTTCAAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGGTCCCTGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAAc    <  1:186545/149‑1 (MQ=255)
                                                                                                                                  tACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGt  <  1:567969/149‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                
AAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTC‑CTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGT  >  CP009273/3204026‑3204304

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