Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,489,1740GT41.8% 101.0 / 31.7 67G315C (GGT→TGT) aldAaldehyde dehydrogenase A, NAD‑linked
Reads supporting (aligned to +/- strand):  ref base G (10/29);  new base T (28/0);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 1.15e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

TGAACTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAGGATATTAT  >  NC_000913/1489029‑1489323
                                                                                                                                                 |                                                                                                                                                      
tGAACTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTa                                                                                                                                                     <  2:793613/149‑1 (MQ=255)
 gAACTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAAcc                                                                                                                                                  >  2:417434/1‑151 (MQ=255)
   aCTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCTATTTTGTTAcccc                                                                                                                                                >  2:790618/1‑151 (MQ=255)
      ggCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGc                                                                                                                                              <  2:559873/150‑1 (MQ=255)
      ggCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCt                                                                                                                                             <  2:957068/151‑1 (MQ=255)
         aGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAATTTGGTAAcccc                                                                                                                                                >  2:406507/1‑145 (MQ=255)
         aGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAAcccccccg                                                                                                                                            >  2:883756/1‑145 (MQ=255)
          gTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAAACCCCCTGAAc                                                                                                                                         >  2:74661/1‑151 (MQ=255)
            cAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAGGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACg                                                                                                                                        <  2:596866/150‑1 (MQ=255)
            cAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAAACCCCCCGaa                                                                                                                                          >  1:233675/1‑148 (MQ=255)
            cAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACg                                                                                                                                        <  1:88942/150‑1 (MQ=255)
            cAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGc                                                                                                                                       <  2:901149/151‑1 (MQ=255)
                 ccATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACg                                                                                                                                   <  2:1282677/150‑1 (MQ=255)
                  cATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTAAATTTGGGAAACCCCCTgga                                                                                                                                          >  1:341816/1‑140 (MQ=255)
                      gTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTTACCCCCCCGAACGCAACGACAtt                                                                                                                              >  2:1057467/1‑150 (MQ=255)
                        tGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGc                                                                                                                            <  2:1130438/150‑1 (MQ=255)
                          aTTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAATTTTGTAAACCCGCTGAACGCAACGACAtt                                                                                                                              >  1:170782/1‑146 (MQ=255)
                             cACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATggg                                                                                                                      <  2:135521/151‑1 (MQ=255)
                              aCGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATgggg                                                                                                                     <  1:1262209/151‑1 (MQ=255)
                                    cATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGtt                                                                                                                <  2:1146631/150‑1 (MQ=255)
                                      ttAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAAATTTGGAACCCCGCTGAACGCAACGACATTGCGGTGGGGGCGTTGAt                                                                                                             >  2:403817/1‑151 (MQ=255)
                                         aTAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTATAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAAACCCCCCGGACGCAACGACATTGCGGTGGGGGCGTTGAtt                                                                                                            >  1:382320/1‑149 (MQ=255)
                                           aGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTaa                                                                                                          >  2:1138568/1‑149 (MQ=255)
                                                 cAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAATTTTGTAAACCCCCCGAACGCAACCGCATTTCCAGGGGGGCGGTGGTTaaaacc                                                                                                      >  1:1081269/1‑143 (MQ=255)
                                                   aGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGATTAATTTTTTAACCCCCCTGAAAGCAAAGACACTGCGGTGGGGCCGGTTATTAACGCCgcgg                                                                                                  >  1:113189/1‑149 (MQ=255)
                                                        gtAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTgg                                                                                           <  2:1219704/150‑1 (MQ=255)
                                                          aaCTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAAACCCGCTTAACGCAAAGACATTTCGGTTGGGGCGGTGATTAACGCCCCGGCGCCgg                                                                                           >  1:451232/1‑148 (MQ=255)
                                                          aaCTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGaa                                                                                         <  2:939884/150‑1 (MQ=255)
                                                          aaCTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGATCAATTTTGTCAACCCTCTGAACGACACCACAATACGATGTGGACGTTCATTAAAGCCgcggcg                                                                                                >  1:329986/1‑144 (MQ=255)
                                                           aCTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAg                                                                                       <  1:308239/151‑1 (MQ=255)
                                                                  aGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAATTTTGGAACCCCGCTGAACGCCACGGCATTGGGGAGGGGGCGTTGGTTAACGCCGCGGCGCGGGGAAGGGGCg                                                                                  >  1:44505/1‑149 (MQ=255)
                                                                  aGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCCATTTTGTAACCCCCCTGAACGCAACGACATTGCGATTGGGCCGTTGATTAACGCCGCGGCGCTGGGAAGGGGCg                                                                                  >  2:1044023/1‑149 (MQ=255)
                                                                  aGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAg                                                                                <  2:399815/151‑1 (MQ=255)
                                                                  aGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGATTAAATTTGTAACCCCGCTGAACGCCACGACAATGCGGTGGGGGCGCTTGTTAACACCCCGGCGCTGGAAAGGGGCg                                                                                  >  2:528014/1‑149 (MQ=255)
                                                                   gAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGc                                                                               <  2:564465/151‑1 (MQ=255)
                                                                    aaCGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATTCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACAGTGCGATGGGGCCGTTGAGTAACGCCGCGGCGCTGGAAAGGGTCGAgga                                                                              <  1:732619/151‑3 (MQ=255)
                                                                      cGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAATTTTGTAAACCCGCTGAACGCAACGCCAATTGGATGGGGCCGCTTATTTACGCCCCCGCGCTGGGAAAggggggg                                                                                 >  1:1240335/1‑143 (MQ=255)
                                                                      cGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAAACCCCCCGAACGCAACGACATTGCGATGGGGGCGTTTATTAACGCCGCGGCGCTGGGAAGGGTGGAGCaa                                                                             >  2:1285318/1‑150 (MQ=255)
                                                                      cGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCaaa                                                                            <  2:551779/151‑1 (MQ=255)
                                                                         gtTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAAACCCCCTGGACCCAACGACATTGCGATGGGGGCGTTGGATAACGCCGCGGCGCTGGGGAGGGTCGAggcaaa                                                                           >  2:1055819/1‑144 (MQ=255)
                                                                          tttATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGGCGAGCaaaaa                                                                          >  2:33858/1‑149 (MQ=255)
                                                                           ttATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCGCTTAACGCAACGACATTGCGATGGGGGCGTTGATTAAAGCCCCGGCGCTGGGAAGGGTCGAGCaaaaa                                                                          >  2:556371/1‑148 (MQ=255)
                                                                                  cAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAATTTGGTAACCCCGCTGAAAGCAACGAAATTTCGGTGGGGGCGGTGGTTAACGCCGCGGCGCTGGGAAGGGTCGGGCAAAAAAAGgcgc                                                                   >  1:724629/1‑148 (MQ=255)
                                                                                   aGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCCCTGAACGCAACGACATTTCGATGGGGCCGTTGATTTACCCCGCGGCGCTGGCGAGGGTCGGGCCAAAAAGGgcgcgcgc                                                               >  2:376396/1‑151 (MQ=255)
                                                                                   aGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGGAAACCCGCCGAAAGCAACGACAATGGGATGGGGCCGGTGGTTAAAGCCGCGGGGCGGGGAAGGGTGGAGCAAAAAGGGgggcggg                                                                >  1:910541/1‑148 (MQ=255)
                                                                                    gAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGGTCAATTTTGTAACCCCGCCGAACGCCACGAAAGTGCGATGGGGCCGTTGATTTACGCCGCGGCGCCGGGAAGGGGGGGGTAAAAAGAGgcg                                                                    >  2:567509/1‑145 (MQ=255)
                                                                                     aaaGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAAATTTGTAACCCCCCTGAACCCAACGACATTTCGGTGGGGGCGTTTATTAACACCGCGGCGCTGGGAAAGGGCGAGCaaaaaagtggg                                                                    >  2:1216127/1‑138 (MQ=255)
                                                                                          cATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAATTTTGTAACCCCGCTGCACGCCACGACAATGCGGTTGGGGCGCTGATTAACGCCCCCGCGCTGGAAAGGGGCGAGCCAAAAGAGgggcgcggcggtg                                                          >  2:222127/1‑143 (MQ=255)
                                                                                            tttATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCCATTTGGTGAACCCCCGCAACGCCAAGACCCTTCGAAGGGGGCGTTGGTTTAAACCCCGGCGGTGGGAAAGGGCGAGCCAAAAGGGGGGCGCGCAGTagaag                                                       >  2:1214566/1‑150 (MQ=255)
                                                                                            tttATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAAACCCGCTGAACGCAACGACATTGCGATGGGGGCGGTGATTAACGCCGCGGCGCTGGGAAGGGGCGAGCAAAAAGGGGGGCGGGCCGt                                                            >  2:336745/1‑145 (MQ=255)
                                                                                             ttATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTagaaga                                                      <  1:162347/150‑1 (MQ=255)
                                                                                             ttATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTagaaga                                                      <  1:1195397/150‑1 (MQ=255)
                                                                                                    cAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCga                                              <  2:294953/151‑1 (MQ=255)
                                                                                                            cAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCg                                      <  2:60298/151‑1 (MQ=255)
                                                                                                             aaTCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGGTGATTAATGCCGCGGCGCTGGGAAAGGTCGAGCAAAAAATGGGGCGCGCCGTAGAAGAAGGGGCGAGAGGGGCg                                      >  2:160660/1‑150 (MQ=255)
                                                                                                             aaTCGGCTGGGTGAAGCGATGCAGGCGGTTCAAATTTGTAACCCCCCTGAACGCCACGACACTGCGGTGGGGCCGGTGGTTTAAGCCCCGGGGCTGGGAAGGGGCGAGCCAAAAACGGCGCCGGCAGGAGAAAAAGGGGGGGGAggggggg                                     >  1:746493/1‑148 (MQ=255)
                                                                                                                  gCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTc                                   <  1:408059/148‑1 (MQ=255)
                                                                                                                   cTGGGTGAAGCGATGCAGGCGGTTTAATTTTGTAACCCCGCTGAACACAACCACATTTCGGAGGGGGCGTTGATTAACGCCCCCGCGCTGGAAAGGGTGGGGTAAAAAATTGCGCCCGCAGTAGAAGAAAGGGGGAGGGTggggggggggg                               >  2:407589/1‑142 (MQ=255)
                                                                                                                     gggTGAAGCGATGCAGGCGGTTCCATTTTGTAACCCCGCTGAACGCCACGGCCATGCGGTGGGTGCGGTGATTAACTCCCCCGCGCTGGGAAGGGGGGAGGCAAAAGTGGCGCGCGCCGTAGAAGAAGGGGGGGGAGAgggggttgggggg                             >  2:546456/1‑140 (MQ=255)
                                                                                                                       gTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCaa                           <  1:1057424/151‑1 (MQ=255)
                                                                                                                         gAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCaaa                          <  1:281336/150‑1 (MQ=255)
                                                                                                                            gCGATGCAGGCGGTTTAAATTTGTAACCCCCCTGAAAGCAACGAAAATGCGATGGGGCCGTTGATTAACGCCCCGGGGCTGGGAAAGGGCGAGCAAAAAATGGCGCGCGCAGTAGAAGAAGGGGCGGGAGGGGGGGTTGGGGGGAAAACgg                      >  1:1106639/1‑151 (MQ=255)
                                                                                                                             cGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGt                     <  2:332274/151‑1 (MQ=255)
                                                                                                                                           cAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAGGa        <  1:883716/150‑1 (MQ=255)
                                                                                                                                           cAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACACCGCGGCGCGGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGGGAGAGTGGCGTTCAGTGGCAAAGCGGTAGAGGGGAAAGGat       <  1:1274136/151‑1 (MQ=255)
                                                                                                                                              tttGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTTGCAGAAGGGGCGAGAGTGGCGTCTGGTGGCAAAGCGGGAGAGGGGAAAGGAtatt    >  2:360284/1‑151 (MQ=255)
                                                                                                                                                tGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAGGAtatta   <  2:1277781/150‑1 (MQ=255)
                                                                                                                                                 ggTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAGGAtattat  <  2:255218/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
TGAACTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAGGATATTAT  >  NC_000913/1489029‑1489323

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