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

AAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGT  >  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)
                                     gTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGaa                                                                                               <  1:1176751/149‑1 (MQ=255)
                                       aagatttggctgtaTTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAATCCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGaaaa                                                                                             <  1:166996/136‑1 (MQ=255)
                                        aacttCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGaaaaa                                                                                            <  1:1107056/146‑1 (MQ=255)
                                         aaTTCCGCACCATTTTGAAATAAGCTGGCGTTGAGGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTCCCCCGGAGCTCTGCGTCGCTTCAAGCGTTCCTGCg                                                                                                 <  2:246186/143‑1 (MQ=255)
                                          aTTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGGGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGa                                                                                                <  2:278142/143‑1 (MQ=255)
                                            tCCGCACCATTTTGAAATAAGCCGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAg                                                                                        <  1:609809/149‑1 (MQ=255)
                                             ccGCAACATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCTAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAgg                                                                                       <  1:1221402/149‑1 (MQ=255)
                                              cGCACCATTTTGAAATAAGCTGGCTTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAATTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGaaaaa                                                                                            <  2:842588/143‑1 (MQ=255)
                                              cGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgt                                                                                      <  1:1080494/149‑1 (MQ=255)
                                                cACCATTTTGAAATACGCGGGCGTTGGTCCCAGCGGCAATCCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGc                                                                                          <  2:104819/143‑1 (MQ=255)
                                                  ccATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAg                                                                                        >  2:1126224/1‑143 (MQ=255)
                                                     ttttTAAATAAGCTGGCGTTGATGCCAGCGGCAAACCAAAATAAACAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGCGCCGTTCGACGTAGCTCTGCCTCGCTTCAAGCAAACAAACAAAAAAGCAGgtg                                                                                     <  2:597222/143‑1 (MQ=255)
                                                     ttttGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTACGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtg                                                                                     >  2:373027/1‑143 (MQ=255)
                                                           aaTAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGCGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTgg                                                                               >  2:329147/1‑143 (MQ=255)
                                                           aaTAAGCCGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAg                                                                         <  1:231450/149‑1 (MQ=255)
                                                            aTAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAACCCAAGGTGAGAGGTACATGCAGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGc                                                                              <  2:862/143‑1 (MQ=255)
                                                                    gCGTTGATGCGAGCGGCAAACCGAATTAATCAAGGGTGAGCGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCCAAAAAGCAGGTGTTCTGGCAGAAGGtc                                                                      <  2:1053210/143‑1 (MQ=255)
                                                                     cGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgtcg                                                               <  1:153161/149‑1 (MQ=255)
                                                                      gTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTcgtcgtcgt                                                              <  1:560621/149‑1 (MQ=255)
                                                                         ggtGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgt                                                                 <  2:762356/141‑1 (MQ=255)
                                                                         gATGCCAGCGGCGAACCGAATGAAGCAAAGGTGAGAGTCGCATGCCGGGAAGTAAGGTACGTAAAAACGAACGATTCTACTAGGGTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTCCTGGCGGAAGTtggtcgc                                                                 >  2:293508/1‑142 (MQ=255)
                                                                          aTGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGAAAGTTCGTCGTCGTGAGt                                                          >  1:627807/1‑149 (MQ=255)
                                                                              cAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGa                                                            <  2:498026/143‑1 (MQ=255)
                                                                                 cGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGGGAGTTCTATGa                                                   >  1:1126171/1‑149 (MQ=255)
                                                                                      aaCCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAAc                                              <  1:605357/149‑1 (MQ=255)
                                                                                      aaCCGAATTAACCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCCACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAAc                                              <  1:203577/149‑1 (MQ=255)
                                                                                       aCCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAAcc                                             >  1:264506/1‑149 (MQ=255)
                                                                                           aaTTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaaaa                                               <  2:661303/143‑1 (MQ=255)
                                                                                            aTTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTa                                        >  1:747766/1‑149 (MQ=255)
                                                                                            aTTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGGGAGTTCTAGGAAAAACCGACTa                                        >  1:683455/1‑149 (MQ=255)
                                                                                              tAAACAAAGGTGAGAGGCACATCCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTAcc                                      <  1:405928/149‑1 (MQ=255)
                                                                                                 tCAAAGGTGAGAGGCACATGCCGGTAATTAAAGCACGTAAAACCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGaa                                   <  1:1084911/149‑1 (MQ=255)
                                                                                                 gccgaGGTGGGAGGCNCATGTCGGTAATGAAAGTAGGGACAAGCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCGGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGaa                                   <  1:1160364/145‑1 (MQ=255)
                                                                                                  cAAAGGTCAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAAc                                  <  1:663875/149‑1 (MQ=255)
                                                                                                  cAAAGGGGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAAc                                  <  1:516597/149‑1 (MQ=255)
                                                                                                                acaTGCCGGTAATTAAAGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTACGCGCGCTaaa                    >  1:499449/1‑149 (MQ=255)
                                                                                                                      cGGTAATTAAAGTAGGGGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCt              <  1:1090099/149‑1 (MQ=255)
                                                                                                                      agttAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCGCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGAAAGAAAAAAAAAAAACCGACTACCGAACGTAAGCGCGCTaaa                    <  2:412058/140‑1 (MQ=255)
                                                                                                                       ggTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAACAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGCCTACCGAACGTAAGCGCGCTAAAg                   >  2:447429/1‑143 (MQ=255)
                                                                                                                        gTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGATCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTCCCGAACGTAAGCGCGCTAAAGCTTCTGc            <  1:849540/149‑1 (MQ=255)
                                                                                                                        gTAATTAAAGTACGTAAAAACGAGCCGTTCGACCTCGGCCTGCGTTGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTACTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGc            >  1:634598/1‑149 (MQ=255)
                                                                                                                             tAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCtc             >  2:217936/1‑142 (MQ=255)
                                                                                                                              aaaGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGaa      <  1:903688/149‑1 (MQ=255)
                                                                                                                               aaGTCCGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGaaa     <  1:1207117/149‑1 (MQ=255)
                                                                                                                               aaGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCa           <  2:230520/143‑1 (MQ=255)
                                                                                                                               aaGTACGTAAAAACGAGCCGTTTGTCTTAGTTCTGCGTCGCTTCTTGGTTCAAAGTAAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGaaa     <  1:364496/149‑1 (MQ=255)
                                                                                                                               aaGTACGTAAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTATTAACTAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGaaa     <  1:364565/149‑1 (MQ=255)
                                                                                                                                aGTCCGTGAAAACGAGCCGTTCGACGTAGCTTTGCGTCGCTTCAAGCGTTCCTTCAAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGGTCCCTGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAAc    <  1:186545/149‑1 (MQ=255)
                                                                                                                                  tACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGt  <  1:567969/149‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                
AAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGT  >  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: