Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 2,529,732 T→C 5.5% E122G (GAG→GGG)  U712_12985 ← Uncharacterized protein yrkH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,529,7320TC 5.5% 221.9 / 3.4 73E122G (GAG→GGG) U712_12985Uncharacterized protein yrkH
Reads supporting (aligned to +/- strand):  ref base T (35/34);  new base C (2/2);  total (37/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.09e-01

TGCACCCTTCCCCTTCGCAAAATCTATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGAAAGCATTTCGGATACCATGACGGATGAACCTTCTTTT  >  CP006881/2529602‑2529869
                                                                                                                                  |                                                                                                                                         
tGCACCCTTCCCCTTCGCAAAATCTATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGc                                                                                                                                     >  1:381880/1‑137 (MQ=255)
tGCACCCTTCCCCTTCGCAAAATCTATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCa                                                                                                                                   >  1:232873/1‑139 (MQ=255)
 gCACCCTTCCCCTTCGCAAAATCTATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCAt                                                                                                                                  >  5:162938/1‑139 (MQ=255)
 gCACCCTTCCCCTTCGCAAAATCTATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCAt                                                                                                                                  >  1:237564/1‑139 (MQ=255)
       ttccccttcGCAAAATCTATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTTCCGATTCGAACAAACTGATACCTCCCACCGCCAGAtttt                                                                                                                             >  6:259008/1‑138 (MQ=255)
        tccccttcGCAAAATCTATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTa                                                                                                                            >  8:357496/1‑138 (MQ=255)
         ccccttcGCAAAATATATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAAAAAACTGATACATCTCACCGCCATCTTTTa                                                                                                                            <  5:186247/137‑1 (MQ=255)
             ttcGCAAAATCTATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACAACTCACCGCCAACTTTTAAATccc                                                                                                                      >  1:202257/1‑139 (MQ=255)
               cGCAAAATCTATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATcccc                                                                                                                     >  3:52347/1‑138 (MQ=255)
               cGCAAAATCTATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCg                                                                                                                    >  2:332197/1‑139 (MQ=255)
               cGCAAAATCTATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCg                                                                                                                    >  2:25812/1‑139 (MQ=255)
               cGCAAAATCTATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCg                                                                                                                    >  6:276181/1‑139 (MQ=255)
               cGCAAAATCTATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCACCCTTTACCGATTCGAACAAACTGATACATCTCAACGCCAGCATTTAAATCCCCg                                                                                                                    >  5:228485/1‑139 (MQ=255)
               cGCAAAATCTATAAAGAAATACACCATACGTGTCGCATCAATAATTGCTGCTACACCATAGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATATCACCGCCCACTTTTAAAAcccag                                                                                                                    >  3:312148/1‑137 (MQ=255)
                 cAAAATCTATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTGCCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCg                                                                                                                    >  4:165611/1‑137 (MQ=255)
                 cAAAATCTATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTGCCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCg                                                                                                                    <  3:165611/137‑1 (MQ=255)
                 cAAAATCTATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCCCACCGCCATCTTTTAAATCCCCGAc                                                                                                                  >  6:126072/1‑139 (MQ=255)
                          taAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTAcc                                                                                                           <  4:1595/137‑1 (MQ=255)
                          taAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTAcc                                                                                                           <  6:162938/137‑1 (MQ=255)
                               aaaTTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTa                                                                                                    <  2:232873/139‑1 (MQ=255)
                               aaaTTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTa                                                                                                    <  8:249030/139‑1 (MQ=255)
                                 aTTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTAc                                                                                                            <  8:344198/129‑1 (MQ=255)
                                 aTTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTAc                                                                                                            >  7:344198/1‑129 (MQ=255)
                                     aCCATACGTGTCGCATCAATAATTGCTGGTTGAGCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGTCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGtt                                                                                              <  4:312148/139‑1 (MQ=255)
                                         tACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCt                                                                                          <  1:332197/139‑1 (MQ=255)
                                         tACGTGTCGCATCAATAATTGCTGCTTCACAATTGGAAATGACCATTTATGACAATCAGCCTTTACCGATTCGAACAAACTGATACATCTGACCGACATCTTTTAAATCGCAGACTTTTACATGATCTAAATGtactct                                                                                          >  3:214188/1‑134 (MQ=255)
                                          aCGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTc                                                                                         >  8:235140/1‑139 (MQ=255)
                                            gtgtCGCATCAATAATTGCTGCTTCACCATTGGAAGTGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCa                                                                                       <  7:178217/139‑1 (MQ=255)
                                            gtgtCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCa                                                                                       <  8:165960/139‑1 (MQ=255)
                                               tCGCATCAATAAATGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTAACGGCTCTAAATGTTCGCTCCATGc                                                                                    >  1:18316/1‑139 (MQ=255)
                                                       ataatTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTc                                                                            <  8:260696/139‑1 (MQ=255)
                                                       ataatTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTc                                                                            <  5:161355/139‑1 (MQ=255)
                                                        taatTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTcc                                                                           <  8:314562/139‑1 (MQ=255)
                                                        taatTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTcc                                                                           <  4:52347/139‑1 (MQ=255)
                                                              ctgctTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCt                                                                      <  2:237564/138‑1 (MQ=255)
                                                                gctTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCa                                                                                       <  6:404660/119‑1 (MQ=255)
                                                                gctTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCa                                                                                       >  5:404660/1‑119 (MQ=255)
                                                                 ctTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTa                                                                    <  4:285914/137‑1 (MQ=255)
                                                                 ctTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCCCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTa                                                                    <  5:259008/137‑1 (MQ=255)
                                                                     aCCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTa                                                               <  5:276181/138‑1 (MQ=255)
                                                                       cATTGGACATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCCCCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGa                                                            <  2:154609/139‑1 (MQ=255)
                                                                        aTTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCCCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGAt                                                           <  5:126072/139‑1 (MQ=255)
                                                                           tgAAATGACCATGTATGACAAGCAGCCTTTACCGATTCTAACAAACTGATACATCTCACCGCCATCTATTAAATCCCCGACTTTTACCCGCTCTAAATGTTCCCTCCATCCTTTAATTACACCGAAAAGGTAGAATAct                                                        >  2:388545/2‑138 (MQ=255)
                                                                           ggAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCttt                                                                                 >  5:251690/1‑114 (MQ=255)
                                                                           ggAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCttt                                                                                 <  6:251690/114‑1 (MQ=255)
                                                                              aaTGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCt                                                                                   >  6:371779/1‑109 (MQ=255)
                                                                              aaTGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCt                                                                                   <  5:371779/109‑1 (MQ=255)
                                                                                 gcccATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTa                                                                    >  8:309158/3‑121 (MQ=255)
                                                                                 gACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCTCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGa                                                  >  2:393881/1‑139 (MQ=255)
                                                                                 gACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTa                                                                    <  7:309158/121‑1 (MQ=255)
                                                                                 gACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGa                                                  >  2:327797/1‑139 (MQ=255)
                                                                                 gACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGa                                                  >  4:240637/1‑139 (MQ=255)
                                                                                      tGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGc                                             >  3:16553/1‑139 (MQ=255)
                                                                                      tGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGc                                             >  7:265861/1‑139 (MQ=255)
                                                                                           gACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCtt                                           >  1:163928/1‑136 (MQ=255)
                                                                                            aCAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACATGCTTCAg                                        >  2:193129/1‑138 (MQ=255)
                                                                                              aaGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGaaa                                     <  7:357496/139‑1 (MQ=255)
                                                                                              aaGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGa                                       >  6:191561/1‑137 (MQ=255)
                                                                                               agcagcCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGAAAg                                    <  8:265861/139‑1 (MQ=255)
                                                                                                gcagcCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGAAAg                                    >  1:46770/1‑138 (MQ=255)
                                                                                                   gcCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGAAAGCAtt                                <  1:390941/139‑1 (MQ=255)
                                                                                                    ccTTTACCGATTCGAACAAACTGATACATCTCTCCGCCAGCTTTTAACTCCCCTACTTTTTCCGGCTCTATATGTTCGCTCCATGCTTTCATTCCACCCGTAAGGTAGGATACACAATGACCTGCTTCAGAAAGCAttt                               >  1:62121/1‑139 (MQ=255)
                                                                                                        tACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGAAAGCATTTCgg                            <  3:332523/138‑1 (MQ=255)
                                                                                                         aCCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGTATACAGAATGACCTGCTTCAGAAAGCATTTCGGa                           <  5:191561/138‑1 (MQ=255)
                                                                                                              ttCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGAAAGCATTTCGGATACCAt                     >  6:178107/1‑139 (MQ=255)
                                                                                                              ttCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGAAAGCATTTCGGATACCAt                     <  3:240637/139‑1 (MQ=255)
                                                                                                                    cAAACAGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTAGCGGCTCTAAATGTTCGCTCGATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGAAAGCATTTCGGATACCATGACGGa               <  6:228485/139‑1 (MQ=255)
                                                                                                                     aaaCTGATACATCTCACCGCCATCTTTTATATCCCCGACTTTTACCGGGTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGAAAGCATTTCGGATACCATGACGGAt              <  2:18316/139‑1 (MQ=255)
                                                                                                                     aaaCTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGAAAGCATTTCGGATACCATGACGGAt              <  1:393881/139‑1 (MQ=255)
                                                                                                                     aaaCTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGAAAGCATTTCGGATACCATGACGGAt              <  3:353682/139‑1 (MQ=255)
                                                                                                                            tACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGAAAGCATTTCGGATACCATGACGGATGAACctt       >  8:37068/1‑139 (MQ=255)
                                                                                                                               aTCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTc                                                                            <  7:268542/67‑1 (MQ=255)
                                                                                                                               aTCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTc                                                                            >  8:268542/1‑67 (MQ=255)
                                                                                                                                 ctcACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGaaa                                     >  6:313116/1‑104 (MQ=255)
                                                                                                                                 ctcACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGaaa                                     <  5:313116/104‑1 (MQ=255)
                                                                                                                                 ctcACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGAAAGCATTTCGGATACCATGACGGATGAACCTTCtttt  >  2:365887/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                         
TGCACCCTTCCCCTTCGCAAAATCTATAAAGAAATTCACCATACGTGTCGCATCAATAATTGCTGCTTCACCATTGGAAATGACCATGTATGACAAGCAGCCTTTACCGATTCGAACAAACTGATACATCTCACCGCCATCTTTTAAATCCCCGACTTTTACCGGCTCTAAATGTTCGCTCCATGCTTTCATTCCACCCTTAAGGTAGGATACAGAATGACCTGCTTCAGAAAGCATTTCGGATACCATGACGGATGAACCTTCTTTT  >  CP006881/2529602‑2529869

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