Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,656,377 T→C 5.8% L228P (CTG→CCG)  U712_18840 → Putative ABC transporter ATP‑binding protein AlbC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,656,3770TC 5.8% 215.0 / 3.6 69L228P (CTG→CCG) U712_18840Putative ABC transporter ATP‑binding protein AlbC
Reads supporting (aligned to +/- strand):  ref base T (33/32);  new base C (2/2);  total (35/34)
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.38e-01

AGTACTGTGAAGCAGGCGGAACCATTTTGTTTTCTTCTCATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCATTGCCGCAGGCCTTGTCATTCTGGCCGTGT  >  CP006881/3656242‑3656515
                                                                                                                                       |                                                                                                                                          
aGTACTGTGAAGCAGGCGGAACCATTTTGTTTTCTTCTCATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTgcg                                                                                                                                         >  5:366288/1‑139 (MQ=255)
    cTGTGAAGCGGGCGGAACCATTTTGTTTTCTTCTCATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGcggcgg                                                                                                                                     >  3:164653/1‑139 (MQ=255)
       tgAAGCAGGCGGAACCATTTTGTTTTCTTCTCATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGaa                                                                                                                                  >  6:48754/1‑139 (MQ=255)
         aaGCAGGCGGAACCATTTTGTTTTCTTCTCATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAAtt                                                                                                                                >  5:301900/1‑139 (MQ=255)
            cAGGCGGAACCATTTTGTTTTCTTCTCATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGa                                                                                                                                   >  4:65246/1‑133 (MQ=255)
            cAGGCGGAACCATTTTGTTTTCTTCTCATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGa                                                                                                                                   <  3:65246/133‑1 (MQ=255)
               gCGGAACCATTTTGTTTTCTTCTCATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTc                                                                                                                            >  5:89326/1‑137 (MQ=255)
               gCGGAACCATTTTGTTTTCTTCTCATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTc                                                                                                                            >  6:134223/1‑137 (MQ=255)
                           tGTTTTCTTCTCATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGcc                                                                                                               <  1:297501/138‑1 (MQ=255)
                           tGTTTTCTTCTCATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGcc                                                                                                               <  8:329430/138‑1 (MQ=255)
                           tGTTTTCTTCTCATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGcc                                                                                                               <  5:325412/138‑1 (MQ=255)
                           tGTTTTCTTCTCATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGcc                                                                                                               <  1:214091/138‑1 (MQ=255)
                           tGTTTTCTTCTCATCTGCTCGATGGCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGcc                                                                                                               <  4:46978/138‑1 (MQ=255)
                             ttttCTTCTCATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACACACAGATCCAAAAGGTCAATCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCAt                                                                                                             <  4:233842/138‑1 (MQ=255)
                               ttcttcTCATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGa                                                                                                           >  1:193809/1‑138 (MQ=255)
                                tcttctCATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGaa                                                                                                          <  2:23954/138‑1 (MQ=255)
                                tcttctCATCTGCTAGATGTCGTGCTAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGGCATTCCGATTGGAGATGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGaa                                                                                                          <  2:149778/138‑1 (MQ=255)
                                    ctcATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGGCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAAc                                                                                                      <  2:898/138‑1 (MQ=255)
                                       aTCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACAt                                                                                                    <  5:48754/137‑1 (MQ=255)
                                       aTCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATaa                                                                                                  >  1:66385/1‑139 (MQ=255)
                                            cTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATccc                                                                                              >  7:107278/1‑138 (MQ=255)
                                            cTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATccc                                                                                              >  1:97117/1‑138 (MQ=255)
                                            cTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATc                                                                                                <  8:161350/136‑1 (MQ=255)
                                            cTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCATCACAACAAACAGATCCAACAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATc                                                                                                <  5:168923/136‑1 (MQ=255)
                                            cTCGATGTCGTGCAAAGATTTTGTGATTATGAGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATGCCt                                                                                             >  1:348099/1‑139 (MQ=255)
                                              cGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTAt                                                                                           >  5:369090/1‑139 (MQ=255)
                                               gttGTCGTGGATAGATTTTATGATTATGCGGCTATTCTTCAGAACAAACAGATCCAAAAGGTCATTCGGATTGGAGAGGAGAGCTATGTGCGGCGGTAATTTTTCGAGGTGATCGGCCATGTATAACATAATCCCTAt                                                                                           <  3:271936/136‑1 (MQ=255)
                                               gATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTAt                                                                                           <  6:255008/138‑1 (MQ=255)
                                                 tGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCAt                                                                                        >  3:318573/1‑139 (MQ=255)
                                                 tGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCAt                                                                                        >  8:118714/1‑139 (MQ=255)
                                                         aaaGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGc                                                                                <  6:301900/139‑1 (MQ=255)
                                                         aaaGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGc                                                                                <  8:107278/139‑1 (MQ=255)
                                                          aaGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCt                                                                               <  7:118714/139‑1 (MQ=255)
                                                           aGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCCGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTg                                                                              >  3:303309/1‑139 (MQ=255)
                                                           aGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCCGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTg                                                                              >  3:298172/1‑139 (MQ=255)
                                                                   tgATTATGCGGCTATTCTTCACAACAAACAGATACAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAAc                                                                       <  7:37872/138‑1 (MQ=255)
                                                                                   cttcACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGaa                                                                                                          <  1:198430/87‑1 (MQ=255)
                                                                                   cttcACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGaa                                                                                                          >  2:198430/1‑87 (MQ=255)
                                                                                      cacaACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGGATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAggg                                                   >  2:350784/1‑139 (MQ=255)
                                                                                      cacaACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGaa                                                                                                          <  1:22877/84‑1 (MQ=255)
                                                                                      cacaACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGaa                                                                                                          >  2:22877/1‑84 (MQ=255)
                                                                                      cacaACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACAt                                                                                                    <  7:103311/90‑1 (MQ=255)
                                                                                      cacaACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACAt                                                                                                    >  8:103311/1‑90 (MQ=255)
                                                                                      cacaACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATc                                                                                                <  7:326133/94‑1 (MQ=255)
                                                                                      cacaACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATc                                                                                                >  8:326133/1‑94 (MQ=255)
                                                                                           caaacaGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGaa                                                                                                          <  2:13717/79‑1 (MQ=255)
                                                                                           caaacaGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGaa                                                                                                          >  1:13717/1‑79 (MQ=255)
                                                                                                  aTCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGc                                        >  1:323258/1‑138 (MQ=255)
                                                                                                   tCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCaaa                                                                        <  1:145785/105‑1 (MQ=255)
                                                                                                   tCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCaaa                                                                        >  2:145785/1‑105 (MQ=255)
                                                                                                       aaaGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCg                                  >  3:22731/1‑139 (MQ=255)
                                                                                                         aGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCa                                <  6:141168/139‑1 (MQ=255)
                                                                                                         aGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCa                                <  2:323258/139‑1 (MQ=255)
                                                                                                         aGGTCATTCCGATTGGAGAGGAGACCGATCCGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCa                                <  4:303309/139‑1 (MQ=255)
                                                                                                         aGGTCATTCCGATTGGAGAGGAGACCGATCCGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCa                                <  4:298172/139‑1 (MQ=255)
                                                                                                           gTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCAtt                              >  1:19723/1‑139 (MQ=255)
                                                                                                            tCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCCAGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCATTg                             >  6:242297/1‑139 (MQ=255)
                                                                                                            tCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCATTg                             >  2:365412/1‑139 (MQ=255)
                                                                                                            tCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCATTg                             >  2:1425/1‑139 (MQ=255)
                                                                                                              aTTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCATTGcc                           >  6:56158/1‑139 (MQ=255)
                                                                                                               ttCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCATTGCCg                          <  4:22731/139‑1 (MQ=255)
                                                                                                                   gATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCATTGCCGCAgg                      <  6:366288/139‑1 (MQ=255)
                                                                                                                   gATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCATTGCCGCAgg                      <  5:242297/139‑1 (MQ=255)
                                                                                                                     ttGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCATTGCCGCAGGcc                    >  4:355733/1‑139 (MQ=255)
                                                                                                                     ttGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAAACGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCAATGCCGCAGGc                     >  1:364732/1‑138 (MQ=255)
                                                                                                                        gagaggagaCCGATCTGGGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCATTGCCGCAGGCCTTg                 <  5:48620/139‑1 (MQ=255)
                                                                                                                                aCCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACAt                                                                                                    <  5:150376/48‑1 (MQ=255)
                                                                                                                                aCCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACAt                                                                                                    >  6:150376/1‑48 (MQ=255)
                                                                                                                                  cGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCATTGCCGCAGGCCTTGTCATTCTGGc       >  7:175937/1‑139 (MQ=255)
                                                                                                                                       tGCGGCGGGAATTTTTCGAGGTGATCGGCCATCAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCATTGCCGCAGGCCTTGTCATTCTGGCCgtgt  >  4:324346/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                          
AGTACTGTGAAGCAGGCGGAACCATTTTGTTTTCTTCTCATCTGCTCGATGTCGTGCAAAGATTTTGTGATTATGCGGCTATTCTTCACAACAAACAGATCCAAAAGGTCATTCCGATTGGAGAGGAGACCGATCTGCGGCGGGAATTTTTCGAGGTGATCGGCCATGAATAACATAATCCCTATCATGTCTTTGCTGTTCAAACAGCTTTACAGCCGGCAAGGGAAAAAGGACGCCATCCGCATTGCCGCAGGCCTTGTCATTCTGGCCGTGT  >  CP006881/3656242‑3656515

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