Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,842,9140AG42.7% 102.1 / 30.7 77M182V (ATG→GTG) gdhAglutamate dehydrogenase, NADP‑specific
Reads supporting (aligned to +/- strand):  ref base A (4/39);  new base G (31/1);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.35e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGTTTT  >  NC_000913/1842765‑1843060
                                                                                                                                                     |                                                                                                                                                      
aTTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGgat                                                                                                                                                       <  2:920498/151‑1 (MQ=255)
aTTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGgat                                                                                                                                                       <  2:915467/151‑1 (MQ=255)
aTTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGgat                                                                                                                                                       <  2:420296/151‑1 (MQ=255)
aTTTCGATCCGAAAGAAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCCCCTGGGCGGGGATACCGACGTCCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTAGGGCGGGgat                                                                                                                                                       <  2:559931/151‑1 (MQ=255)
   tCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGgatga                                                                                                                                                     <  2:1173288/150‑1 (MQ=255)
     gATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGgggatgg                                                                                                                                                  >  2:256009/1‑150 (MQ=255)
     gATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAtgatga                                                                                                                                                  <  2:262490/151‑1 (MQ=255)
        ncGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGNCGNNCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGATGgaa                                                                                                                                                >  2:110765/2‑150 (MQ=255)
           aaaGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGATGaaaaaa                                                                                                                                             >  1:868788/1‑149 (MQ=255)
             aGGAAAAAGCGAAGGTGAAGTGATGCGTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGATGgaaaaaatct                                                                                                                                         >  2:1159158/1‑146 (MQ=255)
               gAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAAGACAAACCcctc                                                                                                                                        >  1:915562/1‑151 (MQ=255)
                aaaaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTcc                                                                                                                                       <  1:289389/151‑1 (MQ=255)
                 aaaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCca                                                                                                                                      <  1:762654/151‑1 (MQ=255)
                 aaaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCca                                                                                                                                      <  1:135576/151‑1 (MQ=255)
                   aaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGGACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATACCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCcaac                                                                                                                                    <  2:385595/151‑1 (MQ=255)
                   aaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCcaac                                                                                                                                    <  2:264130/151‑1 (MQ=255)
                    aGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCcaaca                                                                                                                                   <  1:749406/151‑1 (MQ=255)
                         aGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATAcc                                                                                                                             <  2:1451456/151‑1 (MQ=255)
                             gAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGTGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGGGGGAAAAAATCTTCAACAATaaccccc                                                                                                                         >  1:1344869/1‑145 (MQ=255)
                             gAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGAAAAAAAGCTCCCACAAAACccccc                                                                                                                         >  2:770465/1‑150 (MQ=255)
                              aaGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGAGGAAAAGATATCCAAAAAAAACGcccg                                                                                                                        >  1:1435507/1‑149 (MQ=255)
                              aaGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGAGGAAAAAAATCTCCCACAAAACCGcccg                                                                                                                        >  1:1342075/1‑149 (MQ=255)
                                   aTGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGGGAAAAAAATCTCCCACAAAACCGCCCGCGTCt                                                                                                                  >  2:320528/1‑151 (MQ=255)
                                      cGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGATGAAAAAACCCTCCCACCCAAACGCCTGCGTCttt                                                                                                                >  1:482281/1‑149 (MQ=255)
                                      cGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCa                                                                                                              <  1:823173/151‑1 (MQ=255)
                                      cGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGGCGTGAAGTCGGCTTTATGGCGGGGCGGTAGAAAAAGTTGCCCCAAAATACCggcgccaccttat                                                                                                                >  1:409654/1‑139 (MQ=255)
                                         ttttGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTTTTGAAAAAACTCTCCCACAATACCCCCCCCGTTTTccccc                                                                                                           >  1:971467/1‑150 (MQ=255)
                                         ttttGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCAcc                                                                                                            <  1:808542/150‑1 (MQ=255)
                                                aGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCAGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAggg                                                                                                   <  2:1291618/151‑1 (MQ=255)
                                                    gcTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCtt                                                                                               <  1:368398/151‑1 (MQ=255)
                                                       gatgaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTc                                                                                             <  1:1047482/150‑1 (MQ=255)
                                                       gatgaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCa                                                                                            <  1:573577/151‑1 (MQ=255)
                                                         tgaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGAAAAAACTCTCCCACAAAAACCCCCGCGGCTTTACCCGGAAGGGGCTTTTAtt                                                                                          >  1:4809/1‑151 (MQ=255)
                                                          gaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGGGTAAAAACTCTCCCACAAAACCGCCCGCGTCTTCCCCGGGAAGGGGCTTTCCttt                                                                                         >  1:546280/1‑151 (MQ=255)
                                                          gaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGAAAAAACTCTCCCACAAAAACCCCTGCGTCTTCCCCCGTAAGGGGCCTTCtttt                                                                                         >  2:1020684/1‑151 (MQ=255)
                                                          gaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCAttt                                                                                         <  1:1455927/151‑1 (MQ=255)
                                                              gAACTGTATCGCCACCTGGGCGCGGATCCCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTggc                                                                                      <  1:290535/150‑1 (MQ=255)
                                                                     aTCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGAGGAAAAAAATCTCCAACAATACCGCCCCCGGCTTTTCCCGGAAAGGGCCTTTTTTTTGGGGGAGTCt                                                                              >  1:170839/1‑151 (MQ=255)
                                                                       cGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGAAAAAAATCTCCCACAAAAACGCCCGCGTCTTCCCCCGGGAGGGCCCTTTCtttttgggggggctt                                                                              >  2:806681/1‑139 (MQ=255)
                                                                       cGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCtt                                                                             <  1:225333/150‑1 (MQ=255)
                                                                        gCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGAAAAAAGTCCCCCACAAAACCCCCCCCGGCTTTACCCGGAAGGGGCCTTTATTTTGCGGCAGTCCTAt                                                                           >  1:859556/1‑151 (MQ=255)
                                                                           aCCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGGAAAAACTCTCCCACAAAAACCCCCGCGTCTTCACCCGTAAGGGGCTTTCTTTTTGCGGCAGGCTtttttt                                                                        >  2:402838/1‑149 (MQ=255)
                                                                           aCCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGGAAAAAAACTCCCACAAAAACGCCCCCGGCTTCACCCGGGAGGGCCCTTTATTTTGCGGCAGGCTtttttc                                                                        >  2:61067/1‑149 (MQ=255)
                                                                            ccTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGAAAAAAATCTCCAACAAAAACCCCCCCGTCTTCCCCCGGAAAGGGCTTTCTTTTTGCGGCCGGCTttttttc                                                                       >  2:1021075/1‑148 (MQ=255)
                                                                               gggCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCg                                                                    <  1:1076331/151‑1 (MQ=255)
                                                                                 gcgcGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGa                                                                  <  2:240480/151‑1 (MQ=255)
                                                                                  cgcgGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGaa                                                                 <  2:279322/151‑1 (MQ=255)
                                                                                   gcgGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGAAAAAAATCTCCCACAAAAACCCCCGCGTCTTCACCCGGAAGGGGCTTTTATTTTGCGGCCGGCTTTTTTGccccgg                                                                  >  2:198217/1‑146 (MQ=255)
                                                                                   gcgGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGCGGGGATGAGAAAAACCTTCAACAATAACCGCCCGGGCTTCCACCGTGAGGGGCCCTCCTTTTGGGGGCGACTTAtttgtcccgg                                                                  >  2:199867/1‑141 (MQ=255)
                                                                                     ggATACCGACGTTCCGGCAGGTGATGTCGGGGTTGGTGGTCGTGAAGTCGGCTTTAGGGCGGGGGTGATGAAAAGGCTCTCCAACAATGCCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCt                                                              <  1:1454794/151‑1 (MQ=255)
                                                                                      gATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTa                                                             <  1:709409/151‑1 (MQ=255)
                                                                                          ccGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGAAAAAACTCTCCCACAATAACGCCCGCGTCTTTCCCCGGAAGGGGCTTTCATTTTGGGGGAGGCTTTTTTGCCCCGGAGGTACCgg                                                         >  1:38270/1‑151 (MQ=255)
                                                                                          ccGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGAAAAAGCGCTCCCCCAAAACCGCCCGCGCCTTCACCCGGAAAGGGCTTTTATTTTGGGGGAGGCTTATTTGCCCCGGAGccaaccg                                                         >  2:285468/1‑145 (MQ=255)
                                                                                            gACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCt                                                       <  1:1112648/151‑1 (MQ=255)
                                                                                              cGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGAGAAAAAGAACTTCCACAAAAAAGCCCGCGGCTTTTCCGGGAAGGGGCTTTTATTTTGGGGGGGGCTTTTTTTCCCCGAAGCAAACGGCCAc                                                     >  2:1256902/1‑151 (MQ=255)
                                                                                              cGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGCGTTGGAAAAACCCCCCAACAAAACCGCCACCCCCGTCCCCCCTGAAAGCCTTTGCTTTTGCGGCTGAGTTATTCTCCCCGAAGCTCCCgggtgc                                                     >  1:525013/1‑147 (MQ=255)
                                                                                              cGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTAc                                                     <  1:1140523/151‑1 (MQ=255)
                                                                                                   cGGCAGTTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCt                                                <  2:1220082/151‑1 (MQ=255)
                                                                                                    ggCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTg                                               <  2:1016023/151‑1 (MQ=255)
                                                                                                     gCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTgg                                              <  2:563358/151‑1 (MQ=255)
                                                                                                        ggTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGATGAAAAAAATCTCCAACAAAAACCCCCGCGTCTTCACCCGGAAGGGGCCTTCTTTTTGGGGCAGGCTTTTTTTCCCCGAAGATACCCGCCACGGGCTGGgtt                                           >  2:1184019/1‑151 (MQ=255)
                                                                                                        ggTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGttt                                           <  1:494883/151‑1 (MQ=255)
                                                                                                         gTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTa                                          <  1:69293/151‑1 (MQ=255)
                                                                                                          tGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGAAAAAACTCTCCCACAAAAACGCCTGCGTCTTCACCCGTAAGGGCCCTTTATTTTGCGGCAGGCTTTTTTTCCCGGGAGCTACCGGCTACGGGCTGGgtttt                                         >  1:499447/1‑149 (MQ=255)
                                                                                                              atCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTca                                     <  2:604897/151‑1 (MQ=255)
                                                                                                                cGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTcaca                                   <  2:682009/151‑1 (MQ=255)
                                                                                                                 ggggTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGATGAAAAAAATCTCCCACAATAACGCCCGCGTCTTTACCGGTAAAGGGCTTTTTTTTTGCGGCCGGCTTTTTTTCCCCGGAGCTACCGGCGACGGGCTGGGttttttttc                                    >  1:646767/1‑146 (MQ=255)
                                                                                                                      tggtggTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGc                              <  1:290757/150‑1 (MQ=255)
                                                                                                                        gtggtCGTGAAGTCGGCTTTATGGCGGGGGTGATGAAAAAACTCTCCCACAAAACCCCCCGCGGCTTCACCGGGAAAGGGCCTTTATTTTGGGGGCGGCTTTTTTGCCCGGGAGCTACCCGCTACGGGCTGGGTTTTTTTTCAGAAGCaa                            >  1:985023/1‑150 (MQ=255)
                                                                                                                             cGTGAAGTCGGCTTTATGGCGGGGGTGATGAAAAAAATCTCCCACAAAAACGCCCGCGGCTTTCCCCGGAAGGGCCTTTTATTTGGCGGGGGTCTTTTTTTCCCGGGAGCTACCGGCTCCCGGCTGGGTTTTTTTCCCGAAGCAAAGCTa                       >  1:1237615/1‑150 (MQ=255)
                                                                                                                             cGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTa                       <  2:219504/150‑1 (MQ=255)
                                                                                                                             cGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACAGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTaa                      <  1:256005/151‑1 (MQ=255)
                                                                                                                              gTGAAGTCGGCTTTATGGCGGGGGGGGTGAAAAAACTCTCCAACAAAAACGCCCGCGTCTTCTCCCGGAAGGGCCTTTTATTTTGGGGCAGTCTTATTTTCCCGGGAGCTAACCGCCACGGGCTGGGTTTTTTTTCAGAAACAAAGGTaa                      >  2:878775/1‑150 (MQ=255)
                                                                                                                                gAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAAc                   <  1:806684/151‑1 (MQ=255)
                                                                                                                                 aaGTCGGCTTTATGGCGGGGGGGATGAAAAAACTCTCCAACAAAAACGCCCGCGTCTTCACCCGGGAGGGGCCTTCATTTTGCGGCGGTCTTATTTTCCCGGAAGCTAACGGCGACGGTCTGGGTTTTTTTACAGAAGCAATGCTaaaa                    >  2:693117/1‑149 (MQ=255)
                                                                                                                                   gTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGcc                <  1:1159151/151‑1 (MQ=255)
                                                                                                                                     cGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCACCAATGCCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCa               <  2:1152375/150‑1 (MQ=255)
                                                                                                                                     cGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCa               <  1:320524/150‑1 (MQ=255)
                                                                                                                                             tGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATggg      >  1:680916/1‑151 (MQ=255)
                                                                                                                                              ggCGGGGTTGATGAAAAAGCTCTCCAACAATCCCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGt     <  1:169767/151‑1 (MQ=255)
                                                                                                                                              ggCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGNNGNNNTATTCGCCCGGAAGNNNNNNNCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATggg      <  1:110764/150‑1 (MQ=255)
                                                                                                                                                cGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGttt   <  1:923848/151‑1 (MQ=255)
                                                                                                                                                 ggggATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGtttt  >  1:337241/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
ATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGTTTT  >  NC_000913/1842765‑1843060

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: