Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913535,9180AG48.2% 114.4 / 16.2 83E70G (GAG→GGG) hyihydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base A (11/32);  new base G (40/0);  total (51/32)
Fisher's exact test for biased strand distribution p-value = 6.03e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGC  >  NC_000913/535773‑536068
                                                                                                                                                 |                                                                                                                                                      
gagaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGcgcg                                                                                                                                                   <  1:1146094/151‑1 (MQ=255)
   gaaGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCgg                                                                                                                                                  >  1:771183/2‑149 (MQ=255)
        tGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGATTGcc                                                                                                                                           >  2:1037617/1‑151 (MQ=255)
          cgcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGATTGcc                                                                                                                                           >  2:835642/1‑149 (MQ=255)
             aGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATGGGCGGCGGGGGGGCGCGGTATTGCCTTTAt                                                                                                                                      >  1:311767/1‑151 (MQ=255)
              gtgtGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTAtt                                                                                                                                     <  1:872964/151‑1 (MQ=255)
              gtgtGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGAACTGGCGGCGGGGGGGCGCGGGATTGCCCGTAtt                                                                                                                                     >  2:627921/1‑151 (MQ=255)
              gagtgGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTGGGGGGGCGGGGAGGGCGGTATTGACTgggtt                                                                                                                                     >  1:453455/3‑147 (MQ=255)
                gtgGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAGTGGGGGGGGGGGGGGGGCGGGGTTGGCTGTAtt                                                                                                                                     >  1:693065/1‑149 (MQ=255)
                    ttttCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGGGCGGGATTGCCTGTAtttc                                                                                                                                   >  2:582425/1‑145 (MQ=255)
                     tttCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGc                                                                                                                               <  2:263626/150‑1 (MQ=255)
                      ttCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGATGGGGGGGGGGGGGGGGCGGGATTGCCTGtttttccgg                                                                                                                                >  1:834353/1‑143 (MQ=255)
                       tCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGGGGGGGGGCGCGGGATTGCCTGTATTCCTgg                                                                                                                                >  2:226357/1‑147 (MQ=255)
                        cgcgGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTg                                                                                                                           <  1:118431/151‑1 (MQ=255)
                          cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTgaa                                                                                                                         <  1:468959/151‑1 (MQ=255)
                             cgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTgaagaa                                                                                                                      <  1:232540/151‑1 (MQ=255)
                             cgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGGGCGCGGGATTGCCTGTATTTCTgggc                                                                                                                              >  1:721348/1‑141 (MQ=255)
                             cgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCGGGTTTTCCTGTATTTCttg                                                                                                                                >  1:558251/1‑139 (MQ=255)
                              gTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGaagaag                                                                                                                     <  2:220354/151‑1 (MQ=255)
                               ttGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTCCTGGCCGTggag                                                                                                                        >  2:536809/1‑147 (MQ=255)
                                 gAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAagag                                                                                                                   <  2:501714/150‑1 (MQ=255)
                                 gAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGt                                                                                                                  <  1:1136118/151‑1 (MQ=255)
                                  aaTTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGGGCGCGGGATTTCCTGTATTTCTGGCCGTGAAGAagag                                                                                                                   >  1:1150420/1‑149 (MQ=255)
                                         gTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTTCCGGCCGTGGAGAAGAGGTTCggg                                                                                                            >  2:61992/1‑149 (MQ=255)
                                          tttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGGTg                                                                                                         >  2:829074/1‑151 (MQ=255)
                                           ttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATg                                                                                                         <  2:165874/150‑1 (MQ=255)
                                           ttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGCCGGGGAGAAGAGGTTTGGGGTgg                                                                                                        >  1:869879/1‑151 (MQ=255)
                                            tCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTCGGGGGCGGGGGGGCGGGGGAATGACTGTTTTCCTTGCCCGGCAGCAGAATTTCggg                                                                                                            >  2:396915/1‑146 (MQ=255)
                                            tCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGGGGGGGGGGGCGGGATTTCCTGTATTTCTGGGCGTGGAGAAGAGTTTTGGGGTggg                                                                                                       >  1:248887/1‑150 (MQ=255)
                                             ccTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCg                                                                                                      <  2:961253/151‑1 (MQ=255)
                                                tATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTa                                                                                                    <  2:618811/150‑1 (MQ=255)
                                                   gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGGTGGCGGag                                                                                                   >  2:1170904/1‑148 (MQ=255)
                                                    aCTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGGGGGGGGGCGCGGGATTGCCTGTTTTCCTGGCCGTGAAGAAGAGGTTCGGGGTGGCg                                                                                                      >  2:1167718/1‑144 (MQ=255)
                                                       aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagcagca                                                                                              <  2:1219149/149‑1 (MQ=255)
                                                          aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTCGGCGGCGGGGGGGGGCGGTATTTCCTGTAATCCTGGCCGGGGAGAAGAGGTTTGGGGTGGGGGAGGAGCAGCGa                                                                                          >  2:922935/1‑150 (MQ=255)
                                                           cATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGCGGCGGGGGGGCGCGGGATTGCCTGTTTTTCTGGCCGGGAAGAAGAGGTTTGGGGTGGGGTAGCAGCAGCGAtt                                                                                        >  2:315416/1‑151 (MQ=255)
                                                           cATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGAGCGCGGGATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTTGGGGTGGGGTAGCAGCAGCGGtt                                                                                        >  2:606423/1‑151 (MQ=255)
                                                                aagaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGtt                                                                                    <  1:1173137/150‑1 (MQ=255)
                                                                aagaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTa                                                                                   <  1:259278/151‑1 (MQ=255)
                                                                aagaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTa                                                                                   <  2:604886/151‑1 (MQ=255)
                                                                aagaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTCCTGGCCGTGAAGAAGAGGTTTGGGGTGGCGTAGCAGCAGCGCTTCGtt                                                                                    >  2:358110/1‑150 (MQ=255)
                                                                    aTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGcc                                                                               <  1:1250229/151‑1 (MQ=255)
                                                                     ttAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGCGCGGGATTGCCTTTATTTCTGGCCGGGAAGAAGAGGTTCGGGGGGGGGTAGCCGCAGCGATTTGTTAtttcc                                                                              >  1:853564/1‑147 (MQ=255)
                                                                     ttAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGGGCGGGATTTCCTGTATTTCTTGGCGTGAAGAAGAGGTTTGGGGTGGCGTCGCAGCAGCGAGTCCTTTTTccc                                                                              >  2:496750/1‑151 (MQ=255)
                                                                     ttAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACCGGGGGGCGGGGGGGGGCGGGGATTCCTGTATTTCCGGGCGCGTAGAAGGGGTTCAGGGTGGGGTAGCAGGGGCGAATCGTTATAccc                                                                              >  2:1182474/1‑151 (MQ=255)
                                                                      tAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATGGGGGGCGGGGGGGGGCGGGATTTCCCCTTTTCCCGGCCCGGGAGAAGAGGTTTGGGGTGGGGTAGCAGCAGCGATTCGTTTTGccc                                                                              >  1:1214909/1‑150 (MQ=255)
                                                                        aaaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGGGCGGGATTGCCTGTATTTCTGGCCGTGGAGAAGAGGTTTGGGGTGGGGGAGCAGCAGCGGTTTGTTTTGccc                                                                              >  1:1141878/1‑148 (MQ=255)
                                                                          aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGATGGGGGGCGGGGGGGGGCGGGATTGGCTGTATTCCTGGGCGCGGAGAAGAGTTTCGGGAAGGGGGGGCAGGAGCGATTCGttgtagccgt                                                                            >  2:802278/1‑140 (MQ=255)
                                                                          aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGGATTGCCTGTATTCCTGGCCGTGAAGAAGAGGTTTGGGGTGGCGGAGCAGCAGCGATTCGTTATGCCCGTgcg                                                                         >  2:936636/1‑151 (MQ=255)
                                                                               tgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGCGGCGGGGGGGGGCGGGATTGCCTGTATTTCTTGGCGTGGAGAAGAGGTTTGGGGGGGGGGAGCAGCAGCGGTTCGTTATGCCCGGGGGCTggg                                                                    >  1:1153197/1‑151 (MQ=255)
                                                                                 cTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGGGGCGGGGGGGGGCGGGATTGCCTGTATTCCTGGCCGGGGAGAAGAGGTTTGGGGTGGGGGAGCAGCAGCGATTCGTTATGCCCCGGCGCTggg                                                                    >  1:693102/1‑149 (MQ=255)
                                                                                  tGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGATTTCCTGTATTCCTGGCCGGGAAGAAGAGTTTTGGGGTGGGGGAGCAGCAGCGGTTCGTTATGCCCGTGCGCTggg                                                                    >  2:1207236/1‑148 (MQ=255)
                                                                                  tGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTgg                                                                     <  2:1104880/147‑1 (MQ=255)
                                                                                   ggCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaa                                                                 <  1:1045665/150‑1 (MQ=255)
                                                                                    gCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaata                                                               <  1:402432/151‑1 (MQ=255)
                                                                                     cGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaataa                                                              <  1:1081830/151‑1 (MQ=255)
                                                                                       aGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGGGGGCGGGATTGCCTTTATTCCCGGCCCTGGAGGAGAGTTTCGGGGGGGGGTAGCAGCGACGATTCGTTTTGCCCCCGGGCTGGGGaaaaaa                                                             >  1:135384/1‑150 (MQ=255)
                                                                                         taataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGGATTGCCTGGATTTCTGGCCGGGAAAAAGAGGTTTGGGGTGGGGGAGGAGCAACGATTTGTTATGCCCGTGCGCTGGGGAATaaaaa                                                           >  1:1210999/1‑150 (MQ=255)
                                                                                         taataaACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGGGCGGGATTGACCGTTTTTCTTGCCGTGAAGAAGAGGTTTGGGGTGGGGGAGCAGGGGGGGTTTGTTTTTGCCGGGGGCCGGGGAATaaaaa                                                           >  2:950685/1‑150 (MQ=255)
                                                                                          aataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATaaaaaa                                                          <  2:467769/150‑1 (MQ=255)
                                                                                              aaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGGGGGCGGGATTGCCTTTATTCCTGGGCGGGAAGAAGAGATTTGGGGTGGGGGAGCCGGAGGGATTCGGTTTGGCCGTGCGCTGGGGAATaaaaaaa                                                         >  1:129523/1‑147 (MQ=255)
                                                                                              aaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTaa                                                     <  2:902479/151‑1 (MQ=255)
                                                                                                  cGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATCCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACt                                                   <  1:613169/149‑1 (MQ=255)
                                                                                                  cGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGGGGGGGGGGGCGGGATTTGCTGTATTTCTGGCCGGGGAGAAGAGGTTTGGGGTGGCGTAGGAGCAGCGAGTCGTTTTGCCCGTGCGCTGGGGAATAAAAAAATTAAATGt                                                 >  1:1235668/1‑151 (MQ=255)
                                                                                                   gAACATACGCTGCACAATTTACCGGCGGGTGTCTGGGGGGGGGGGGGGGGCGGGATTGCCTGTATTTCTGGGCGTGGGGGAGAGGTTTTGGGGGGGGGAGCAGCAGCGAGTTGGTATGGCCGCGGGCTGGGGAATAAAAAAAATaa                                                     >  1:26116/1‑146 (MQ=255)
                                                                                                   gAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTc                                                <  1:976921/151‑1 (MQ=255)
                                                                                                   gAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGGGGGGGGGGGCGGGATTGCCTGTATTTCTGGCCGTGGAGAAGAGTTTTGGGGTGGGGGAGCAGCAGCGGTTCGTTTTGTCCGTGCGCTGGGGGATAAAAAAATTAAATGt                                                 >  2:976927/1‑150 (MQ=255)
                                                                                                     aCATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGGGGGGGAGCGCGGGATTTCCTTTATTTCTGGGCGGGGAGAAGAGGTTTGGGGGGGGGTAGCAGCAGCGATTTGTTTTGGCCGCGCGCTGGGGAAAAAAAAAAATAACtttat                                               >  1:131584/1‑146 (MQ=255)
                                                                                                      cATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTgg                                             <  1:382510/151‑1 (MQ=255)
                                                                                                       aTACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGATTGCCCGTATTCCCGGCCGGGGAGAAGAGGTTCGGGGTGGCGGAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGGAATAAAAAAAATAAATTTCTgg                                             >  2:285528/1‑150 (MQ=255)
                                                                                                           gCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGGTTGCCTGTATTTCTGGCCGGGGAGAAGAGGTTCGGGGGGGCGGAGCAGGAGGGGTTTGGTTTTCCCCTGCGGTGGGGAATAAAAAAAATAACTGTTTGGGCg                                          >  1:422344/1‑149 (MQ=255)
                                                                                                              gCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTa                                       <  2:136980/149‑1 (MQ=255)
                                                                                                                 caATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAAc                                   <  1:102486/150‑1 (MQ=255)
                                                                                                                   aTTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGcc                                <  1:404521/151‑1 (MQ=255)
                                                                                                                            cggGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGttt                       <  1:248273/151‑1 (MQ=255)
                                                                                                                            cggGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGttt                       <  1:1236622/151‑1 (MQ=255)
                                                                                                                                    tGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTg               >  2:529117/1‑151 (MQ=255)
                                                                                                                                     gggCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTg               <  2:702112/150‑1 (MQ=255)
                                                                                                                                     gggCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTg               >  2:353509/1‑150 (MQ=255)
                                                                                                                                       gcggcgGGGGAGCGCGGGATTGCCTGTATTCCTGGCCGTGAAGAAGAGGTTCGGGGTGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGGAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGGTTCAGCAGTGaa             >  1:499146/1‑150 (MQ=255)
                                                                                                                                               ggAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCAc    >  2:1227097/1‑151 (MQ=255)
                                                                                                                                                gAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACg   >  2:263340/1‑151 (MQ=255)
                                                                                                                                                gAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACg   >  2:316563/1‑151 (MQ=255)
                                                                                                                                                 aGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGc  <  2:826767/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
GAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGC  >  NC_000913/535773‑536068

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