Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 345,302 C→A 17.2% intergenic (+311/‑102) yahL → / → yahM uncharacterized protein/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913345,3020CA17.2% 188.0 / 3.9 73intergenic (+311/‑102)yahL/yahMuncharacterized protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (46/13);  new base A (5/9);  total (51/22)
Fisher's exact test for biased strand distribution p-value = 3.64e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTATATCTGTATGCCAACTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTATTCAGGATGGGGCAGACCAATGAAAAAAGCCCAC  >  NC_000913/345063‑345523
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                               
ttttttttttgtttGCCACTTGGCTTTCCACTTTTTTTTTGCCCGATATTTTGTTGGGATCGGAGACCCTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGAc                                                                                                                                                                                                                                 <  1:1202164/230‑1 (MQ=255)
   gagttctgtatGCCAACTGGCATTCCATCTCTTGTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAg                                                                                                                                                                                                                               <  2:1231343/235‑1 (MQ=255)
         tttatGCCAACTGGCATTCCATCTCTTTTTTCGCCGATTTTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTTGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTTACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAata                                                                                                                                                                                                                         <  2:784560/237‑1 (MQ=255)
           tttggcAACCGTCATTCCAGCTTATTTTTCGCCGATATTTTGTTGTGATCGGCGACAATTCGCTGAAAAAAGCACCAGTATGGGTTTCGCAGCCATGCCGTGTATAAAAAATGATCTCATTCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATg                                                                                                                                                                                                                      <  1:798118/234‑1 (MQ=255)
           tatGCCAACTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCCTGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATg                                                                                                                                                                                                                      >  2:663971/1‑239 (MQ=255)
           tatGCCAACTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAatat                                                                                                                                                                                                                       >  2:109582/1‑238 (MQ=255)
           tatGCCAACTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAatat                                                                                                                                                                                                                       >  2:1055832/1‑238 (MQ=255)
           tatGCCAACTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATg                                                                                                                                                                                                                      <  2:352806/239‑1 (MQ=255)
           tatGCCAACTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACAAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGCGATTCATATGAAAAAAATCATAATTCCATCATGTTTGTTTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGCAGGAAAatat                                                                                                                                                                                                                       >  2:1059400/1‑238 (MQ=255)
             ttccAACTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGa                                                                                                                                                                                                                              <  1:186734/228‑1 (MQ=255)
             tGCCAACTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTTTGGGGTCTTTTCGTAATTTTCCGCATCGCACCAAAGCCTGACAGAAAAAATGcg                                                                                                                                                                                                                    >  1:352806/1‑238 (MQ=255)
             tGCCAACTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGa                                                                                                                                                                                                                              >  2:186734/1‑230 (MQ=255)
             tGCCAACTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGGGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCAGGTTTGTGTGGGGCCTTTTCTGTATCTTACGCATCGAACTCAAGCCTCACAGAAAAAATGCt                                                                                                                                                                                                                    >  1:436389/1‑239 (MQ=255)
             tGCCAACTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGGTTATCAAAAATGATCTCATGCAGCTGTTTTGTGACTGAGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGGGGGGTCTTTTCGGTATCTTACGCATCCCACTCAACCCCGAAAGAAAAAATgta                                                                                                                                                                                                                    >  1:1280758/1‑237 (MQ=255)
             tGCCAACTGGCATTCCATCTCTTTTTCTGCAGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATACAAAATGATCTCATGCAGATGTTTTGTAAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGAGGTGGGTCTTTTCAGTCTCTTATGCCTCGCACTAAAGCATGCAAGAAACTATgtt                                                                                                                                                                                                                    >  2:570650/1‑237 (MQ=255)
               ccAACTGGCATTCCATCCTTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGt                                                                                                                                                                                                                  <  2:555196/239‑1 (MQ=255)
                   cTGGCATTCCATCTTTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAata                                                                                                                                                                                                                         <  1:790654/229‑1 (MQ=255)
                   cTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAata                                                                                                                                                                                                                         >  2:790654/1‑229 (MQ=255)
                   cTGCCACTCAATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGAGGTTTCGCCGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAAACATAATTCCATCACGTTTGGTTGGGGCCTTTTCTGTAACTTAGCCATACCACTAAAACCTACAAGaaaaaaagcgtaaga                                                                                                                                                                                                                >  1:849977/1‑227 (MQ=255)
                                  ttttttCGCTTTTTTTTTTTTGTTATGGCGCACACTTCGTTTAAAAAAGCACTAGTGGGGTTTCTCCAGCCATGCGGTGTATAAAAAATGATCTCTGGCATATGTTTTGTGAAGGTGTTGTTTGACTTTAATATGAAAAAAATCATAATTCCATAATGTTTGTGTGGGGTTTTTTCTGTATCTTCCGCACCGCACTCAAGCCTAACAAAAAATATGTTGTAAGCCACATATCAAAAccc                                                                                                                                                                                              <  1:1055832/239‑1 (MQ=255)
                                    ttttCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCg                                                                                                                                                                                             >  2:702989/1‑238 (MQ=255)
                                      ttCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCa                                                                                                                                                                                          >  2:1110004/1‑239 (MQ=255)
                                                          aTCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCa                                                                                                                                                                      >  2:861301/1‑239 (MQ=255)
                                                          aTCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGAAAGAAAATATGCTGTACGGCTCATATCAAAAACCGCCATTAGCTCATCTGGCAGAGCa                                                                                                                                                                      >  2:836255/1‑239 (MQ=255)
                                                                    tctccgataaaaaaaGCCCCAGTTGTGGTTTCGCAGCAATGCGGTTAAAAAAAAATGATCTCATGCAGTTTTTTGTGGAATGTGTTGTTTGACATTAATATAAAAAAATTAATAATCCCATCATGTTTGGTGGGGGTTTTTTCTGTATTTTACGCACCGCACTCAAGCCTAACAAAAAATAGGTTGTAAGGTTAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc                                                                                                                                                             <  2:487286/230‑1 (MQ=255)
                                                                        cGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGCCAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAagag                                                                                                                                                                        >  2:288317/1‑224 (MQ=255)
                                                                        cGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAgg                                                                                                                                                         >  2:494063/1‑239 (MQ=255)
                                                                         gataaaaaaaGCACAGGTGGTGTTTTCGCAGCAGGGCGGTGAAAAAAAAATAATCTCAGGAAGGTTTTTTGGAATTTTGTGGGTTGACATTCATAAAAAAAAAATCATAATCCAATCATGTTTGGGTGGGTTTTTTTCTTTATTTTCCGCATCGCACTAAACCCTGAAAAAAAATTTGCTGTAAGGTTAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGa                                                                                                                                                        <  2:1091340/235‑1 (MQ=255)
                                                                                    cACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATgg                                                                                                                                             >  2:121147/1‑239 (MQ=255)
                                                                                    cACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGACGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGAATgg                                                                                                                                             >  2:1000481/1‑239 (MQ=255)
                                                                                                  tCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAg                                                                                                                                >  2:28941/1‑238 (MQ=255)
                                                                                                                           aTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAACTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTAACTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTACGACTCATATCATAAACCGCCATTAGCCCAACAGGAAGAGCAAACGACAACCAGAATTGTTTTATGGTACGGGGTCAGCGGACCCGCTGGCGGGCCCACTTTTCa                                                                                                       >  1:791503/1‑239 (MQ=255)
                                                                                                                           aTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCa                                                                                                       >  2:781792/1‑239 (MQ=255)
                                                                                                                           aTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAAACGCCATTAGCTCATCAGGAAGAGCAGAGAACAAACAGGATTATTCAATGGTACGGAGTTGGAGACCTCGGAGGGGGCACAACttttgg                                                                                                       >  1:561249/1‑237 (MQ=255)
                                                                                                                             gATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaa                                                                                                     >  2:1212210/1‑239 (MQ=255)
                                                                                                                             gATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaa                                                                                                     >  2:660982/1‑239 (MQ=255)
                                                                                                                             gATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCATCGATGGCGGGCCAACTTTTCaaa                                                                                                     >  2:634772/1‑239 (MQ=255)
                                                                                                                             gATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGAACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaa                                                                                                     >  2:745367/1‑239 (MQ=255)
                                                                                                                             gATCTCATGCAGATGTTTTGTGAATGTGCTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCCCTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGAATGGTATGGGGTTCGAGGCCGCGATGGACGCACAACTTTTCaaa                                                                                                     >  1:67209/1‑239 (MQ=255)
                                                                                                                              acctttgggagtttttttGGGAATGGGTGGGTTGACATTCAAATAAAAAAAATCATAACCCCATCTGTTTGGTGGGGGGTTTTTCCTGTATCGTACGAAGCCCACCAAATCCATAAAAAAAATTTGCTGTAGGGCTAAAAACAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaaa                                                                                                    <  2:567879/225‑1 (MQ=255)
                                                                                                                                   aggaagttGTTTGGTAAATGGTTTGGTTGACATTAATATGAAAAAAATCATAATTCCATATTGTTGGTGTGGTGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGAAAAAAAATATGCGGTAAGGCTAATATCAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtctt                                                                                               <  2:149535/232‑1 (MQ=255)
                                                                                                                                         aTGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAGGGAACGGGGTTCGAGGCCTCGATGCAGGTCCAACTTTTCAAAACTCTTCTGaa                                                                                          >  2:1223663/1‑238 (MQ=255)
                                                                                                                                         aTGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTACTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATCTCACAACCCGCCATTAGCTCATCAGGTAGAGCAGACGCAAACCAGGATAGTTGTATGGAACGGGCTTGGAGGACTCGATGGCGGTCCCACTTTTTAAAACACTTCTGAAt                                                                                         >  1:567879/1‑239 (MQ=255)
                                                                                                                                             tgtgtgAAGTTGTGGTTTCACATTAGAATAAAAAAACTCATATACCACTTATGTGTGTGGGGGGTCTTTCCAGTATCGTACGCACCTCACTCAAGCCTGACAGAAAATATGGTGTAAGGCCAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTACAAAACTCTTCTGAATCAGa                                                                                     <  2:78058/237‑1 (MQ=255)
                                                                                                                                                 ttaaaTGTGGTGGTTATAATTAATAAAAAAAAAATAATATTCCATTTTTGTTGTTGGGGGTTTTTTTTGTTTTTTTCCGCCTCGCACTCAACCCAAAAAAAAAATATGCTGAAAGCCAAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATcc                                                                                  <  2:1249930/236‑1 (MQ=255)
                                                                                                                                                    aaTGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCATGAAGAGCAGACGTCAACCAGGATTGATGTATGGTACGGGGATCGAAGCCCAGATGGCGGCCCAACTTTTAaaaaactatcgtaaacaggaacacc                                                                               >  1:638031/1‑218 (MQ=255)
                                                                                                                                                    aaTGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGCATCAGATCCCCt                                                                               >  2:1144980/1‑239 (MQ=255)
                                                                                                                                                    aaTGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGAGGGCGGCCCAACTTTTCAAAAACATTCTGAATAAGATCCCCt                                                                               >  1:459055/1‑239 (MQ=255)
                                                                                                                                                          tttgttctattttttataAAAAAAACTATAACCCCATCTGGTTTGGGGGGGTTTTTTCTTGTATCTGACACACCAACACAAACCTAGAAAAAAAATATGAGAAAGGGCACAAATAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTtattat                                                                         <  1:660982/225‑1 (MQ=255)
                                                                                                                                                              tggactttaaaaaaaaaaaaaTATTATGTCCATCGTGTTGGTGGGGGTTCTTTTCTGTATCTTCCGCCCCGCACTCAAGCCGAAAAAAAAATTTGTTGTGAGGCCAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatcat                                                                      <  2:776416/224‑1 (MQ=255)
                                                                                                                                                                            aaaaaaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTTTTGTATGGTACGGGGTTCGAGGTCTCGACGGCGGCCCAACTTTTCAAAAATATTCTGAATACGATCCCCTTATTCTCATCCCTGCAAAGGGGa                                                        >  2:933934/1‑238 (MQ=255)
                                                                                                                                                                            aaaaaaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATACCTGCAAAGTGGaa                                                       >  2:1240350/1‑239 (MQ=255)
                                                                                                                                                                            aaaaaaaTAAAAATCCCATGTTTTTGGGGGGGGTTTTTTTATGTATCTCACGCATCGCACACAAGCCTGACAAAAAATAGGCTGAAGGGCTAATAACAAAGCCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGaa                                                       <  2:496371/239‑1 (MQ=255)
                                                                                                                                                                              aaaaaTCATAATTCCATCATGTTTGGGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGCCGACAACCAGGATTGTTGAATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTTTTATCATCCCTGCAAAGTGGcaca                                                     >  2:1254052/1‑239 (MQ=255)
                                                                                                                                                                                        aTTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGTCGGTCCAACTTTTCAAAACTCTTTTGATTCAGATCCCCTTTTTATCATTCCTGCAAAAGGGAAAAtttctata                                             >  1:496371/1‑231 (MQ=255)
                                                                                                                                                                                             tatatgttttGGGTGGGTTTTTTTCGTTATCTAAAGCATCGCACTCAAGCCTACAAGAAAATAGGCGTGAAGCTTCTTTATAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCTACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTgg                                       <  2:397666/232‑1 (MQ=255)
                                                                                                                                                                                               tcatGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCt                                     >  1:129795/1‑239 (MQ=255)
                                                                                                                                                                                                    tttGTGGTGGGGTGTTTTCTGTATTCTCCCGTACGCACCCAAGCCCGAAAAAAAATATTGGTGAAGGCCTAAATAAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTAtt                                  <  2:1280758/239‑1 (MQ=255)
                                                                                                                                                                                                      tgtgtgGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTACCAACTCTTCTGATTCAGATCCCCTTATTATCATCCCTCCAAAGTGGAAACTTGCCacatgacccttatacacc                              >  1:591606/1‑222 (MQ=255)
                                                                                                                                                                                                        gtgtgGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTTTTTGGCTATTCAGGa                             >  2:854768/1‑239 (MQ=255)
                                                                                                                                                                                                        gtgggGGGTTTTTTCTTTATCATACGCATCGACCTAAAGCTTAAAAAAAATGTTGCTTAAAGTCTAATAAAAAAAGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTATTCAGGa                             <  2:436389/239‑1 (MQ=255)
                                                                                                                                                                                                              gggTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTATTCAGGATgggg                        >  2:335390/1‑239 (MQ=255)
                                                                                                                                                                                                               ggTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCACTTGTTTGGCTATTAAGGATGGGGc                       >  1:13173/1‑239 (MQ=255)
                                                                                                                                                                                                                   ttttCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTATTCAGGAGGGGGCAGa                    >  2:994348/1‑238 (MQ=255)
                                                                                                                                                                                                                    tttCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGGTCGAGGCCTCGAGGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTTCAAAGAGGCACATTTCAATTGTTTGGCTATTCAGAATGGGGCAGAca                  >  1:1296015/1‑238 (MQ=255)
                                                                                                                                                                                                                    tttCTGTATCTTACGCATCGCACTCAAACCTGACAGAAAATATTCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCACCAGGAATAGCAGACGACAACCAGGCTTGTTGTCTGGTACGGGGTTCGCAACATCGATGGCGGTCCACCTTTTCAAAACTCTTCTGAATACGATCCCCTTCTTATCATCCCCGCACCGTGTAACCTTGCCATTGTTTgctgtttccggagagggacgacc                  >  1:979040/1‑217 (MQ=255)
                                                                                                                                                                                                                    tttCTGTAACTTACGCATCGCTCTCAGGCCTGACAGAAAATATGCTGTAAAGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGAAGACGACAAACAGGATTGTTGTGTGGTACAGGGTTCCAGTGCTCGTAGACGCACCAACCATTGAAAACACTCCGCAGTCAGATCCCCTCATTTTCATCCCTGCAACCAGGAACCTTGCAATGGATTGGCTATTACGGATGGGGCAGAc                   >  1:657936/1‑238 (MQ=255)
                                                                                                                                                                                                                             aTCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTCCGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGa                                                                                     >  2:752703/1‑165 (MQ=255)
                                                                                                                                                                                                                               tCTTACGCATCGCACTCAAGCCTGACAGAAAATATg                                                                                                                                                                                                                      <  2:740202/36‑1 (MQ=255)
                                                                                                                                                                                                                                  tACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTCCAAAGTGGACCATTGCCCTTGTTTTGCCCTTCAGGTATGGGCAAACCAAGAaaaaaa        >  1:874407/1‑238 (MQ=255)
                                                                                                                                                                                                                                       acccgacttcaaaccataCAAAAAATATTGTGAAAGGCTAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAGGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTATTCAGGATGGGGCAGACCAATGAAAAAAGCCca   <  2:791503/222‑1 (MQ=255)
                                                                                                                                                                                                                                       aTCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTATTCAGGa                             >  2:201823/1‑212 (MQ=255)
                                                                                                                                                                                                                                       aTCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATAAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCCATTCAGGATGGGGCAGACCAGTGAAAAAAGACcac  >  2:525790/1‑239 (MQ=255)
                                                                                                                                                                                                                                         cGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCt                                                                                                                                                                                                            >  1:740202/1‑36 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                               
CTGTATATCTGTATGCCAACTGGCATTCCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCCATTGTTTGGCTATTCAGGATGGGGCAGACCAATGAAAAAAGCCCAC  >  NC_000913/345063‑345523

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