Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 2,863,520 T→A 5.7% T289S (ACG→TCG)  U712_14775 ← Protein AroA(G)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,863,5200TA 5.7% 219.5 / 2.7 70T289S (ACG→TCG) U712_14775Protein AroA(G)
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base A (2/2);  total (34/36)
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.19e-01

TGCGACTGACGGATCAGGGTGAACCTCAGCCATTACGCCATCAGCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCTTGTGACATGATGTATTCAGCAGCATTGATGAA  >  CP006881/2863386‑2863655
                                                                                                                                      |                                                                                                                                       
tGCGACTGACGGATCAGGGTGAACCTCAGCCATTACGCCATCAGCACCGATCGCTAAAGCGGCTGTAGCGGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTACCATCAACAAAGACTGGCAAATGCGTTTCt                                                                                                                                     >  7:46185/1‑139 (MQ=255)
  cGACTGACGGATCAGGGTGAACCTCAGCCATTACGCCATCAGCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTg                                                                                                                                   >  4:142899/1‑139 (MQ=255)
  cGACTGACGGATCAGGGTGAACCTCAGCCATTACGCCATCAGCACCGATCGCTAAAGCGGCTTTAGCGGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACCTCCACAAAGACTGGCAACTGCGTTTCTTg                                                                                                                                   >  3:74794/1‑139 (MQ=255)
          ggATCAGGGTGAACCTCAGCCATTACGCCATCAGCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCaaaa                                                                                                                           >  5:66872/1‑139 (MQ=255)
           gATCAGGGTGAACCTCAGCCATTACGCCATCAGCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAAt                                                                                                                          <  1:199476/139‑1 (MQ=255)
               aGGGTGAACCTCAGCCATTACGCCATCAGCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCgg                                                                                                                       <  7:152424/138‑1 (MQ=255)
               aGGGTGAACCTCAGCCATTACGCCATCAGCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCgg                                                                                                                       <  4:147105/138‑1 (MQ=255)
                cggTGAAGATCAGCCATTAAGCAATCATCACCGATCGCTAATGAGTCTTTAGCTGTCGGAAGGAAGAGTTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGGAAATGCGTTTCTTGTTTCAAAATCGGca                                                                                                                     <  2:149881/138‑1 (MQ=255)
                   tGACCCTCAGCCACTACGCCAACAGCACCGATCGCCAAAGCGGCTATAGCGGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGAAACATCAACAAAGACTGGCAAAGGCGTTTCTTGTTTCAAAATCGGCACAg                                                                                                                  <  6:206778/139‑1 (MQ=255)
                     aaCCTCAGCCATTACGCCATCAGCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTAGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCt                                                                                                                >  6:296173/1‑139 (MQ=255)
                      aacTCAGCCATTACGCCATCAGCACCGATCGCTAAAGCGTCTTTAGATGTAGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCt                                                                                                                >  2:12297/3‑138 (MQ=255)
                      aCCTCAGCCATTACGCCATCAGCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCt                                                                                                                >  3:43559/1‑138 (MQ=255)
                      aCCTCAGCCATTACGCCATCAGCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGATTCTTGTTTCAAAATCGGCACAGCt                                                                                                                >  7:291492/1‑138 (MQ=255)
                      aCCTCAGCCATTACGCCATCAGCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGATTCTTGTTTCAAAATCGGCACAGCt                                                                                                                >  7:286734/1‑138 (MQ=255)
                           aGCCATTACGCCATCAGCACCGATCGCTAACGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAata                                                                                                          >  5:206778/1‑139 (MQ=255)
                           aGCCATTACGCCATCAGCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAata                                                                                                          >  6:130279/1‑139 (MQ=255)
                                 tACGCCATCAGCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCgg                                                                                                                       <  8:84633/120‑1 (MQ=255)
                                 tACGCCATCAGCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCgg                                                                                                                       >  7:84633/1‑120 (MQ=255)
                                 tACGCCATCAGCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGc                                                                                                    >  1:18023/1‑139 (MQ=255)
                                        tCAGCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTtc                                                                                             >  6:176364/1‑139 (MQ=255)
                                           gCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGATTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCg                                                                                          <  8:286734/139‑1 (MQ=255)
                                           gCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGATTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCg                                                                                          <  8:291492/139‑1 (MQ=255)
                                              ccGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGTATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTg                                                                                       <  5:170206/139‑1 (MQ=255)
                                              ccGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTg                                                                                       >  7:215902/1‑139 (MQ=255)
                                              ccGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTg                                                                                       <  3:305006/139‑1 (MQ=255)
                                                gATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCCCGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGAGTTTCTTTTTTCAAAATCGGCACAGCTGATATATCCAGCGAGTGTCTCGTTGCt                                                                                     >  7:288805/1‑139 (MQ=255)
                                                 aTCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGca                                                                                                                     <  3:12115/106‑1 (MQ=255)
                                                 aTCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGca                                                                                                                     >  4:12115/1‑106 (MQ=255)
                                                 aTCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTg                                                                                    >  2:23020/1‑139 (MQ=255)
                                                 aTCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTg                                                                                    >  8:30394/1‑139 (MQ=255)
                                                 aTCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATAAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTg                                                                                    >  2:147364/1‑139 (MQ=255)
                                                  tCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGt                                                                                   <  7:134888/139‑1 (MQ=255)
                                                        aaGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCata                                                                             <  8:215902/139‑1 (MQ=255)
                                                        aaGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCata                                                                             <  5:176364/139‑1 (MQ=255)
                                                                       tCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATc                                                                                                        >  7:286629/1‑97 (MQ=255)
                                                                       tCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATc                                                                                                        <  8:286629/97‑1 (MQ=255)
                                                                             gCAAGATGGCACGGCTGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATGCCGCGCTcaca                                                        <  1:84994/139‑1 (MQ=255)
                                                                             gCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTcaca                                                        <  7:270016/139‑1 (MQ=255)
                                                                             gCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTcaca                                                        <  5:209650/139‑1 (MQ=255)
                                                                             gCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTcaca                                                        <  3:142899/139‑1 (MQ=255)
                                                                               aaGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGaa                                                                                                             <  5:294143/84‑1 (MQ=255)
                                                                               aaGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGaa                                                                                                             >  6:294143/1‑84 (MQ=255)
                                                                                       aCGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAAt                                               <  7:30394/138‑1 (MQ=255)
                                                                                       aCGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAAt                                               <  3:276606/138‑1 (MQ=255)
                                                                                         ggcggcCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATCCCTTTATTGTTTCAAAATCGGCACAGCTCAAATATCCAGAGTGTTTCTCGTTGCTGTTTCATATGATCTGATTCCGCGCTCACACAGGATAATTTg                                            >  3:111325/1‑139 (MQ=255)
                                                                                                     aaTGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCtt                                 <  6:306639/138‑1 (MQ=255)
                                                                                                       tGCGTAAGATCAAGAAAGACTGGCAATTGCGTTTCTTGTTTCAATATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTGGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCTtgtg                              <  4:111325/139‑1 (MQ=255)
                                                                                                       tGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAg                                                                                                                  >  4:114667/1‑55 (MQ=255)
                                                                                                       tGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAg                                                                                                                  <  3:114667/55‑1 (MQ=255)
                                                                                                       tGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCTtgt                               >  7:25516/1‑138 (MQ=255)
                                                                                                         cGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCTTGTGAc                            <  5:130279/139‑1 (MQ=255)
                                                                                                         cGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCTTGTGAc                            <  7:312894/139‑1 (MQ=255)
                                                                                                         agTACCATCAACAAAGACTGGCACATGCGTTTCTTGTGTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTGTCATATGGTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCTTGTGAc                            <  8:46185/138‑1 (MQ=255)
                                                                                                            aaCATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTg                                                                                       >  5:218314/1‑77 (MQ=255)
                                                                                                            aaCATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTg                                                                                       <  6:218314/77‑1 (MQ=255)
                                                                                                            aaCATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCTTGTGACat                          >  7:168203/1‑138 (MQ=255)
                                                                                                             aCATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCTTGTGACatg                         <  2:18023/138‑1 (MQ=255)
                                                                                                                   acaaAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTcaca                                                        >  2:270721/1‑101 (MQ=255)
                                                                                                                   acaaAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTcaca                                                        <  1:270721/101‑1 (MQ=255)
                                                                                                                   acaaAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCTTGTGACATGATGTAtt                  >  2:85913/1‑139 (MQ=255)
                                                                                                                      aaGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCTTGTGACATGATGTATTcag               <  1:23020/139‑1 (MQ=255)
                                                                                                                      aaGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCTTGTGACATGATGTATTcag               <  1:147364/139‑1 (MQ=255)
                                                                                                                         aCTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTTCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAAtt                                              <  2:48798/105‑1 (MQ=255)
                                                                                                                         aCTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAAtt                                              >  1:48798/1‑105 (MQ=255)
                                                                                                                            ggCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCTTGTGACATGATGTATTCAGCAGCAt         >  7:56719/1‑139 (MQ=255)
                                                                                                                             gCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCTTGTGACATGATGTATTCAgcagca          <  5:296173/137‑1 (MQ=255)
                                                                                                                              cAAATTCGTTTCTTGTTTCAAAATCGGCACAGATGAAATATCCAGCGTGTATCTCGTTGATGTATCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTCGTAATTTCCTTGTGACATGATGTATTCAgcagca          >  4:342316/1‑136 (MQ=255)
                                                                                                                                aaTGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCTTGTGACATGATGTATTCAGCAGCATtgat     <  1:12297/139‑1 (MQ=255)
                                                                                                                                  tGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCTTGTGACATGATGTATTCAGCAGCATtgatg    >  5:233338/1‑138 (MQ=255)
                                                                                                                                   gCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCTTGTGACATGATGTATTCAGCAGCATTGATGaa  >  8:200908/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                       
TGCGACTGACGGATCAGGGTGAACCTCAGCCATTACGCCATCAGCACCGATCGCTAAAGCGGCTTTAGCTGTCGGAAGCAAGAGGTCACGGCGGCCTGTTGAATGCGTAACATCAACAAAGACTGGCAAATGCGTTTCTTGTTTCAAAATCGGCACAGCTGAAATATCCAGCGTGTTTCTCGTTGCTGTTTCATATGTTCTGATTCCGCGCTCACAAAGGATAATTTGGTCATTTCCTTGTGACATGATGTATTCAGCAGCATTGATGAA  >  CP006881/2863386‑2863655

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