Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	364,526	0	C	A	41.7%	79.9 / 21.9	60	G594C (GGC→TGC)	lacZ	beta‑D‑galactosidase
Reads supporting (aligned to +/- strand):  ref base C (23/12);  new base A (0/25);  total (23/37)
Fisher's exact test for biased strand distribution p-value = 5.09e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGG‑AAAAACTGCTGCTGGTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTG  >  NC_000913/364379‑364671                                                                                                                                                       |                                                                                                                                                  aTTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGG‑AAAAACTGCTGCTGGTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCaa                                                                                                                                                 >  2:967558/1‑150 (MQ=255) aTTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGG‑AAAAACTGCTGCTGGTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCaa                                                                                                                                                 >  2:1103444/1‑150 (MQ=255)    cGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGG‑AAAAACTGCTGCTGGTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAAt                                                                                                                                              <  1:574263/150‑1 (MQ=255)       tGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGG‑AAAAACTGCTGCTGGTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCAcc                                                                                                                                          >  2:206428/1‑151 (MQ=255)            aCTTCGATGGTTTGCCCGGATAAACGGAACTGG‑AAAAACTGCTGCTGGTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGt                                                                                                                                     >  2:1045580/1‑151 (MQ=255)              ttCGATGTTTTGCCCGGAAAAACGGAACTGG‑AAAACCTGCTGCTGGTGTTTTG‑CTCCCGTCAGCGCGGGATGCGGCGTGCGGTCGGCAAAGACCAGCCCGTTCATACAGAACGGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGt                                                                                                                                     <  1:357608/149‑1 (MQ=255)                 gATGGTTTGCCCGGATAAACGGAACTGG‑AAAAACTGCTGCTGGTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCAAAAATCACCGCCGTAAGc                                                                                                                                 <  2:383641/150‑1 (MQ=255)                    ggTTTCCCCGGA‑AAAAGGAACCGGAAAAAACTGCTGCGGGTGTTTTG‑CTCCCGTCAGCGCTGGTTGCGGCGTGCGGCCGACAAAGACCAGCCCTTCCATACAGACCGGGCGAT‑CGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGAc                                                                                                                             <  1:241149/151‑1 (MQ=255)                           cggaaaaaaCGGACGGGG‑AAAATCTGCTGC‑GGTTTTTTGTCTTCCGTCAGCGCTGGAAGCGCCGGGGGGTCGGCAAAGACCAGCCCGTCCACACAG‑ACCGGCGGTGCGTCCGGGGTCCCGCAAAATTCTCCGCCGTAAGCCGACCACGGGt                                                                                                                      <  2:1042611/145‑1 (MQ=255)                           ccGGATAAACGGAACTGG‑AAAAACTGCTGCTGGTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGt                                                                                                                      >  1:1101807/1‑151 (MQ=255)                             ggATAAACGGAACTGG‑AAAAACTGCTGCTGGTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGAAAAGACCAGACCGTTCATACAAAACTGGCAAT‑CGTTCGGGGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGtt                                                                                                                     <  2:28339/150‑1 (MQ=255)                               aaaaaCGGAACTGA‑AAAAACTGCTGCTGGTTTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGCCGGAAAAGACCAGACCGTTCATACAGAACGGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTg                                                                                                                    <  1:1051430/147‑1 (MQ=255)                                 aaaCGGACCTGG‑AAAACCTGCTGCTGGTGTTTTC‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGCCCGTTCATACAGAACTGGCGAT‑GGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGt                                                                                                                <  1:530277/151‑1 (MQ=255)                                     ggAACTGG‑AAAAGCTGCTGCTGGTGTTTTG‑TTTCCGGCAGGGCGGGAGGGGGCGTGCGGCGGGAAAAGCCCAACCCGTTAAAAAAGAACTGCCGTT‑GGTTCGGCGTATCGCAAAAATCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTtc                                                                                                            <  1:355608/151‑1 (MQ=255)                                     ggAACTGG‑AAAAACTGCTGCTGGTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTtc                                                                                                            >  2:187635/1‑151 (MQ=255)                                           gg‑AAAACCTGCTGCTGGTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGAAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTACCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTcatcat                                                                                                       <  2:1092061/150‑1 (MQ=255)                                               aaaaCTGCTGTTGGTGTTTTG‑CTTCCGTCGGCGCTGGATGCGGGGTGGGGTCGGAAAAGACCAGTCCGTTCATACAGATCTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATAttt                                                                                                   <  2:459812/151‑1 (MQ=255)                                                       tgctgGTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGc                                                                                            >  2:677809/1‑150 (MQ=255)                                                        gctgGTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCGCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGc                                                                                            <  2:995766/149‑1 (MQ=255)                                                        gctgGTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCg                                                                                           >  1:850878/1‑150 (MQ=255)                                                        gctgGTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGAAAAGCCCAGCCCGTTCATACAAACCGGGGGAT‑CGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAg                                                                                             <  1:538701/148‑1 (MQ=255)                                                        gctgGTGTTTTG‑CTTCCGTCAGCGCGGGATGGGGCGTGGGGTCGGAAAAGCCCAGACCGTTCATACAGAACTGGCGAT‑GGTTCGGCTTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGa                                                                                          <  1:556679/151‑1 (MQ=255)                                                          tgGTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACt                                                                                        >  2:70617/1‑151 (MQ=255)                                                           ggTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGAAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCTTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGa                                                                                          <  2:602294/148‑1 (MQ=255)                                                           ggTGTTTTG‑CTTCCGCCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGa                                                                                          <  2:94213/148‑1 (MQ=255)                                                           ggTGTTTTG‑CTCCCGCCAGCGCGGGATGGGGCGTGCGGTCGGAAAAGACCAGCCCGTCAATACAGAACTGGCGAT‑CGTTCGGCGTATCGCAAAAACCTCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGc                                                                                            <  2:57458/146‑1 (MQ=255)                                                                tttG‑CTGCCGGCAGCGTGGGAGGGGGCGTGCGGTGGACAAAGACCAGACCGTTCATACAAAACTGGAGAG‑CTTGCGGGGTATCCCAAAAGTCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATcc                                                                                   <  2:494188/150‑1 (MQ=255)                                                                        ccGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGACCTGGCGAT‑CGTTCGGCGTACCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTc                                                                           <  2:29102/151‑1 (MQ=255)                                                                        ccGTCAGCGCGGGATGCGGCGTGGGGTCGGCAAAGACCAGACCTTTCATACAGAACGGGCGAT‑CGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCa                                                                              <  2:734833/148‑1 (MQ=255)                                                                         cgccAGCGCTGGATGCGGCGTGCGGCCGGCAAAGCCCAGCCCGTTCATACAAACCTGGCGTT‑CTTCCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCa                                                                              <  2:671690/144‑1 (MQ=255)                                                                          gTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTccc                                                                         >  2:202599/1‑151 (MQ=255)                                                                               cgcgggaggcgggtggCGTTCGGAAAAGACCAGACCATCCAAAAAAAACTGGCGTT‑CGTTCGGCTTACCCCAAAAACCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGAc                                                                     <  1:1093355/135‑1 (MQ=255)                                                                                gcgggagggggggTGCGGTCGGAAAAGACAAGCCCGTTAATAAAAAACGGGCGAT‑CGTTCGGCGTATCCCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACg                                                                    <  1:116783/138‑1 (MQ=255)                                                                                   ggATGCGGCGTGCGGTCGGCAAAGCCCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTACCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGa                                                                      <  2:270883/145‑1 (MQ=255)                                                                                   ggATGCGGCGGGCGGCGGGAAAAGACCAGACCGTTCATACAGACCTGGCGAT‑CGTGCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGaa                                                                  <  2:246766/149‑1 (MQ=255)                                                                                         ggcgTGCGGTCGGCAAAGCCCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGccc                                                           <  2:887385/150‑1 (MQ=255)                                                                                          gcgTGCGGTCGGAAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTg                                                         <  2:668188/151‑1 (MQ=255)                                                                                             tGCGGCCGGAAAAGACCAGCCCTTGCATCCAGAACGGGCGAT‑CGTTCGGCGTACCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTg                                                         <  1:1013053/148‑1 (MQ=255)                                                                                              gCGGTCGGCAAAGCCCAGACCGTCCATACAGACCTGGCGAT‑CGTTCGGCGTATCCCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTaaa                                                     <  1:702137/151‑1 (MQ=255)                                                                                              gCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTaa                                                      <  2:778024/150‑1 (MQ=255)                                                                                              gCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTaa                                                      >  1:475053/1‑150 (MQ=255)                                                                                               cGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAAc                                                    >  1:586499/1‑151 (MQ=255)                                                                                                  tCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACACAGTCCCAGACGAAGCCGCCCTGTAAACggg                                                 >  2:627550/1‑151 (MQ=255)                                                                                                    ggCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACgggg                                                >  1:628255/1‑150 (MQ=255)                                                                                                         agaccgcccctttaaTACAAACCGGGCGCT‑GGTTCGGGGTACCCCCAAAAGCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGT‑CCAGACGAAGCCGCCCTGTAAACGGGGAt                                              <  2:353190/132‑1 (MQ=255)                                                                                                            ccagaccGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGAc                                        <  2:185558/150‑1 (MQ=255)                                                                                                             cagaccGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACg                                       >  1:20680/1‑150 (MQ=255)                                                                                                              agcccGTTCATACAGACCTGGCGAT‑CGTTCGGCGTATCGCAAAAACCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGAc                                        <  1:898574/145‑1 (MQ=255)                                                                                                                 ccGTTCATACAGACCTGGCGTT‑CTTTCGGCGTATCGCCAAAACCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACg                                  <  1:380941/151‑1 (MQ=255)                                                                                                                           aaaaCGGGCGAT‑CGTTCGGGGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAg                          <  2:206489/147‑1 (MQ=255)                                                                                                                           aGACCTGGCGTT‑CGTTCGGCGTATCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTa                        <  2:94863/151‑1 (MQ=255)                                                                                                                               cTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTAttt                     >  2:897491/1‑150 (MQ=255)                                                                                                                                 ggCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGc                  >  2:427496/1‑151 (MQ=255)                                                                                                                                 ggCGAT‑CGTCCGGCTTACCGCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGc                             <  1:1075310/140‑1 (MQ=255)                                                                                                                                  gCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCg                 >  1:570795/1‑151 (MQ=255)                                                                                                                                         gTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCg           >  2:106881/1‑151 (MQ=255)                                                                                                                                             ggCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCa        >  1:538167/1‑150 (MQ=255)                                                                                                                                                 tATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACt   >  2:637498/1‑151 (MQ=255)                                                                                                                                                 tATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACt   >  2:437210/1‑151 (MQ=255)                                                                                                                                                  aTCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTg  >  2:550365/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                  ATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGG‑AAAAACTGCTGCTGGTGTTTTG‑CTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGAT‑CGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTG  >  NC_000913/364379‑364671

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