Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,443,2770AC41.9% 114.9 / 13.2 62V274G (GTA→GGA) secYpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (36/0);  new base C (0/26);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 4.77e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAAT  >  NC_000913/3443130‑3443417
                                                                                                                                                     |                                                                                                                                            
gACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTAccc                                                                                                                                           >  2:939626/1‑151 (MQ=255)
   gCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCg                                                                                                                                        >  1:581346/1‑151 (MQ=255)
      tAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGcc                                                                                                                                      >  2:186767/1‑150 (MQ=255)
            aaCACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTc                                                                                                                                >  2:738923/1‑150 (MQ=255)
                  tAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTc                                                                                                                          >  2:692812/1‑150 (MQ=255)
                   aaggCGGTCCCCCGGGCTGAAAATACAGCAAATTTGTTGCCAGCCATTCCCAACCAGTCCCCCCCCCGACCCATGGCGGGTGGGTCGCGGGAAAAAAAAAAATACTGGAAGGAAAGTTT‑GCGGGA‑ATCCCCCCCGCCATATTCACTTTCAg                                                                                                                        <  2:589955/148‑1 (MQ=255)
                     aGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCg                                                                                                                      >  2:684532/1‑151 (MQ=255)
                       cGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGt                                                                                                                    >  2:930203/1‑151 (MQ=255)
                             cccAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgt                                                                                                              >  1:970374/1‑151 (MQ=255)
                             cccAGGCTAAAAAAACAGCAAATTTGTTGCCACCCTTTCCCAACAAGTCCCCCCCCCAAACCAGGGGGCGGTGGCCCCCGGAAAAAAAATAATCCTGAAAGCGAAAATG‑GCCGGG‑ATTCCCCCCGCCATATTCACTTTCAGCGGTAAAt                                                                                                                <  2:780825/149‑1 (MQ=255)
                              ccAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtg                                                                                                             >  1:543190/1‑151 (MQ=255)
                                aGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtgtg                                                                                                           >  1:1070310/1‑151 (MQ=255)
                                    tGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                       >  1:1004990/1‑151 (MQ=255)
                                       aaaTACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                                    >  1:306589/1‑151 (MQ=255)
                                        aaTACAGCAAATTTGTTGTCACCCATTCCCAACCAGTCCCCCCCCCGAACCATGGCGCGATGGTCGCGGGAAACAGAAAAATACGGGAAGCGAAAATG‑GCGGGA‑ATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc                                                                                                   <  2:970338/151‑1 (MQ=255)
                                         ataaagaaaaaTTGTTGTCACCCAGTCCCACCCAGTACCCCCCCCGAACCAGGACGGGATGGTCGCCGGGAACAAAATAATACGGAAAGCGAAGATG‑GCGGGG‑ATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgca                                                                                                  <  1:60902/143‑1 (MQ=255)
                                         aTACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc                                                                                                   >  1:21561/1‑150 (MQ=255)
                                            cAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcagc                                                                                                >  1:856913/1‑150 (MQ=255)
                                              gCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                              >  1:381006/1‑150 (MQ=255)
                                               cGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                              >  2:720166/1‑149 (MQ=255)
                                               cGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                              >  2:779285/1‑149 (MQ=255)
                                                 aaaTTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGac                                                                                          >  1:13438/1‑151 (MQ=255)
                                                 aaaTTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGac                                                                                          >  2:354596/1‑151 (MQ=255)
                                                       tttcccgCCATTCCCACCCACTCCCCCCCCCGAACCAGGCCGCGAGGGTCGCCGGAAAAAAAAAAATACTGGAAGCGAAGATT‑GCGGGG‑ATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgac                                                                                     <  2:142226/144‑1 (MQ=255)
                                                            aGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                               >  1:484457/1‑151 (MQ=255)
                                                              ccattcccACCCAGCCCCCCCCCCGAACCATGCCGCGATGGTCGCCGGAAACAAAATAATACTGGAAGCGAAGTTG‑GCGGGG‑ATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTg                                                                              <  2:1048830/144‑1 (MQ=255)
                                                                        ccagccccccccccAAACCATGACGCGATGGTCGCGGGAAACAAAAAAATACTGGAAGCAAAGATT‑GCCGGA‑ATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg                                                                   <  1:659533/142‑1 (MQ=255)
                                                                        ccAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg                                                                   >  2:939984/1‑151 (MQ=255)
                                                                        ccAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg                                                                   >  2:788168/1‑151 (MQ=255)
                                                                             aCCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                >  1:1113209/1‑149 (MQ=255)
                                                                             aCCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                >  2:352281/1‑149 (MQ=255)
                                                                              ccccccccAAACCAGGGCGCGTTGGTCGCGGGAAACAAAATAATACGGGAAGCGAAGATT‑GGCGGG‑TTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTa                                                               <  1:949789/146‑1 (MQ=255)
                                                                              ccGCACCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGt                                                             >  2:600374/1‑151 (MQ=255)
                                                                                         ccATGGCGCGTTGGTCGCGGGGAAAAAAATAATACTGAAGGCGAAGATG‑GCCGGG‑ATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAAt                                                   <  1:571907/150‑1 (MQ=255)
                                                                                               cgcgGTGGCCGCGGAGAAAAAAAAAATACGGAAAGGGAATTTT‑GCGGGG‑TTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcg                                                <  2:668908/147‑1 (MQ=255)
                                                                                               cgcgAGGGTCGCCGGGAAAAAAAAAATACTGGAAGAGAAGATG‑GCCGGG‑ATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAAt                                                   <  2:945150/144‑1 (MQ=255)
                                                                                                gcgATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGtt                                           >  2:1100845/1‑151 (MQ=255)
                                                                                                  gATGGTCGCGGGAAACAAAATAATCCTGAAAGCAAAGATT‑GCCGGA‑ATGCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcg                                                <  1:790026/144‑1 (MQ=255)
                                                                                                  gATGGCCGCCGGAAACAAAAAAATCCTGGAAGCGAA‑ATT‑GGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcgg                                               <  2:981200/145‑1 (MQ=255)
                                                                                                         gCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGc                                   >  1:236949/1‑150 (MQ=255)
                                                                                                            ggaaaCAAAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTca                                <  1:1018701/147‑1 (MQ=255)
                                                                                                               aaaaaaatactgttGGAGGCGAAGTTG‑GCGGGG‑ATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaacaa                            <  1:54256/138‑1 (MQ=255)
                                                                                                               aaaaaaataatACTGGAAG‑GAAAATTGGCCGGG‑TTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaacaa                            <  2:843801/146‑1 (MQ=255)
                                                                                                               aaCAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaaca                             >  2:642472/1‑150 (MQ=255)
                                                                                                                aaaaaaaaatcCTGGAAGCGAAGATT‑GCGGGG‑ATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaaca                             <  1:306057/142‑1 (MQ=255)
                                                                                                                  aaaataatAGTGGAAGGGAAGATT‑GCGGGG‑ATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                         <  2:571264/149‑1 (MQ=255)
                                                                                                                  aaaaaaataCGGAAAGGGAAGATG‑GCGGGA‑ATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                          <  2:909514/145‑1 (MQ=255)
                                                                                                                   gAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac                        >  2:342075/1‑151 (MQ=255)
                                                                                                                   gAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac                        >  2:1069149/1‑151 (MQ=255)
                                                                                                                    aataatCCTGGATGATAATATT‑GCCGGC‑ATACCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaacaa                            <  2:156355/146‑1 (MQ=255)
                                                                                                                       aataCTGAAAGCGAAGATT‑GCCGGA‑ATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaaca                    <  1:978394/151‑1 (MQ=255)
                                                                                                                            tGGAAGCGAAGATT‑GCCGGG‑ATTAGCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGa                >  2:36489/1‑150 (MQ=255)
                                                                                                                             gaaaGCGAAGTTT‑GCGGGT‑TTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAg                 <  1:866298/146‑1 (MQ=255)
                                                                                                                               aaGCGAAGATT‑GCCGGG‑ATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGaa               <  1:594347/148‑1 (MQ=255)
                                                                                                                                   gAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt        >  2:1007295/1‑151 (MQ=255)
                                                                                                                                   gAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt        >  1:1032042/1‑151 (MQ=255)
                                                                                                                                    aaGATT‑GCCGGG‑ATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt        <  1:834839/150‑1 (MQ=255)
                                                                                                                                     aGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGc      >  2:780233/1‑151 (MQ=255)
                                                                                                                                       aTT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa     >  1:336023/1‑150 (MQ=255)
                                                                                                                                         t‑GCCGGG‑ATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGc      <  1:1002812/147‑1 (MQ=255)
                                                                                                                                           gCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaat  >  1:321463/1‑150 (MQ=255)
                                                                                                                                            gcgggtaTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaat  <  1:1015625/144‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATT‑GCCGGG‑ATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAAT  >  NC_000913/3443130‑3443417

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