Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,4410AC49.3% 87.1 / 32.1 69G384G (GGT→GGGyqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (32/3);  new base C (0/34);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.50e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGA  >  NC_000913/2984292‑2984583
                                                                                                                                                      |                                                                                                                                              
gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCacac                                                                                                                                               <  1:376554/151‑1 (MQ=255)
      gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTgg                                                                                                                                          <  2:1034371/150‑1 (MQ=255)
      gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGc                                                                                                                                         <  2:182830/151‑1 (MQ=255)
        cATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGccc                                                                                                                                       >  2:410739/1‑151 (MQ=255)
        cATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGccc                                                                                                                                       >  1:854452/1‑151 (MQ=255)
                     aaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACaa                                                                                                                           >  2:695330/1‑150 (MQ=255)
                      aaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACaa                                                                                                                           >  1:641721/1‑149 (MQ=255)
                           gggAGCCAAGAAAATTTATGTCGACTTGCCAGCCATCATTATTAAAAAGATTAAGGCCGGAAGATATTTATGGGAAGGGAGAAAATGCAGAGCTATTCGTCACGTTAAAGGGTAAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTg                                                                                                                    <  1:1042774/151‑1 (MQ=255)
                           gggAGCCAAAAAAATTTATGTCGACTCGCCAGCCATTATTATTCAACAGATAAATGCCGGCAGATATTTATCGGTAAGGAGAAAAGGCAGAGTTTTTCGTCACGTTAAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTg                                                                                                                    <  1:913173/151‑1 (MQ=255)
                           gggACCCAAAAAAATGTATGTCGACTTCCCACCAATCTTTTTTAAAAAAATAAATGACGGCAGATATTTACGGGAAAGGAGAAAATGCAGAGTTTTTCGCCACGTTAAAGGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTg                                                                                                                    <  1:938884/151‑1 (MQ=255)
                           gagAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGtt                                                                                                                     >  2:374610/3‑150 (MQ=255)
                            ggAGCCAAGAAAATGTATGTCGACTTGCCAGCCATCATTATTAAACAGATTAATGACGACAGATATTTAGCGGTAAGGAGGAAAGGCAGAGCTATTCGCAACGTTCAATGGTCAAGGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGc                                                                                                                   <  1:326971/151‑1 (MQ=255)
                            gcaccccAAAAAATGTATGCCGCCTGCCCGTCAATCTTTATACAAAAGATTATAGACGGAATTTTTTTATCGGTAAGGAAAAAATGAAGACCTTTTCGCCACGTTAAAGGTTAAATCCCACCCCCCGCCCCCCGCCGATACAAAGCGtt                                                                                                                     <  1:718937/145‑1 (MQ=255)
                             gAGCCAAAAAAATGTATGCCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTAGGGGTAGGGAGAAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGcc                                                                                                                  <  2:886493/151‑1 (MQ=255)
                               gCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGcc                                                                                                                  >  1:103342/1‑149 (MQ=255)
                                  cAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                             >  1:836522/1‑151 (MQ=255)
                                      aattgtatggCGATTTCCCAGCCATCATTATCAAACAGATAAATGACGGCAAATATTTATCGGAAAGGAGAAAATGCGGACTTATTCGCCACTTTCAAGGGCAAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                         <  1:33227/148‑1 (MQ=255)
                                      aaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                         >  2:162112/1‑151 (MQ=255)
                                      aaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                         >  2:250034/1‑151 (MQ=255)
                                       aaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                        >  1:433957/1‑151 (MQ=255)
                                               tCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                 >  2:1105891/1‑150 (MQ=255)
                                                cGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc                                                                                               >  1:994695/1‑151 (MQ=255)
                                                       ccagccaTCTTTATACAAAAAATAAATGAGGGAAAATATTTATGGGAAAGGAGGAAATGCAGAGCTTTTCGTCACTTTCAATGGTAAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                          <  1:606174/149‑1 (MQ=255)
                                                       ccagccaTCATTATTCAACAATTTAATGACGGCAGATATTTACGGATAAGGAGGAAATGCAGAGCTATTCGTCACTTTAAATGGACAATGCCCCACCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                          <  2:967486/149‑1 (MQ=255)
                                                             atcatTATACAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGCCACTTTCAAGGGTCAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCa                                                                                      <  2:597067/147‑1 (MQ=255)
                                                              tcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGATCTTTTCGTCCCGTTCAATGGTCAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt                                                                                    <  2:658891/148‑1 (MQ=255)
                                                               catTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                >  2:646038/1‑151 (MQ=255)
                                                                attattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                               >  1:998376/1‑151 (MQ=255)
                                                                 ttattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTTTTCGTCCCTTTAAATGGTCAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGt                                                                                             <  2:635361/136‑1 (MQ=255)
                                                                   attCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                            >  1:60376/1‑151 (MQ=255)
                                                                       aaaaGATTAATGACGGAAAATTTTTATCGGTAGGGAGAAAAGGCAGAGTTATTCGCCACTTTCAATGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                            <  1:1120799/144‑1 (MQ=255)
                                                                       aaCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                           >  1:247057/1‑148 (MQ=255)
                                                                        aCATATTAATGACGGCAAATTTTTGTCGGAAAGAAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTAAAGGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCagag                                                                       <  2:182281/151‑1 (MQ=255)
                                                                         cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                         >  1:496039/1‑148 (MQ=255)
                                                                          aGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga                                                                      >  1:39202/1‑150 (MQ=255)
                                                                               aaTGACGGCAGATTTTTATGGGTAAGGAGAAAATGCAGGGCTTTTCGTCGCGTTAAATGGTCATTGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGTGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                 <  2:812346/150‑1 (MQ=255)
                                                                               aaTGACGGAAGATATTTATCGGAAAGGAGAAAATGAAGAGCTTTTCGTCACTTCAAAGGGTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                  <  1:783947/149‑1 (MQ=255)
                                                                                  gagggCATATATTTATCGGAAAGGAGAAAAGGAAGAGCTTTTCTTCCCGTAAAAGGGTAAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                               <  2:125053/146‑1 (MQ=255)
                                                                                  gACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt                                                             >  1:377148/1‑151 (MQ=255)
                                                                                   aCGGCAGATATTAAGCGGAAGGGAGAAAAGGAAGGTCTTTTCGCCGCGTTAAATGGTCAAGGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt                                                             <  1:306636/150‑1 (MQ=255)
                                                                                     ggaaaatatTTATCGGAAGGGAGAAAAGGCAGAGCTTTTCGCCACGTTAAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                          <  1:1021414/146‑1 (MQ=255)
                                                                                     ggaaaatatTTAGCGTTAGGGAGAAAATGCAAAGCTTTTCGCCCCTTTAAATGGTAATCGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa                                                              <  1:475713/142‑1 (MQ=255)
                                                                                     ggCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTTTTCGTCACTTTCAATGGTCAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAtt                                                            <  1:326616/149‑1 (MQ=255)
                                                                                           tatTTATGGGTAAGGAGAAAATGCAGAGCTTTTCGTCCCTTTCAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg                                                         <  1:410700/146‑1 (MQ=255)
                                                                                           tatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                    >  2:195442/1‑151 (MQ=255)
                                                                                           tatTTATCGGAAAGGAGAAAATGCAGAGCTATTCGTCACTTTCAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAc                                                     <  2:43493/150‑1 (MQ=255)
                                                                                           tatTTACCGGTAAGGAGAAAAGGCAGAGTTTTTCGTCACTTTAATTGGTCAAGCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt                                                             <  2:945660/142‑1 (MQ=255)
                                                                                             tttttgggTAAGGAGAAAATGCAGAGTTTTTCTTCCGTTTCAATGGTCAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCaa                                                      <  2:869220/141‑1 (MQ=255)
                                                                                                  cGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                             >  1:1054477/1‑151 (MQ=255)
                                                                                                      aaGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGGCAATGCCACACCATCGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                          >  1:1076003/1‑150 (MQ=255)
                                                                                                      aaGGAGAAAATGCAGAGCTATTCGTCACGTTAAATGGACAACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc                                            <  2:322600/148‑1 (MQ=255)
                                                                                                       aggaggAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                          >  1:623623/1‑133 (MQ=255)
                                                                                                       aggaggAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTa                                         >  1:408979/1‑150 (MQ=255)
                                                                                                             aaaGGCAGAGCTATTCGTCACGTTCAAGGGTCAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                  <  1:694333/151‑1 (MQ=255)
                                                                                                              aaggCAGAGTTTTTCGTCGCGTTAATTGGTCACGCCCACCCCCCGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGa                                 <  1:602673/148‑1 (MQ=255)
                                                                                                               aTGCAGACTTTTTCGTCACGTTCAAGGGTAAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAc                                <  1:273060/151‑1 (MQ=255)
                                                                                                                   agagCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                              >  2:678032/1‑149 (MQ=255)
                                                                                                                    gagCTATTCGTCCGGTTAAAGGGTCAACCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt                           <  1:65098/151‑1 (MQ=255)
                                                                                                                        tATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                        >  2:620159/1‑150 (MQ=255)
                                                                                                                         ttttctccacttCAAAGGACAATTGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                        <  1:696495/140‑1 (MQ=255)
                                                                                                                           tCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                    >  2:974177/1‑151 (MQ=255)
                                                                                                                             gTCACTTTAAATGGCAAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                    <  2:775691/149‑1 (MQ=255)
                                                                                                                             gTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgcg                  >  2:655866/1‑151 (MQ=255)
                                                                                                                              tCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                 >  2:757101/1‑151 (MQ=255)
                                                                                                                                 ctttaaagggTAAATGCCCCCCCCGGCCCTCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcca               <  2:2229/142‑1 (MQ=255)
                                                                                                                                 cGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccac              >  2:528238/1‑151 (MQ=255)
                                                                                                                                      aaTGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg         >  1:909665/1‑151 (MQ=255)
                                                                                                                                         ggTCAACGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgac       <  2:982553/150‑1 (MQ=255)
                                                                                                                                            cAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg   >  1:83032/1‑151 (MQ=255)
                                                                                                                                                  ccACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattga  <  1:981656/147‑1 (MQ=255)
                                                                                                                                                   caccCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGTCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg   <  2:840509/145‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGA  >  NC_000913/2984292‑2984583

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