Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,188,518 2 bp→CA coding (840‑841/1083 nt) U712_10965 ← Histidinol‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,188,5180TC97.6% 264.5 / NA 83T281A (ACA→GCA) U712_10965Histidinol‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (38/43);  minor base A (0/2);  total (38/45)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*CP0068812,188,5190TA97.6% 254.3 / NA 82K280N (AAA→AATU712_10965Histidinol‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base A (38/42);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTTTTTCCAGCAGAGCCTGAAACAGTTCATCGGAAGGGCGGTTAAAATCAATTAACACAAAGTTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTCTTCGCAAAATCATAATATTGCTGAAGGCCTGTTTTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCGATCGCCGCAGCCTGACCTAGACGGCTCGTGTTAAATGGCTCTCTCGCCGGCTCAATCTGGCGAATCAGGTTCTCGTCAG  >  CP006881/2188387‑2188649
                                                                                                                                   ||                                                                                                                                  
gCCTTTTTCCAGCAGAGCCTGAAACAGTTCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCattatt                                                                                                                               <  8:92592/138‑1 (MQ=255)
 ccTTTTTCCAGCAGAGCCTGAAACAGTTCATCAGCAGGAAGTTTAAAGTCAATTACCACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGACTGCATTAtttt                                                                                                                             >  3:429241/1‑139 (MQ=255)
     tttCCAGCAGAGCCTGAAACAGTTCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCt                                                                                                                          <  4:177280/138‑1 (MQ=255)
      ttCCAGCAGAGCCTGAAACAGTTCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTc                                                                                                                         >  4:172552/1‑138 (MQ=255)
      ttCCAGCAGAGCCTGAAACAGTTCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTc                                                                                                                         >  7:137706/1‑138 (MQ=255)
      ttCCAGCAGAGCCTGAAACAGTTCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTc                                                                                                                         >  5:383547/1‑138 (MQ=255)
      ttCCAGCAGAGCCTGAAACAGTTCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTc                                                                                                                         >  1:58681/1‑138 (MQ=255)
         cagcagAGCCTGAAACAGTTCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGaca                                                                                                                     <  3:172552/139‑1 (MQ=255)
         cagcagAGCCTGAAACAGTTCATCAGCAGGACGTTTAAAGTCAAATAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGaca                                                                                                                     <  2:292025/139‑1 (MQ=255)
           gcagAGCCTGAAACAGTTCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGaca                                                                                                                     >  2:232433/1‑137 (MQ=255)
             agagCCTGAAACAGTTCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAgg                                                                                                                 <  7:173281/139‑1 (MQ=255)
             agagCCTGAAACAGTTCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAgg                                                                                                                 <  2:317393/139‑1 (MQ=255)
              gatccTGAAACAGTTCATCAGCAGGACGTTTAAAGTCAATTATCACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGa                                                                                                                <  8:328767/136‑1 (MQ=255)
                gCCTGAAACAGTTCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAg                                                                                                              <  2:162368/139‑1 (MQ=255)
                gCCTGAAACAGTTCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAg                                                                                                              <  4:204401/139‑1 (MQ=255)
                 ccTGAAACAGTTCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCCCTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAg                                                                                                              >  5:386856/1‑138 (MQ=255)
                 ccTGAAACAGTTCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAg                                                                                                              >  6:98630/1‑138 (MQ=255)
                           ttCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCt                                                                                                                          >  4:258324/1‑116 (MQ=255)
                           ttCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCt                                                                                                                          <  3:258324/116‑1 (MQ=255)
                           ttCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACATGAAGCAATAAACGCt                                                                                                   <  1:396504/139‑1 (MQ=255)
                           ttCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCt                                                                                                   >  8:264992/1‑139 (MQ=255)
                           ttCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCt                                                                                                   >  8:152736/1‑139 (MQ=255)
                           ttCATCAGCAGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCt                                                                                                   >  4:120030/1‑139 (MQ=255)
                              atcagcaGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGa                                                                                                <  6:210752/133‑1 (MQ=255)
                                cagcaGGACGTTTAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGat                                                                                               >  5:350304/5‑138 (MQ=255)
                                     ggacgtttAAAGTCAATTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGatcatc                                                                                           >  8:52298/7‑137 (MQ=255)
                                                   aaTTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTACACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATcgcc                                                                           >  4:161427/1‑139 (MQ=255)
                                                   aaTTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATcgcc                                                                           >  7:287361/1‑139 (MQ=255)
                                                    aTTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATcgccg                                                                          <  4:359996/139‑1 (MQ=255)
                                                    aTTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATcgccg                                                                          <  6:350304/139‑1 (MQ=255)
                                                    aTTAACACAAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATcgccg                                                                          >  8:121487/1‑139 (MQ=255)
                                                        acacAAAATTTGTCTGTGAAGGGTAGCACTGTAAACAATGAGTGTTCTCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTTGACACAGGAAGCAATAAACGCTTGATCATCACGTGCTGCAATCGCCGCatc                                                                      <  8:345203/139‑3 (MQ=255)
                                                        acacAAAATTTGTCTGTGAAGGGTAGCACATTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGc                                                                      <  6:386856/139‑1 (MQ=255)
                                                         cacaAAATTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGcc                                                                     <  4:218567/139‑1 (MQ=255)
                                                              aatttGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGcc                                                                <  2:66890/137‑1 (MQ=255)
                                                              aatttGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGcc                                                                <  4:393931/137‑1 (MQ=255)
                                                              aatttGTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGcc                                                                <  4:57921/137‑1 (MQ=255)
                                                                   gTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATcgccg                                                                          <  1:358752/124‑1 (MQ=255)
                                                                   gTCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATcgccg                                                                          >  2:358752/1‑124 (MQ=255)
                                                                    tCTGTGAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACg                                                          >  3:251263/1‑139 (MQ=255)
                                                                      ggtgAAGGGGAGCACTTTAAACCATTAGTTTTCGCAAAATCATAAAACGGCGGAAGGCCTGCATAATTTTGCGCGACACAGGAAGCAAAAACCGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGc                                                        <  2:178434/138‑1 (MQ=255)
                                                                        tgAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAggcggc                                                        <  8:154664/137‑1 (MQ=255)
                                                                        tgAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAg                                                             <  2:177247/132‑1 (MQ=255)
                                                                        tgAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAg                                                             >  1:177247/1‑132 (MQ=255)
                                                                        tgAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGc                                                        <  7:422397/137‑1 (MQ=255)
                                                                        tgAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGc                                                        <  5:193672/137‑1 (MQ=255)
                                                                        tgAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGc                                                        <  2:58681/137‑1 (MQ=255)
                                                                        tgAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGc                                                        <  7:264992/137‑1 (MQ=255)
                                                                        tgAAGGGTAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCtg                                                      >  7:167841/1‑138 (MQ=255)
                                                                        tgAAGGGTAGCACTTTAAACCATGAGTATTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGc                                                        <  6:165051/137‑1 (MQ=255)
                                                                        tgAAGGGTAGCAATTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGAAATCGCCGCAGCCTGGCCCAGACGGc                                                        <  5:249848/137‑1 (MQ=255)
                                                                          aaGGGGAGCACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCtgg                                                     <  7:152736/138‑3 (MQ=255)
                                                                           aGGGTAGCACTTTAAACCTTGAGTTTTCGAAAAATCATAATCCTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCTGgta                                                   <  4:307515/139‑2 (MQ=255)
                                                                                  cACTTTAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCTGGTATTaaacg                                            >  3:307515/1‑137 (MQ=255)
                                                                                       tAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGTTGCAATCGCCGCAGCCTGGCCCAGACGGCTGGTATTAAACGGctct                                       >  2:147354/1‑139 (MQ=255)
                                                                                       tAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCt                                                       <  3:313956/123‑1 (MQ=255)
                                                                                       tAAACCATGAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCt                                                       >  4:313956/1‑123 (MQ=255)
                                                                                            cATGAGTTTTCGCAAAATCAAAATACTGCTGAAGGCCGGCATAATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCAACAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCTGGTATTAAACGGCTCTCTCGc                                  <  4:325462/139‑1 (MQ=255)
                                                                                               gAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCt                                                                                                                          <  8:138266/48‑1 (MQ=255)
                                                                                               gAGTTTTCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCTGGTATTAAACGGCTCTCTCGCCgg                               >  4:276871/1‑139 (MQ=255)
                                                                                                    ttCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGc                                                        >  5:374901/1‑109 (MQ=255)
                                                                                                    ttCGCAAAATCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCTGGTATTAAACGGCTCTCTCGCCGGCTCaa                          >  3:54180/1‑139 (MQ=255)
                                                                                                    ttCGCAAAATCATAATACTGCTGAAGGCCTGAATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGc                                                        <  6:374901/109‑1 (MQ=255)
                                                                                                         aaaaTCATAATACTGCTGAAGGCCTGCATTATTTTGCGCGACCCAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCGGGCCCAGACGGCTGGTATTAAACGGCTCTCTCGCCGGCTCAATCTgg                     <  4:172778/139‑1 (MQ=255)
                                                                                                          aaaTCATAATACTGCTGATGGCCTGCATTATTTTGCTCGACAACGGAAGCAATAAACGCTTGATAATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCTGGTATTAAACGGCTCTCTCGCCGGCTCAATCTGGc                    <  8:167841/139‑1 (MQ=255)
                                                                                                          aaaTCATAATACTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGTTGCAATCGCCGCAGCCTGGCCCAGACGGCTGGTATTAAACGGCTCTCTCGCCGGCTCAATCTGGc                    <  1:147354/139‑1 (MQ=255)
                                                                                                                  atacTGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCTGGTATTAAACGGCTCTCTCGCCGGCTCAATCTGGCGAATCAgg            <  8:137706/139‑1 (MQ=255)
                                                                                                                     ctGCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCTGGTATTAAACGGCTCTCTCGCCGGCTCAATCTGGCGAATCAGGTtc         >  4:399570/2‑139 (MQ=255)
                                                                                                                       gCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCTGGTATTAAACGGCTCTCTCGCCGGCTCAATCTGGCGAATCAGGTtctc       >  1:110133/1‑139 (MQ=255)
                                                                                                                       gCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCTGGTATTAAACGGCTCTCTCGCCGGCTCAATCTGGCGAATCAGGTtctc       >  5:446827/1‑139 (MQ=255)
                                                                                                                       gCTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCTGGTATTAAACGGCTCTCTCGCCGGCTCAATCTGGCGAATCAGGTtctc       >  5:191828/1‑139 (MQ=255)
                                                                                                                        cTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGc                                                        <  1:380511/89‑1 (MQ=255)
                                                                                                                        cTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGc                                                        >  2:380511/1‑89 (MQ=255)
                                                                                                                        cTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCTGGTATTAAACGGCTCTCTCGCCGGCTCAATCTGGCGAATCGGGTTCTCg      >  7:231121/1‑139 (MQ=255)
                                                                                                                        cTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCTGGTATTAAACGGCTCTCTCGCCGGCTCAATCTGGCGAATCAGGTtctc       <  6:383547/138‑1 (MQ=255)
                                                                                                                        cTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCTGGTATTAAACGGCTCTCTCGCCGGCTCAATCTGGCGAATCAGGTTCTCg      >  2:400204/1‑139 (MQ=255)
                                                                                                                        cTGAAGGCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCTGGTATTAAACGGCTCTCTCGCCGGCTCAATCTGGCGAATCAGGTTCTCg      >  7:50672/1‑139 (MQ=255)
                                                                                                                        cTGAAGGCCTGCATTATTTTGATCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCACCGACGGCTGGTATTAAACGGCTCTCTCGCCGCCTCAATCTGGCGAATCAGGTtctct      >  7:373836/1‑138 (MQ=255)
                                                                                                                             ggCGTGCTTTTTTTTGCTCGAAATATTAATAAATATAGGCTTTATCAGAAAGTGTTGTATTGGCCGCAGCCTGGCCCAGAGGGGTGGTATTAAACGGCTCTCTCGCCGGCTCAATCTGGCGAATCAGGTTCTCGTCAg  <  8:373836/138‑1 (MQ=255)
                                                                                                                             ggCCTTATTTATTTTGCTCGACACTGTAAGCAAAAAAGGCTTGATCATCATGTGCTCAAAGCGCCGCAGCCTGGCCCAGACGGCTGGTATTAAACGGCTCTCTCGCCGGCTCAATCTGGCGAATCAGGTTCTCGTCAg  <  2:225105/138‑1 (MQ=255)
                                                                                                                             ggCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGtgctgcaa                                                                                <  2:294794/60‑3 (MQ=255)
                                                                                                                             ggCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGtgctgcaa                                                                                >  1:294794/1‑58 (MQ=255)
                                                                                                                             ggCCTGCATTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCAATCGCCGCAGCCTGGCCCAGACGGCTGGTATTAAACGGCTCTCTCGCCGGCTCAATCTGGCGAATCAGGTTCTCGTCAg  >  3:414548/1‑138 (MQ=255)
                                                                                                                                   ||                                                                                                                                  
GCCTTTTTCCAGCAGAGCCTGAAACAGTTCATCGGAAGGGCGGTTAAAATCAATTAACACAAAGTTTGTCTGTGAAGGGTAGCACTTTAAACCATGAGTCTTCGCAAAATCATAATATTGCTGAAGGCCTGTTTTATTTTGCTCGACACAGGAAGCAATAAACGCTTGATCATCAAGTGCTGCGATCGCCGCAGCCTGACCTAGACGGCTCGTGTTAAATGGCTCTCTCGCCGGCTCAATCTGGCGAATCAGGTTCTCGTCAG  >  CP006881/2188387‑2188649

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: