Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB143169,7600AC20.8% 144.9 / 35.9 73intergenic (+724/+2963)PP_RS08750/fhuA_1tRNA‑Val/Ferrichrome outer membrane transporter/phagereceptor
Reads supporting (aligned to +/- strand):  ref base A (30/27);  new base C (6/9);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACGAGAAACTTTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGG  >  ACB143/169654‑169893
                                                                                                          |                                                                                                                                     
aCGAGAAACTTTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCgtg                                                                                                         >  1:21032/1‑137 (MQ=11)
aCGAGAAACTTTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCgtg                                                                                                         >  3:204022/1‑137 (MQ=11)
     gtatcctGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGt                                                                                                    <  1:199575/131‑1 (MQ=11)
     aaaCTTTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGt                                                                                                    <  1:111221/137‑1 (MQ=32)
     aaaCTTTGTATGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACAGATAGTGAACCAGAACCTTGAGGGAAAGGCTAAAAGAACCCCGGAGAGGGCAGTGAAATAGAACATCAAACCGTGTTCGt                                                                                                    >  4:49199/1‑137 (MQ=11)
       aCTTTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTAc                                                                                                  <  8:350207/137‑1 (MQ=32)
       aCTTTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTAc                                                                                                  <  8:267596/137‑1 (MQ=32)
          ttGTCTGAATATGGGGGGACCATCCTCCATGGCTAAATACTACTGACTTACTGTAAGTGAACGAGTTCCGTGAGGGAAAGTAGAATAGAACCCCGGTGAGTGGAGTGTAATAGAACCTGAAACCGTGTGCGTTCAAg                                                                                               <  2:262885/137‑1 (MQ=11)
          ttGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAg                                                                                               <  2:21032/137‑1 (MQ=32)
          ttGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACATGAAACCGTGTGAGTACAAg                                                                                               >  2:57125/1‑137 (MQ=17)
               ttaaTATGGTGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTg                                                                                          <  8:156138/135‑1 (MQ=11)
                 aaTATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTggg                                                                                        <  7:91518/137‑1 (MQ=32)
                         gggACCATCCTCCAAGGCTAAATACTGCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAATGGCGAAAAGAACCCCGGCGTGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGcacgct                                                                                <  1:56596/137‑6 (MQ=11)
                          ggACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGt                                                                                                    >  6:83220/1‑116 (MQ=11)
                          ggACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGt                                                                                                    <  5:83220/116‑1 (MQ=11)
                                tcctccAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAg                                                                         >  3:26852/1‑137 (MQ=32)
                                  ctccAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGtt                                                                           >  6:67934/1‑133 (MQ=32)
                                  ctccAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGtt                                                                           <  5:67934/133‑1 (MQ=32)
                                   tccAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTg                                                                      >  8:318782/1‑137 (MQ=32)
                                       aGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTg                                                                  <  6:340800/137‑1 (MQ=32)
                                        ggCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGc                                                                 >  2:52608/1‑137 (MQ=32)
                                         gCTAAATACTACTGACTGACCGATAGTGAAACAGTACCGTGAGGTAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCg                                                                <  2:48344/137‑1 (MQ=17)
                                           tAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTa                                                              >  3:281422/1‑137 (MQ=32)
                                             aaTACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCgtgtac                                                                                                      >  8:313004/1‑93 (MQ=11)
                                             aaTACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCgtgtac                                                                                                      <  7:313004/95‑3 (MQ=11)
                                              aTACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTgg                                                                                         >  8:261972/1‑107 (MQ=32)
                                              aTACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTgg                                                                                         <  7:261972/107‑1 (MQ=32)
                                                actactGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCttt                                                         >  6:296576/1‑137 (MQ=32)
                                                     tgactgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGtat                                                    >  4:211516/1‑137 (MQ=32)
                                                      gactgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGtata                                                   >  4:322027/1‑137 (MQ=32)
                                                      cactcacCGATAGGGAACCAGTACCGGGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGtata                                                   <  6:95522/136‑1 (MQ=12)
                                                       actgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCAGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATaa                                                  >  2:286345/1‑137 (MQ=17)
                                                         tgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATg                                                >  7:77207/1‑137 (MQ=32)
                                                          gacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTg                                                                  <  2:342206/118‑1 (MQ=32)
                                                          gacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTg                                                                  >  1:342206/1‑118 (MQ=32)
                                                          gacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATgg                                               >  3:199004/1‑137 (MQ=32)
                                                          gacCGATAGTGAACCAGAACCGTCAGGCAAAGGCGAAAAGAACCCCGGAGAGGGGAGAGAAATAGAACCTCAAACCGTGTCCGTACAAGCAGTGGCAGCCTACATTCTTAGGTCACTGCGAACCTTTTGTATAATgg                                               >  8:205935/1‑137 (MQ=11)
                                                            ccGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGt                                             >  1:93293/1‑137 (MQ=32)
                                                             tgttAGTGAACCAGTACCGTGAGTGAAATGCGAAAAGAACCCCGTAGAGGGGAGTGAAATAGATCCTGAAACCGTGTGCGTACAAGGAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                            <  5:79142/134‑1 (MQ=11)
                                                             cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                            <  2:203822/137‑1 (MQ=32)
                                                             cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                            <  7:161940/137‑1 (MQ=32)
                                                             cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                            <  5:129512/137‑1 (MQ=32)
                                                             cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                            <  5:336469/137‑1 (MQ=32)
                                                             cGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                            <  5:376007/137‑1 (MQ=32)
                                                              gATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCa                                           >  6:277448/1‑137 (MQ=32)
                                                               aTAGTGAACCAGTACCGTGAGTGAAAGTCGAAAAGCACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAg                                          <  2:119191/137‑1 (MQ=12)
                                                               aTAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAg                                          >  2:161736/1‑137 (MQ=11)
                                                                ttgTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTTGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGc                                         <  7:287602/135‑1 (MQ=17)
                                                                tAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGc                                         <  4:26852/137‑1 (MQ=32)
                                                                      aCCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTta                                   >  4:13105/1‑137 (MQ=11)
                                                                       ccAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGATAGGGGAGTTAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTTCGTACCTTTTGTATAATGGGTCAGCGACTtat                                  <  3:49199/137‑1 (MQ=12)
                                                                        cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTc                                            <  7:159917/126‑1 (MQ=11)
                                                                        cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTc                                            >  8:159917/1‑126 (MQ=11)
                                                                             ccGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCTTGTGACTGCGTACCTTTTGTATAATGTGTCAGCGACTTATATTctt                            >  8:161906/1‑135 (MQ=11)
                                                                             ccGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTg                                                                                          <  2:216361/75‑1 (MQ=34)
                                                                             ccGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTg                                                                                          >  1:216361/1‑75 (MQ=34)
                                                                             ccGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAg                            >  2:206380/1‑137 (MQ=32)
                                                                             ccGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAg                            >  6:187859/1‑137 (MQ=32)
                                                                              cGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                           >  6:114310/1‑137 (MQ=32)
                                                                                 gAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGc                        <  1:161736/137‑1 (MQ=11)
                                                                                   gggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGa                      >  7:278077/1‑137 (MQ=32)
                                                                                     gAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGc                    >  6:40599/1‑137 (MQ=32)
                                                                                     gAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGc                    >  7:201661/1‑137 (MQ=32)
                                                                                      aaaGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAggt                   <  3:89821/137‑3 (MQ=11)
                                                                                      aaaGGCGAAAAGAACCCCGGAGAGGTGAGTGAAATAGAACCTGAAACCGTGTGCGTACTAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCt                   <  7:110224/137‑1 (MQ=17)
                                                                                      aaaGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCt                   <  1:117462/137‑1 (MQ=32)
                                                                                         ggCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTaa                <  1:244218/137‑1 (MQ=11)
                                                                                           cGAAAAGAACCCCGGAGAGGTGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAAcc              <  5:27686/137‑1 (MQ=17)
                                                                                           cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAAcc              >  6:56771/1‑137 (MQ=32)
                                                                                            gAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCg             <  1:286345/137‑1 (MQ=32)
                                                                                                aGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATa         <  5:40599/137‑1 (MQ=32)
                                                                                                       cGGAGAGGGGAGTGAAATAGAACCTTAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAgg  >  5:306581/1‑137 (MQ=17)
                                                                                                       cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAgg  >  1:343494/1‑137 (MQ=32)
                                                                                                          |                                                                                                                                     
ACGAGAAACTTTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGG  >  ACB143/169654‑169893

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: ‑