Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi2,663,4660GT23.7% 228.9 / 63.4 97V9L (GTG→TTG) suhBinositol monophosphatase
Reads supporting (aligned to +/- strand):  ref base G (40/34);  new base T (10/13);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTCTTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAAT  >  NC_000913_3_pae_pgi/2663325‑2663613
                                                                                                                                             |                                                                                                                                                   
gTCTTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgca                                                                                                                                             >  1:671654/1‑150 (MQ=255)
 tctttcttCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcag                                                                                                                                            >  1:149353/1‑150 (MQ=255)
 tctttcttCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcag                                                                                                                                            >  4:506607/1‑150 (MQ=255)
         cATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAg                                                                                                                                    <  3:12913/150‑1 (MQ=255)
           acTAAAGATTATTCACGCAGCTTAGCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCGGAACATCGCCTTGCGCGCTGCGCGCACGGc                                                                                                                                  <  2:874621/149‑1 (MQ=255)
            cTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCg                                                                                                                                 >  1:326848/1‑150 (MQ=255)
             tAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCga                                                                                                                                >  1:275116/1‑149 (MQ=255)
             tAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCga                                                                                                                                <  2:725323/150‑2 (MQ=255)
              aaaGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCggg                                                                                                                               >  1:196259/1‑150 (MQ=255)
              aaaGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCggg                                                                                                                               <  3:209758/150‑1 (MQ=255)
              aaaGACTATTCAGGCATCTTATCATACAACGAAGACAGATGCCGAGCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCggg                                                                                                                               <  4:450450/150‑1 (MQ=255)
                aGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCga                                                                                                                                >  2:765938/1‑146 (MQ=255)
                aGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCga                                                                                                                                >  4:320868/1‑146 (MQ=255)
                aGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCga                                                                                                                                <  3:320868/147‑2 (MQ=255)
                aGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCga                                                                                                                                <  1:765938/147‑2 (MQ=255)
                aGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGta                                                                                                                             >  3:77585/1‑150 (MQ=255)
                 gATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaa                                                                                                                            >  1:674331/1‑150 (MQ=255)
                  attattCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCgatgttc                                                                                                                           <  2:275116/150‑7 (MQ=255)
                  attattCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACACCGCCTTGCGCGCTGCGCGCACGGCgatgttc                                                                                                                           <  1:364833/150‑7 (MQ=255)
                   ttattCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCgatgttca                                                                                                                          <  2:326848/150‑8 (MQ=255)
                       tCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaat                                                                                                                      <  2:149353/150‑1 (MQ=255)
                         aCGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTg                                                                                                                    >  3:372383/1‑150 (MQ=255)
                            cATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCggg                                                                                                                               >  2:288191/1‑136 (MQ=255)
                            cATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCggg                                                                                                                               <  1:288191/136‑1 (MQ=255)
                                 tATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatt                                                                                                                          >  1:361825/1‑136 (MQ=255)
                                 tATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatt                                                                                                                          <  2:361825/136‑1 (MQ=255)
                                        aaaCGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGaaa                                                                                                     >  1:525055/1‑150 (MQ=255)
                                              aGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGTCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCgg                                                                                               <  1:686638/150‑1 (MQ=255)
                                                  aGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAACTGAACTGCCAAAAACTATGAAACCCCGGACGc                                                                                           >  3:278027/1‑150 (MQ=255)
                                                    aTGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCga                                                                                                                                >  2:73023/1‑110 (MQ=255)
                                                    aTGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCga                                                                                                                                <  1:73023/111‑2 (MQ=255)
                                                    aTGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcag                                                                                                                                            >  2:198407/1‑99 (MQ=255)
                                                    aTGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcag                                                                                                                                            <  1:198407/99‑1 (MQ=255)
                                                         gATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGaa                                                                                    >  1:494407/1‑150 (MQ=255)
                                                           tctcGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGc                                                                                  <  2:196259/150‑1 (MQ=255)
                                                           tctcGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCACGGACGCTGTGGAAGc                                                                                  <  2:525055/150‑1 (MQ=255)
                                                             tcGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGa                                                                                <  4:413205/150‑1 (MQ=255)
                                                               gctgctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGt                                                                                        <  2:506001/140‑1 (MQ=255)
                                                               gctgctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGt                                                                                        >  1:506001/1‑140 (MQ=255)
                                                               gctgctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGc                                                                              <  1:784267/150‑1 (MQ=255)
                                                               gctgctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGc                                                                              >  4:561465/1‑150 (MQ=255)
                                                                         tctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGc                                                                     >  2:579424/1‑149 (MQ=255)
                                                                         tctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGc                                                                     <  1:579424/149‑1 (MQ=255)
                                                                            gcgcCGTTTTCCCGTTCTTTAACAGCCATTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTa                                                                 <  3:23220/150‑1 (MQ=255)
                                                                              gcCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGa                                                                                                       <  4:66062/110‑1 (MQ=255)
                                                                              gcCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGa                                                                                                       >  3:66062/1‑110 (MQ=255)
                                                                                        cGTTCTTTAACATCCACTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                     >  2:527386/1‑150 (MQ=255)
                                                                                             tttAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAAc                                                <  3:561465/150‑1 (MQ=255)
                                                                                              ttAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAcccc                                                                                                 >  2:221002/1‑100 (MQ=255)
                                                                                              ttAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAcccc                                                                                                 <  1:221002/100‑1 (MQ=255)
                                                                                               tAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGt                                              <  3:428456/150‑1 (MQ=255)
                                                                                                agcATCCAGTGAGAGAGACCGATGCATCCGGTGCTGAACATCGCCGTGCGCGCTGCGCGCAAGGCGGGTAATCTAATTGCCAAAACCTATGAAACCGCGAACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACCATTCCGTGACCAACGTa                                             <  4:132538/148‑1 (MQ=255)
                                                                                                aaCATCCAGGGAGAGAGACCGATTCATCCGATGCTGAACATCGCCTTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACCATGAACCCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTa                                             <  2:812928/150‑1 (MQ=255)
                                                                                                   aTCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCaa                                                                                                                                     >  1:9291/1‑59 (MQ=255)
                                                                                                   aTCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCaa                                                                                                                                     <  2:9291/59‑1 (MQ=255)
                                                                                                   aTCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGAt                                          >  4:179635/1‑150 (MQ=255)
                                                                                                     ccAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATaa                                        >  1:445157/1‑150 (MQ=255)
                                                                                                       aGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTAt                                                                                                         <  1:754594/83‑1 (MQ=255)
                                                                                                       aGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTAt                                                                                                         >  2:754594/1‑83 (MQ=255)
                                                                                                        gTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGc                                     >  4:230268/1‑150 (MQ=255)
                                                                                                         tGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTg                                                                                         >  2:611341/1‑97 (MQ=255)
                                                                                                         tGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTg                                                                                         <  1:611341/97‑1 (MQ=255)
                                                                                                         tGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCCCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATa                                         <  4:55712/145‑1 (MQ=255)
                                                                                                         tGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCCCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATa                                         >  3:55712/1‑145 (MQ=255)
                                                                                                            gagagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATAAAAAAGCAAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGcc                                 <  4:77585/150‑1 (MQ=255)
                                                                                                            gagagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAAATTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAgccgcc                                 >  2:369564/1‑150 (MQ=255)
                                                                                                               agaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGaa                              >  3:201517/1‑150 (MQ=255)
                                                                                                               agaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGaa                              <  3:506607/150‑1 (MQ=255)
                                                                                                                gaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                     <  1:164466/126‑1 (MQ=255)
                                                                                                                gaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                     >  2:164466/1‑126 (MQ=255)
                                                                                                                    gATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGt                         <  2:494407/150‑1 (MQ=255)
                                                                                                                       gCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGAt                      >  4:111911/1‑150 (MQ=255)
                                                                                                                        cATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGAtt                     >  1:994611/1‑150 (MQ=255)
                                                                                                                                  tGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGcc                                 <  2:70727/128‑1 (MQ=255)
                                                                                                                                  tGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGcc                                 >  1:70727/1‑128 (MQ=255)
                                                                                                                                  tGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGa           >  4:310273/1‑150 (MQ=255)
                                                                                                                                  tGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGa           >  2:118913/1‑150 (MQ=255)
                                                                                                                                  aattacccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGa           >  3:345652/8‑150 (MQ=255)
                                                                                                                                   gAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGAt          >  3:559541/1‑150 (MQ=255)
                                                                                                                                   gAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGAt          <  3:111911/150‑1 (MQ=255)
                                                                                                                                      cATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCgg                          >  3:580576/1‑131 (MQ=255)
                                                                                                                                      cATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCgg                          <  4:580576/131‑1 (MQ=255)
                                                                                                                                       cccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCa                                                                    >  2:1025507/3‑88 (MQ=38)
                                                                                                                                       cccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCa                                                                    <  1:1025507/86‑1 (MQ=38)
                                                                                                                                       aTCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGt      >  4:427920/1‑150 (MQ=255)
                                                                                                                                        ccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTg                                                                                         <  2:480076/65‑1 (MQ=25)
                                                                                                                                        ccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTg                                                                                         >  1:480076/2‑66 (MQ=25)
                                                                                                                                        ccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGt                                                       <  2:227652/99‑1 (MQ=39)
                                                                                                                                        ccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGt                                                       >  1:227652/2‑100 (MQ=39)
                                                                                                                                         cGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCgg                          >  2:112798/1‑128 (MQ=255)
                                                                                                                                         cGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCgg                          <  1:112798/128‑1 (MQ=255)
                                                                                                                                           ccGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATc                  >  2:175481/1‑134 (MQ=255)
                                                                                                                                           ccGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATc                  <  1:175481/134‑1 (MQ=255)
                                                                                                                                           ccGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAAt  >  3:356536/1‑150 (MQ=255)
                                                                                                                                           ccGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAAt  <  3:179635/150‑1 (MQ=255)
                                                                                                                                            cGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTg                                                                                         >  2:425471/1‑62 (MQ=255)
                                                                                                                                            cGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTg                                                                                         <  1:425471/62‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GTCTTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAAT  >  NC_000913_3_pae_pgi/2663325‑2663613

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