Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,073,2020TC43.7% 16.8 / 115.0 87V555A (GTC→GCC) fluCP4‑44 prophage, antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base T (26/23);  new base C (19/19);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 8.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAA  >  NC_000913/2073065‑2073337
                                                                                                                                         |                                                                                                                                       
gAAAGCCGTCAACCCGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCTATTGACCCCACGAATGTc                                                                                                                                        <  1:340023/139‑1 (MQ=255)
 aaaGCCGTCAACCCGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCa                                                                                                                                       >  1:141890/1‑139 (MQ=255)
   aGCCGTCAACCCGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCAct                                                                                                                                     <  1:7017/139‑1 (MQ=255)
       gTCAACCCGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCg                                                                                                                                 >  2:37806/1‑139 (MQ=255)
               gaatgaaGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGt                                                                                                                         <  1:454888/139‑1 (MQ=255)
                aatgaaGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTg                                                                                                                        <  1:412671/139‑1 (MQ=255)
                aatgaaGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTg                                                                                                                        <  2:62263/139‑1 (MQ=255)
                aatgaaGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTg                                                                                                                        <  1:533175/139‑1 (MQ=255)
                aatgaaGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTg                                                                                                                        <  1:550095/139‑1 (MQ=255)
                aatgaaGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTg                                                                                                                        <  1:28800/139‑1 (MQ=255)
                     aGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCAcc                                                                                                                   <  1:37806/139‑1 (MQ=255)
                         acgctgacgctgAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGa                                                                                                               <  2:159911/139‑1 (MQ=255)
                              gacgctgaACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATc                                                                                                          <  2:483273/139‑1 (MQ=255)
                                 gctgaACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCa                                                                                                                                       >  1:495420/1‑107 (MQ=255)
                                 gctgaACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATcccc                                                                                                       >  1:500990/1‑139 (MQ=255)
                                 gctgaACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATcccc                                                                                                       >  1:279563/1‑139 (MQ=255)
                                 gctgaACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATcccc                                                                                                       >  1:482950/1‑139 (MQ=255)
                                     aaCGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATa                                                                                                   >  1:342004/1‑139 (MQ=255)
                                             tACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGt                                                                                                                         >  1:52175/1‑109 (MQ=255)
                                             tACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGt                                                                                                                         <  2:52175/109‑1 (MQ=255)
                                               ccGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATcc                                                                                                         >  1:553287/1‑123 (MQ=255)
                                               ccGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATcc                                                                                                         <  2:553287/123‑1 (MQ=255)
                                                cGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTg                                                                                        <  1:387611/139‑1 (MQ=255)
                                                    accacGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGc                                                                                                                       <  2:367204/104‑1 (MQ=255)
                                                    accacGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGt                                                                                    <  2:132429/139‑1 (MQ=255)
                                                    accacGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACACCACGAATGTCACTCTCGCCTCCGGTGc                                                                                                                       >  1:367204/1‑104 (MQ=255)
                                                       acGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCgg                                                                                 <  2:429434/139‑1 (MQ=255)
                                                         ggATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCggtg                                                                               >  1:524573/1‑139 (MQ=255)
                                                         ggATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCggtg                                                                               >  2:80075/1‑139 (MQ=255)
                                                             gTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTTCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGtggtgg                                                                           >  1:152026/1‑139 (MQ=255)
                                                             gTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGa                                                                                                               >  1:490046/1‑103 (MQ=255)
                                                             gTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGa                                                                                                               <  2:490046/103‑1 (MQ=255)
                                                              tCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGa                                                                          <  2:517986/139‑1 (MQ=255)
                                                              tCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGa                                                                          >  1:517986/1‑139 (MQ=255)
                                                               cATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGAt                                                                         <  2:141890/139‑1 (MQ=255)
                                                                 ttGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGa                                                                       >  1:531259/1‑139 (MQ=255)
                                                                  tGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGAc                                                                      >  2:562229/1‑139 (MQ=255)
                                                                  tGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGAc                                                                      <  1:568365/139‑1 (MQ=255)
                                                                   gCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGaa                                                                                                              >  1:227843/1‑98 (MQ=255)
                                                                   gCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGaa                                                                                                              <  2:227843/98‑1 (MQ=255)
                                                                       aGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAg                                                                 <  1:62077/139‑1 (MQ=255)
                                                                       aGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAg                                                                 <  2:279563/139‑1 (MQ=255)
                                                                         cgcgGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGcc                                                               <  2:152026/139‑1 (MQ=255)
                                                                         cgcgGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGcc                                                               <  2:428954/139‑1 (MQ=255)
                                                                                 aGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGa                                                       <  2:512705/139‑1 (MQ=255)
                                                                                       gAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGAtt                                                 <  1:562229/139‑1 (MQ=255)
                                                                                        aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGa                                                                          >  1:78083/1‑113 (MQ=255)
                                                                                        aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGa                                                                          <  2:78083/113‑1 (MQ=255)
                                                                                        aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACGGATTc                                                >  1:235507/1‑139 (MQ=255)
                                                                                        aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTc                                                >  2:330446/1‑139 (MQ=255)
                                                                                        aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTc                                                >  1:139559/1‑139 (MQ=255)
                                                                                        aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTc                                                >  1:5750/1‑139 (MQ=255)
                                                                                        aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTc                                                >  1:557735/1‑139 (MQ=255)
                                                                                        aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGTCTCCGGTGCCACCTGGa                                                                                                               <  1:350126/76‑1 (MQ=255)
                                                                                        aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGa                                                                                                               >  2:350126/1‑76 (MQ=255)
                                                                                        aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTc                                                >  2:271230/1‑139 (MQ=255)
                                                                                             gACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTc                                           >  2:213977/1‑139 (MQ=255)
                                                                                              aCCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCa                                          >  2:369853/1‑139 (MQ=255)
                                                                                              aCCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCa                                          >  1:237070/1‑139 (MQ=255)
                                                                                                cGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCAcc                                        <  1:330446/139‑1 (MQ=255)
                                                                                                cGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCAcc                                        <  2:285707/139‑1 (MQ=255)
                                                                                                        cTGTGCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCg                                <  2:450194/139‑1 (MQ=255)
                                                                                                          gtgCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGca                              <  1:404752/139‑1 (MQ=255)
                                                                                                          gtgCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGca                              <  2:359238/139‑1 (MQ=255)
                                                                                                           tgCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGcac                             >  1:371021/1‑139 (MQ=255)
                                                                                                           tgCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGa                                                                                                               <  1:50263/57‑1 (MQ=255)
                                                                                                           tgCTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGa                                                                                                               >  2:50263/1‑57 (MQ=255)
                                                                                                             cTGAACGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAg                           >  1:13970/1‑139 (MQ=255)
                                                                                                              tGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAgg                          <  2:342004/139‑1 (MQ=255)
                                                                                                                aaCGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCTCGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGa                        >  2:79789/1‑139 (MQ=255)
                                                                                                                aaCGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGt                                                                                    >  1:2956/1‑79 (MQ=255)
                                                                                                                aaCGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGt                                                                                    <  2:2956/79‑1 (MQ=255)
                                                                                                                aaCGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGAt                                                                         >  2:99739/1‑90 (MQ=255)
                                                                                                                aaCGGTGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGAt                                                                         <  1:99739/90‑1 (MQ=255)
                                                                                                                  cGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAg                      >  2:555713/1‑139 (MQ=255)
                                                                                                                     tGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGcc                                                                                                                      <  1:547946/40‑1 (MQ=255)
                                                                                                                     tGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGcc                                                                                                                      >  2:547946/1‑40 (MQ=255)
                                                                                                                     tGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTc                   >  2:119256/1‑139 (MQ=255)
                                                                                                                     tGCCATTGACCCCACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTc                   >  2:27021/1‑139 (MQ=255)
                                                                                                                               cccACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCaa         >  1:586825/1‑139 (MQ=255)
                                                                                                                               cccACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCaa         >  2:337794/1‑139 (MQ=255)
                                                                                                                               cccACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCAc                                   <  1:334117/113‑1 (MQ=255)
                                                                                                                               cccACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCAc                                   >  2:334117/1‑113 (MQ=255)
                                                                                                                               cccACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCaa         >  1:515356/1‑139 (MQ=255)
                                                                                                                               cccACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCaa         >  1:235955/1‑139 (MQ=255)
                                                                                                                                 cACGAATGCCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAAcc       >  2:470103/1‑139 (MQ=255)
                                                                                                                                      aTGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGaa  <  1:369853/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
GAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAA  >  NC_000913/2073065‑2073337

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