Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,443,277	0	A	C	42.0%	169.1 / 25.7	88	V274G (GTA→GGA)	secY	preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (51/0);  new base C (0/37);  total (51/37)
Fisher's exact test for biased strand distribution p-value = 1.15e-25
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTG  >  NC_000913/3443128‑3443427                                                                                                                                                           |                                                                                                                                                       cAGACGCATAGAGTAACACATAAAGCGGTTGCACAGGTTGCAAACACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑a                                                                                                                                                        <  4:848741/150‑1 (MQ=255)  aGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑Acc                                                                                                                                                      >  1:854427/1‑151 (MQ=255)    aCGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑Acccc                                                                                                                                                    >  1:1208719/1‑151 (MQ=255)      gCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑Accccc                                                                                                                                                   >  3:66333/1‑150 (MQ=255)      gCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCg                                                                                                                                                  >  1:1050158/1‑151 (MQ=255)      gCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCg                                                                                                                                                  >  1:557112/1‑151 (MQ=255)          agagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCat                                                                                                                                              >  1:274187/1‑151 (MQ=255)          agagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCat                                                                                                                                              >  4:1230146/1‑151 (MQ=255)               aaCACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTc                                                                                                                                          >  4:1087088/1‑150 (MQ=255)                    aTAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTc                                                                                                                                    >  1:625087/1‑151 (MQ=255)                     tAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTc                                                                                                                                    >  4:191506/1‑150 (MQ=255)                             ttGCCCAGCCTGCAAATACAGCAAAATTGTTGTCAGCCGGTCCCACCCAG‑CCCCGCCCCCGACCAAGGC‑CGGGTGGGCCGCGGGGAACAAAATAATACGGGAA‑GCGAAGATTGCCGG‑GATC‑CCCCCCGCCATATTCACTTTCAGCGGTaa                                                                                                                            <  4:1123259/150‑1 (MQ=255)                               gCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAAtg                                                                                                                         >  1:1286924/1‑151 (MQ=255)                                cccAGGCTGAAAATACAGCAAATTTGTTGCCAGCCATTCCCACCCAG‑TCCCCCCCCCAAACCATGA‑CGGGTTGGTCGCCGGAAACAAAAAAATACTGAAA‑GCGAAGATTGCGGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAAtgt                                                                                                                        <  1:827314/151‑1 (MQ=255)                                 ccAGGCTGAAAATACAGAAAAATTGTTGTCGGCCATTTCCACCCAG‑TCCCCCCCCCGAACCATGA‑CGCGATGGTCGCCGGAAACAAAATAATACTGGAA‑GAGAAGATTGCCGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAAt                                                                                                                          <  2:30121/148‑1 (MQ=255)                                   aGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAAtgtgtg                                                                                                                     >  1:390615/1‑151 (MQ=255)                                    ggctcaaaaaACAGCGAATTTTTTGGCGGCCAG‑CCCCACCAGCGCCCGCCCCCGACCCGGGC‑GGCGATGGGCGCCGGAAACAGAAAAATATTGGAG‑GGAAAGTTTGCGGG‑GTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGc                                                                                                                    <  3:63912/145‑1 (MQ=255)                                    ggctaaaaaaACAGCAAATTTGTTGCCAGCCAGTCCCAACCAG‑TCCCGCCCCCGACCCAGGA‑CGGGGTGGTCGCGGGAAACAAAATAATACTGGAA‑GCGAAGATTGCGGG‑AATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAAtgtgt                                                                                                                      <  1:543163/143‑1 (MQ=255)                                     gCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                                   >  1:925306/1‑151 (MQ=255)                                        gCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtg                                                                                                                >  1:30235/1‑151 (MQ=255)                                          aaaTACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                                              >  1:1444699/1‑151 (MQ=255)                                              aCAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcagc                                                                                                          >  1:1371217/1‑151 (MQ=255)                                                aGCAAAATTGTTTTAAGCCAGTCCCACCCAG‑TCCCCCCCCCGACCCAGGG‑CGCGGTGGTCGCGGGAAACAAAATAATACTGGAG‑GCGAAGTTTGCGGG‑GTTC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcagca                                                                                                         <  4:1233621/150‑1 (MQ=255)                                                  aaaattttttgccggCCATTCCCACCCAG‑CCCCCCCCCCAACCCAAGC‑CGAGGTGGTCGCGGGAAAAAAAATAATAGTGGAA‑GTGAAGATTGCGGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAg                                                                                                      <  1:286437/137‑1 (MQ=255)                                                   gAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                        >  1:134584/1‑148 (MQ=255)                                                    aaaTTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGac                                                                                                    >  4:847329/1‑151 (MQ=255)                                                    aaaTTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGac                                                                                                    >  1:1077422/1‑151 (MQ=255)                                                     attttttgcCAGCCAGTCCCACCCAG‑CCCCCCCCCCAACCCAGGA‑CGCGTTGGTCGCGGGGAACAAAATAATACTGGAG‑GCGAAGTTGGCCGG‑GTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGaca                                                                                                   <  1:973897/146‑1 (MQ=255)                                                       tttttgcCAGCCAGTCCCAACCAG‑TCCCGCCCCCGAACCATGA‑CGCGATGGTCGCGGGAAACAGAATAATACTGGAA‑GCGAAAATTGCGGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGaca                                                                                                   <  4:836369/146‑1 (MQ=255)                                                       tttttgcCACCCAGTCCCACCCAG‑CCCCCCCCCCAACCCATGA‑CGGGATGGTCGCGGGGAACAAAATAATACTGGAA‑GCGAAGATTGCCGG‑GTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGac                                                                                                    <  4:537329/145‑1 (MQ=255)                                                            gtCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACgacgac                                                                                            >  1:273233/1‑151 (MQ=255)                                                            gtCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACgacgac                                                                                            >  4:122962/1‑151 (MQ=255)                                                                 ccAGTCACAACAAT‑TCCCCCCCCCGAACCATGA‑CGTGATGGCCGCGGGGAACAAATAAATACGGAAC‑GCGAA‑ATTGCCGGTGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                         <  4:478273/149‑1 (MQ=255)                                                                     ttCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGAc                                                                                   >  1:1302849/1‑151 (MQ=255)                                                                     ttCCAACCAG‑TACCGCCACCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGa                                                                                    >  3:62753/1‑150 (MQ=255)                                                                      tCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACg                                                                                  >  4:586378/1‑151 (MQ=255)                                                                         aaCCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGttt                                                                               >  1:548634/1‑151 (MQ=255)                                                                          acccagccccccccccGAACCATGA‑CGCGATGGTCGCGGGGAAAAAAATAATACGGGAA‑GCGAAGTTTGCGGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGt                                                                                 <  1:1115379/138‑1 (MQ=255)                                                                             aG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg                                                                           >  1:646299/1‑151 (MQ=255)                                                                                 aCCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                          >  4:1141718/1‑149 (MQ=255)                                                                                  ccGCCCCCGAACAATGA‑CGGGATGGTCGCGGGGAAAAAAATAATACGGGAA‑GCGAAGATTGCCGG‑AATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTc                                                                              <  4:903994/144‑1 (MQ=255)                                                                                  ccGCCCCCAAACCAGGA‑GGGGATGGTCGCGGGAAACAAAAAAATTCTGGAG‑GGAAAGATTGCGGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTa                                                                         <  4:1167071/149‑1 (MQ=255)                                                                                   cGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGtt                                                                      >  1:967032/1‑151 (MQ=255)                                                                                     cccccGAACAATGC‑CGGGTGGGTCGCCGGAAACAAAATAATACTGGAG‑GCAAAGTTTGCGGT‑TTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGttt                                                                     <  4:1466845/150‑1 (MQ=255)                                                                                      ccccCAACCAAGAGCGGGATGGTCGCCGGAAAAAAAAAAATGTTGAAA‑GCGAATGTGGCCGG‑TTTC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGttt                                                                     <  4:863900/150‑1 (MQ=255)                                                                                        ccaaccaaggA‑CGCGATGGCCGCCGGAAAAAAAAAAATGCTGGAG‑GGAAAGTTGCCCGG‑AATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACca                                                                 <  2:66526/142‑1 (MQ=255)                                                                                             ccATGA‑CGCGATGGTCGCCGGAAACAGAAAAATACTGGAA‑GCAAAGATTGCCGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAAt                                                             <  1:683162/150‑1 (MQ=255)                                                                                              ccgtgccggggtgttCGCCGGAAACAAAATAATTCGGGAG‑GCGAATTTTGCGGG‑TATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACaa                                                              <  4:1041061/135‑1 (MQ=255)                                                                                                tGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcgg                                                         >  3:139138/1‑151 (MQ=255)                                                                                                aggcggcgATGGTCGCCGGAAAAAGAATAATCGTGGAA‑GAGAATTTTGCCGG‑TTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgc                                                           <  4:1202927/145‑1 (MQ=255)                                                                                                  gaggggATGGTCGCGGGAAACAAAAAAAT‑GTGGAAGGCGAAGTTTGCGGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggc                                                        <  4:424576/146‑1 (MQ=255)                                                                                                     gcgATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGtt                                                     >  1:589670/1‑151 (MQ=255)                                                                                                       gATGGTGGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTg                                                    >  4:552454/1‑150 (MQ=255)                                                                                                       gATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTg                                                    >  4:747566/1‑150 (MQ=255)                                                                                                          ggccgcgGGAAACAAAAAAATACGGGAA‑GCGAAGATTGCGGG‑TATA‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGa                                                   <  1:1087172/145‑1 (MQ=255)                                                                                                              gCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGc                                             >  4:106633/1‑150 (MQ=255)                                                                                                                cGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTc                                           >  4:793036/1‑150 (MQ=255)                                                                                                                 gggAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTca                                          >  4:1163662/1‑150 (MQ=255)                                                                                                                 gggAACAAAAAAATCCGGGAA‑GGGAAGATTGCCGG‑GTTC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTc                                           <  1:206894/149‑1 (MQ=255)                                                                                                                      cAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                                    >  4:191939/1‑151 (MQ=255)                                                                                                                      cAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                                    >  1:1105692/1‑151 (MQ=255)                                                                                                                      cAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                                    >  1:714843/1‑151 (MQ=255)                                                                                                                       aaaaaaataCTGGAA‑GCAAAGATGGCCGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                                   <  2:147892/146‑1 (MQ=255)                                                                                                                        gAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac                                  >  1:1409000/1‑151 (MQ=255)                                                                                                                         aataatACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac                                  >  1:1022965/1‑150 (MQ=255)                                                                                                                          ataataCTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATaca                                 >  1:797498/1‑150 (MQ=255)                                                                                                                           taataCTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaac                               >  1:575552/1‑151 (MQ=255)                                                                                                                            aatccgggAA‑GCAAAGTTTGCCGG‑GATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaaca                              <  1:853217/145‑1 (MQ=255)                                                                                                                             ataCTGGAG‑GGGAAGTTTGCCGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                             <  1:1336436/151‑1 (MQ=255)                                                                                                                              taCTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                             >  1:412789/1‑150 (MQ=255)                                                                                                                                 tGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGa                          >  1:1052474/1‑150 (MQ=255)                                                                                                                                  ggAA‑GCGAAGTTGGCCGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                             <  1:1154614/146‑1 (MQ=255)                                                                                                                                     aaggaaaGATTGCCGG‑GATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGt                      <  1:1138930/146‑1 (MQ=255)                                                                                                                                     a‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGCAAATGTGTGCTCTGCGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc                     >  4:103708/1‑151 (MQ=255)                                                                                                                                       ggaaaGTTTGCGGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaaca                              <  4:553039/138‑1 (MQ=255)                                                                                                                                        cGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCAc                   >  2:34850/1‑151 (MQ=255)                                                                                                                                         gAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt                  >  2:105122/1‑151 (MQ=255)                                                                                                                                          aatttggCCGA‑GTTT‑CCCACCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                             <  1:713613/133‑1 (MQ=255)                                                                                                                                          aaGATTGCGGG‑GTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt                  <  1:1356011/150‑1 (MQ=255)                                                                                                                                          aaGATTGCCGG‑GATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCa                    <  1:574574/148‑1 (MQ=255)                                                                                                                                          aaGATTG‑CGG‑GATTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTg                 <  1:940131/151‑1 (MQ=255)                                                                                                                                              ttGCGGG‑AATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa               <  1:503182/149‑1 (MQ=255)                                                                                                                                                gcggg‑gATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGc                <  2:70483/143‑1 (MQ=255)                                                                                                                                                 cggg‑gATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCaa              <  1:295313/145‑1 (MQ=255)                                                                                                                                                  cGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact         >  4:241022/1‑151 (MQ=255)                                                                                                                                                   gg‑gATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatac          <  1:958557/149‑1 (MQ=255)                                                                                                                                                   gg‑gATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact         >  4:840794/1‑150 (MQ=255)                                                                                                                                                       aTT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaat      >  4:371241/1‑150 (MQ=255)                                                                                                                                                        tt‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaatac    >  1:1264631/1‑151 (MQ=255)                                                                                                                                                           aCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTg  >  1:97922/1‑151 (MQ=255)                                                                                                                                                           |                                                                                                                                                       CAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGA‑CGCGATGGTCGCCGGGAACAGAATAATACTGGAA‑GCGAAGATTGCCGG‑GATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTG  >  NC_000913/3443128‑3443427

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