Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,109	0	A	C	44.1%	81.5 / 26.7	61	G219G (GGT→GGG)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (27/6);  new base C (0/26);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 2.32e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537960‑4538258                                                                                                                                                        |                                                                                                                                                      aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCa                                                                                                                                                       <  1:692830/150‑1 (MQ=255) aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCAc                                                                                                                                                      <  1:624639/151‑1 (MQ=255)   ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCt                                                                                                                                                    >  2:137635/1‑151 (MQ=255)       gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAt                                                                                                                                                >  1:490900/1‑151 (MQ=255)       gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAt                                                                                                                                                <  2:135069/151‑1 (MQ=255)       gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAt                                                                                                                                                >  2:102676/1‑151 (MQ=255)        cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAtt                                                                                                                                               >  2:33548/1‑151 (MQ=255)          aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTc                                                                                                                                              >  2:572361/1‑150 (MQ=255)           aaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTca                                                                                                                                             >  2:177290/1‑150 (MQ=255)              cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAc                                                                                                                                         >  1:177804/1‑151 (MQ=255)                  aGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAccgcc                                                                                                                                     >  1:894424/1‑151 (MQ=255)                      aCGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCAc                                                                                                                                 <  1:237639/151‑1 (MQ=255)                        gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCAcc                                                                                                                                >  2:896926/1‑150 (MQ=255)                           aTCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGcc                                                                                                                            >  2:639025/1‑151 (MQ=255)                                     gagaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                                  >  1:1153789/1‑151 (MQ=255)                                        aCGGGAGCAAGTTTTTTCCTTTTTAATTTATTCCCGGTAAAATAAATTAAAAATCCGGCTTCTTTTCGC‑AACCC‑TGTTTTGGCTCCCCCATAAATAAGCCAGGTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTcc                                                                                                                       <  1:922462/143‑1 (MQ=255)                                          gggAGCAAGTTTTTTCCTTTTTAAATTATTCCCGTTAAAATAAAGTAAAAATCCGCCTCCCTTCCCC‑AACCC‑TGTTTTGGCTCCCCCATTAAAAACCCAGTTTTTATCCCTTTTTTCCCCAACCGCCGCACCAGCCGTTCCGTACCAgg                                                                                                               <  1:585989/149‑1 (MQ=255)                                          gggAGAAAGTTTATCCCTTTTAAAATTATTCCCGGTAAAATAAATTAAAAATCCGGCACCCGTCCGC‑AACCC‑TGTTTGGGCTCCCCCATAAAAAACCCGGTTTTTCTCCCTTTTATTCAAACCCGCCGCACCAGCCGTTCCGTACCa                                                                                                                 <  1:910814/147‑1 (MQ=255)                                                aaGCTTATTCCTTTTAAAATTATTCCCGGTAAAATAAAGTTCAAATACGGCATCCGTCCGA‑AATCC‑TGGTTTGGCTCCGCCATTAAAAAGCCAGGTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                         <  2:931284/149‑1 (MQ=255)                                                    ttATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     >  1:276255/1‑149 (MQ=255)                                                      aTTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                                                                                 >  1:1103565/1‑151 (MQ=255)                                                         ccATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                               >  2:702739/1‑150 (MQ=255)                                                         ccATTTAAAATTATCCCGGGTAAAATCAATTTAAAAT‑CGGCCACCCTTCCCAAATCC‑TGTTTTGGTTTCGCTATTAATAACCCGGGTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                                <  1:1068073/149‑1 (MQ=255)                                                              ttAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                          >  1:96948/1‑150 (MQ=255)                                                              taaattttttcccGGTAAAATAAATTTCAAATCCGGCTTCCGTTCGC‑ACTCC‑GGTTTTGGTTCCCCCTTAAAAAACCCAGGTTTTACCCCCTTTATCCCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                          <  1:1117177/139‑1 (MQ=255)                                                                aatttattaCCGGTAAAATCAAGTTAAAACGCGGCATCCGTTCCA‑ATTCC‑GGTTTTGGCTTCCCCATTAATAACCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                       <  2:396492/148‑1 (MQ=255)                                                                aaaTTTTTCCCGGT‑AAAACAAGTTCAAAATCCGCCTCCCTTCCG‑AATCCTTGTTTTGGTTTCCCAATTAAAAACCCAGTTTTTACCCCTTTTATCCCAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAt                                                                                        <  1:636320/150‑1 (MQ=255)                                                                 aaTTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTg                                                                                      >  2:660537/1‑151 (MQ=255)                                                                      ttaCGGGTGAAATCAATTAAAAATACGGAATCAGTCGCC‑ACTCC‑TGGTTCGTCCCCCCTATTAAAAAGCCGGTTTTTATCCCCTTTCTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt                                                                                    <  1:243758/148‑1 (MQ=255)                                                                       taCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                                >  2:716584/1‑151 (MQ=255)                                                                         ccGGTAAAATAAAGTTCAAATACGGCACCCGTTCGC‑AATCC‑TGGTTTGGTTCCGCCATTAATAAGCCAGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                              <  2:839689/151‑1 (MQ=255)                                                                           ggTAAAATAAAGTCAAAATCCGCCACCCTTTCGC‑AATCC‑TGTTTTGGCTTCCCAATTAATAAGCCGGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAg                                                                            <  2:157721/151‑1 (MQ=255)                                                                            gTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg                                                                           >  2:1071112/1‑151 (MQ=255)                                                                                   cAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                     >  2:276365/1‑150 (MQ=255)                                                                                    aaGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                    >  1:1144289/1‑150 (MQ=255)                                                                                     aGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                    >  1:875505/1‑149 (MQ=255)                                                                                       taaaaaTACGGCACCCGTTCGC‑AATCC‑TGTTTTGGCTCCCCCATTAATAACCCAGTTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGa                                                                 <  2:902248/147‑1 (MQ=255)                                                                                          aaaTCCGCCATCGTTTCGA‑AATCT‑TGGTTGGGTTTCCCCATAAATAAGCCAGTTTTTATCTCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa                                                              <  2:632522/150‑1 (MQ=255)                                                                                          aaaTACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                             >  1:390166/1‑151 (MQ=255)                                                                                           aaTCCGGCACCCGTTCGC‑AACCC‑TGTTTTGGCTCCGCCATAAATAAGCCAGGTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaaca                                                            <  2:773271/151‑1 (MQ=255)                                                                                            atccGGCACCCGTTCGA‑AACCC‑GGGTTGGGTTTCGCCTTAAATAAGCCAGGTTTTACCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAg                                                           <  2:354836/148‑1 (MQ=255)                                                                                              aCGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCATAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                          >  1:153487/1‑150 (MQ=255)                                                                                                    tCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAt                                                   >  1:619339/1‑151 (MQ=255)                                                                                                          ccc‑AATCC‑TGTTTTGGCTCCCCCATTAAAAACCCGGTTTTTATCCCCTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaa                                               <  1:162626/147‑1 (MQ=255)                                                                                                            cgcacccttgttttGGTTTCCCCATAAAAAAGCCAGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                             <  1:49020/140‑1 (MQ=255)                                                                                                              aaccc‑TGGTTTGGCTTCCCCATTAATAAGCCAGTTTTTACCCCCTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                             <  1:41369/145‑1 (MQ=255)                                                                                                               atccttgttttGGTTTCGCCTTTAATAGCCCGGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaagaa                                            <  2:378784/141‑1 (MQ=255)                                                                                                                       tttGGTTCCGCCTTTAAAAAGCCGGTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                                   <  2:1151084/150‑1 (MQ=255)                                                                                                                       tttGGTTCCGCCATTAAAAAGCCAGTTTTTACCACCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                                   <  1:858915/150‑1 (MQ=255)                                                                                                                        ttGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                 >  1:643684/1‑151 (MQ=255)                                                                                                                          ggttcccccATAAAAAACCCAGGTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                                <  1:821488/143‑1 (MQ=255)                                                                                                                             ttCGCCATTAATAAGCCAGGTTTTATCCCCTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                 <  2:617780/146‑1 (MQ=255)                                                                                                                             ttCCCCATTAAGAGGCCTGGTTTTACCGCCTTAATGCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtt                              <  2:1040518/149‑1 (MQ=255)                                                                                                                               cGCCATTAATAACCCAGTTTTTACCCCCTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                          <  1:121882/151‑1 (MQ=255)                                                                                                                                  cATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                       >  2:759411/1‑151 (MQ=255)                                                                                                                                   ataaataccccggtttttAACCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                      <  2:856760/137‑1 (MQ=255)                                                                                                                                       ataaGCCAGTTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                  <  2:428482/151‑1 (MQ=255)                                                                                                                                         aggCCAGGTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa                 <  1:660524/148‑1 (MQ=255)                                                                                                                                               gttttaATCCCTTTTATTCCCCCCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc           <  1:1084055/148‑1 (MQ=255)                                                                                                                                                gTTTTATCACCTTTATTCACAACCGCCGCACAAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGTTTTTTTGTCAAAATAGTGAGCa          >  1:33013/1‑150 (MQ=255)                                                                                                                                                 ttttATCCCCTTCATCCCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt         <  1:228014/150‑1 (MQ=255)                                                                                                                                                     aTCACCTATATTCACAGCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGa      <  2:550617/149‑1 (MQ=255)                                                                                                                                                        aCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  1:186201/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                      ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGC‑AATCC‑TGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537960‑4538258

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