Predicted mutation
evidence seq id position mutation annotation gene description
RA AP012306 1,881,975 A→C intergenic (‑74/+38) yojL ← / ← ompC predicted thiamine biosynthesis lipoprotein/outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,881,9750AC99.7% 1200.0 / ‑5.3 290intergenic (‑74/+38)yojL/ompCpredicted thiamine biosynthesis lipoprotein/outer membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (132/157);  total (132/158)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

Only 100 of 298 total aligned reads displayed.

TGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTG  >  AP012306/1881828‑1882121
                                                                                                                                                   |                                                                                                                                                  
tgtgGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGaaacaa                                                                                                                                                  <  2:3985330/150‑1 (MQ=255)
  tgGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAgg                                                                                                                                                <  2:871368/150‑1 (MQ=255)
      ttttGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGccc                                                                                                                                            <  2:2269016/150‑1 (MQ=255)
      ttttGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGccc                                                                                                                                            <  1:396508/150‑1 (MQ=255)
           aTCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAgg                                                                                                                                       >  2:3028847/1‑150 (MQ=255)
            tCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGc                                                                                                                                      >  1:1043629/1‑150 (MQ=255)
             cGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGcc                                                                                                                                     <  2:1890821/150‑1 (MQ=255)
               caaccaacAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAGGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCt                                                                                                                                   <  2:545700/150‑1 (MQ=255)
               caaccaacAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCACGCAGGCCCt                                                                                                                                   <  1:2265726/150‑1 (MQ=255)
                aaccaacAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCtt                                                                                                                                  >  1:3863457/1‑150 (MQ=255)
                 accaacAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCttt                                                                                                                                 <  2:667868/150‑1 (MQ=255)
                  ccaacAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTg                                                                                                                                <  1:1877519/150‑1 (MQ=255)
                    aacaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGtt                                                                                                                              <  1:1550709/150‑1 (MQ=255)
                     acaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTc                                                                                                                             <  1:2947485/150‑1 (MQ=255)
                     acaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTc                                                                                                                             <  1:4077954/150‑1 (MQ=255)
                       aaaGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGa                                                                                                                           >  1:3334992/1‑150 (MQ=255)
                        aagaagAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGat                                                                                                                          <  1:3252842/150‑1 (MQ=255)
                        aagaagAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGat                                                                                                                          >  2:321509/1‑150 (MQ=255)
                           aagaGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATc                                                                                                                       <  2:3854594/150‑1 (MQ=255)
                            agagCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCa                                                                                                                      <  1:773284/150‑1 (MQ=255)
                            agagCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCa                                                                                                                      >  2:3957383/1‑150 (MQ=255)
                             gagCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCaa                                                                                                                     <  2:1533382/150‑1 (MQ=255)
                             gagCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCaa                                                                                                                     <  1:543920/150‑1 (MQ=255)
                              agCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAAt                                                                                                                    <  2:2498429/150‑1 (MQ=255)
                               gcgcGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAGTc                                                                                                                   <  2:2170142/150‑1 (MQ=255)
                                   ggCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCgaga                                                                                                               >  2:1329111/1‑150 (MQ=255)
                                      agccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTa                                                                                                            <  1:2528760/150‑1 (MQ=255)
                                       gccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAg                                                                                                           <  2:3513300/150‑1 (MQ=255)
                                        ccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTTAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGa                                                                                                          <  2:1386493/150‑1 (MQ=255)
                                         cagcagTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGaa                                                                                                         <  1:2072003/150‑1 (MQ=255)
                                           gcagTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACt                                                                                                       >  1:1491102/1‑150 (MQ=255)
                                           gcagTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACt                                                                                                       <  2:2313333/150‑1 (MQ=255)
                                            cagTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTg                                                                                                      >  1:1792205/1‑150 (MQ=255)
                                               tGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCttt                                                                                                                                 >  2:2727927/1‑120 (MQ=255)
                                                   aCGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAAcc                                                                                               <  2:1911729/150‑1 (MQ=255)
                                                      cgGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGa                                                                                            >  2:712729/1‑150 (MQ=255)
                                                       gggTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGAc                                                                                           <  1:2775293/150‑1 (MQ=255)
                                                           aaaGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCag                                                                                       >  2:1722738/1‑150 (MQ=255)
                                                            aaGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCaga                                                                                      <  2:4061435/150‑1 (MQ=255)
                                                            aaGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCaga                                                                                      >  1:2735263/1‑150 (MQ=255)
                                                              gCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGc                                                                                    <  2:1071856/150‑1 (MQ=255)
                                                              gCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGc                                                                                    >  2:736025/1‑150 (MQ=255)
                                                              gCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGc                                                                                    <  2:3583968/150‑1 (MQ=255)
                                                                   tttCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGa                                                                               >  2:2980206/1‑150 (MQ=255)
                                                                      ccATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGt                                                                            >  1:2137801/1‑150 (MQ=255)
                                                                      ccATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGt                                                                            >  2:282709/1‑150 (MQ=255)
                                                                       cATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGtt                                                                           >  1:829182/1‑150 (MQ=255)
                                                                        aTACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTa                                                                          >  1:3423837/1‑150 (MQ=255)
                                                                        aTACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTa                                                                          <  2:2869993/150‑1 (MQ=255)
                                                                        aTACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTa                                                                          <  2:3668564/150‑1 (MQ=255)
                                                                         tACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTAt                                                                         >  2:934093/1‑150 (MQ=255)
                                                                         tACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTAt                                                                         <  1:3752530/150‑1 (MQ=255)
                                                                           cGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCa                                                                       <  2:1804670/150‑1 (MQ=255)
                                                                             tGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAgt                                                                     >  1:2673531/1‑150 (MQ=255)
                                                                               aTTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAgtgt                                                                   >  1:635044/1‑150 (MQ=255)
                                                                               aTTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAgtgt                                                                   <  1:1725632/150‑1 (MQ=255)
                                                                                ttATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGtt                                                                  >  1:2880704/1‑150 (MQ=255)
                                                                                ttATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGtt                                                                  >  2:2295718/1‑150 (MQ=255)
                                                                                   tCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGAt                                                               <  1:3570447/150‑1 (MQ=255)
                                                                                   tCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGAt                                                               <  1:3517302/150‑1 (MQ=255)
                                                                                     cttaTGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGc                                                             <  2:2577894/147‑1 (MQ=255)
                                                                                     cTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGc                                                             >  2:294684/1‑150 (MQ=255)
                                                                                     cTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGc                                                             <  2:2416631/150‑1 (MQ=255)
                                                                                      tCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGcc                                                            <  1:368596/150‑1 (MQ=255)
                                                                                       cATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCa                                                           <  2:1233573/150‑1 (MQ=255)
                                                                                       cATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCa                                                           <  1:1392701/150‑1 (MQ=255)
                                                                                        aTGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAg                                                          >  1:3113355/1‑150 (MQ=255)
                                                                                         tGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGc                                                         <  2:4110689/150‑1 (MQ=255)
                                                                                          gCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCg                                                        >  2:2376160/1‑150 (MQ=255)
                                                                                          gCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCg                                                        >  1:69715/1‑150 (MQ=255)
                                                                                                ggTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTATAACTGGTAAACCAGACCCAGAGCTACGATGATATCAGTGTTCATGCCAGCGTCACGa                                                  >  2:3440984/1‑150 (MQ=255)
                                                                                                  tCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGt                                                <  2:2461924/150‑1 (MQ=255)
                                                                                                      aaGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAATCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAAc                                            <  2:3456678/150‑1 (MQ=255)
                                                                                                         agTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTgg                                         >  2:1109361/1‑150 (MQ=255)
                                                                                                          gTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGt                                        >  1:1087930/1‑150 (MQ=255)
                                                                                                          gTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGt                                        >  1:601623/1‑150 (MQ=255)
                                                                                                           tACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGtt                                       <  1:871611/150‑1 (MQ=255)
                                                                                                           tACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGtt                                       <  2:1149056/150‑1 (MQ=255)
                                                                                                           tACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGtt                                       <  2:162475/150‑1 (MQ=255)
                                                                                                                   aaaaCTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCca                               <  1:1923699/150‑1 (MQ=255)
                                                                                                                    aaaCTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcag                              <  1:482844/150‑1 (MQ=255)
                                                                                                                    aaaCTGAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcag                              >  1:1260302/1‑150 (MQ=255)
                                                                                                                         gAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGt                         >  2:3495295/1‑150 (MQ=255)
                                                                                                                         gAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGt                         <  1:4143367/150‑1 (MQ=255)
                                                                                                                         gAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGt                         <  1:236124/150‑1 (MQ=255)
                                                                                                                         gAGTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGt                         >  1:3073599/1‑150 (MQ=255)
                                                                                                                          aGTTTGTACGCTGAAAACAATAAAACAAGTTCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTTACGAGTGAAGTGGTTGTCGTCCAGCAggat                        >  2:83935/1‑148 (MQ=255)
                                                                                                                           gTTTGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTg                       >  2:470748/1‑150 (MQ=255)
                                                                                                                             ttGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGAt                     <  2:1794220/150‑1 (MQ=255)
                                                                                                                             ttGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGAt                     >  2:421024/1‑150 (MQ=255)
                                                                                                                              tGTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGAtt                    >  2:1461706/1‑150 (MQ=255)
                                                                                                                               gTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGAttt                   <  2:1039394/150‑1 (MQ=255)
                                                                                                                               gTACGCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGAttt                   <  1:1161688/150‑1 (MQ=255)
                                                                                                                                   gCTGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTa               >  2:592346/1‑150 (MQ=255)
                                                                                                                                     tGAAAACAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGatat                                                                         >  1:2856899/1‑90 (MQ=255)
                                                                                                                                       aaaaCAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCa           <  1:3714977/150‑1 (MQ=255)
                                                                                                                                          aCAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACg        <  1:2908972/150‑1 (MQ=255)
                                                                                                                                          aCAATGAAACAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACg        >  2:1673425/1‑150 (MQ=255)
                                                                                                                                                aaacaaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTg  <  1:3194844/150‑1 (MQ=255)
                                                                                                                                                aaacaaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTg  <  2:3330013/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTG  >  AP012306/1881828‑1882121

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: