Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,421,650 (C)7→8 coding (532/858 nt) ydaU → Rac prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,421,6431.C94.6% 319.1 / 8.6 92coding (525/858 nt)ydaURac prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base C (41/46);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

CAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAAC  >  NC_000913/1421496‑1421787
                                                                                                                                                    |                                                                                                                                                  
cAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑‑cc                                                                                                                                                 <  2:223586/150‑1 (MQ=255)
   aaaGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑‑ccccc‑c                                                                                                                                            <  1:602124/151‑1 (MQ=255)
    aaGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑‑ccccc‑cc                                                                                                                                           <  1:806800/151‑1 (MQ=255)
         aGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑‑CCCCC‑CCcgggg                                                                                                                                      <  1:1030876/151‑1 (MQ=255)
         aGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑‑CCCCC‑CCcgggg                                                                                                                                      >  1:112021/1‑151 (MQ=255)
           cAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑‑CCCCC‑CCcggggg                                                                                                                                     >  2:1227169/1‑149 (MQ=255)
                gTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCg                                                                                                                                 <  1:1041543/150‑1 (MQ=255)
                  tCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAg                                                                                                                               <  1:40388/150‑1 (MQ=255)
                    aaacaaaacaAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGgtg                                                                                                                            >  2:1284335/1‑151 (MQ=255)
                    aaacaaaacaAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGgtg                                                                                                                            <  2:1396895/151‑1 (MQ=255)
                     aacaaaacaaGAGAGAAATTTAACAGGTGATCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGgtgt                                                                                                                           >  2:1684298/1‑151 (MQ=255)
                      acaaaacaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTc                                                                                                                          <  2:484251/151‑1 (MQ=255)
                       caaaacaaGATAAAAAACGAACAAGTGGTCAAACATATGTTGAATTAGGGGTTGAACATGATGTAAACACAAATTCAAAGAAAAAAGACACAGATAAAGATCTAAAAACAGACCCCCGCCGAAAA‑‑CCCCC‑CATGCGGAATCGGGGTGTc                                                                                                                          <  1:778949/149‑1 (MQ=255)
                       caaaacaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCa                                                                                                                         <  1:1110871/151‑1 (MQ=255)
                       caaaacaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATACCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGgtgt                                                                                                                           >  1:1277213/1‑149 (MQ=255)
                        aaaacaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCa                                                                                                                         >  2:377085/1‑150 (MQ=255)
                        aaaacaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCa                                                                                                                         >  2:433192/1‑150 (MQ=255)
                           acaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCaaaa                                                                                                                      >  2:1013552/1‑150 (MQ=255)
                           acaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTCTCaaaaa                                                                                                                     <  1:1340796/151‑1 (MQ=255)
                              agagagAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGtt                                                                                                                  >  1:1658421/1‑151 (MQ=255)
                               gagagaAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGttt                                                                                                                 >  1:449770/1‑151 (MQ=255)
                                 gagaAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTg                                                                                                                >  1:1263015/1‑150 (MQ=255)
                                   gaAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGAcc                                                                                                             <  2:88398/151‑1 (MQ=255)
                                    aaaTTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGAccc                                                                                                            >  1:1066265/1‑151 (MQ=255)
                                        ttAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCctct                                                                                                         <  1:102392/150‑1 (MQ=255)
                                          aaCAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTg                                                                                                        >  1:1557679/1‑149 (MQ=255)
                                             aGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGatat                                                                                                   <  2:26986/151‑1 (MQ=255)
                                               gtgtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTa                                                                                                 <  1:1402710/151‑1 (MQ=255)
                                                   tCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACtt                                                                                              <  2:1312984/150‑1 (MQ=255)
                                                    cAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTg                                                                                            >  2:48635/1‑151 (MQ=255)
                                                    cAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTg                                                                                            <  1:730393/151‑1 (MQ=255)
                                                     aaaCAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGc                                                                                           >  1:1078664/1‑151 (MQ=255)
                                                     aaaCAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGc                                                                                           >  1:1505755/1‑151 (MQ=255)
                                                       aCAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCg                                                                                         <  1:1079846/151‑1 (MQ=255)
                                                         aGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGaa                                                                                       <  2:1107752/151‑1 (MQ=255)
                                                                gAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGAtt                                                                                <  2:991587/151‑1 (MQ=255)
                                                                gAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGAtt                                                                                <  1:1435369/151‑1 (MQ=255)
                                                                 aaGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGAttt                                                                               <  2:1435828/151‑1 (MQ=255)
                                                                  aGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTc                                                                              <  2:1409572/151‑1 (MQ=255)
                                                                    tggtgTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCt                                                                             >  1:1485001/1‑151 (MQ=255)
                                                                    tggtgTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTg                                                                            >  1:76854/1‑151 (MQ=255)
                                                                    tggtgTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAAATGGATTTCt                                                                             >  1:1689487/1‑150 (MQ=255)
                                                                     ggtgTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGt                                                                           >  2:1358403/1‑151 (MQ=255)
                                                                          ttGAACATTATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGtctc                                                                        <  1:1351655/149‑1 (MQ=255)
                                                                             aaCATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCttt                                                                   >  1:129716/1‑151 (MQ=255)
                                                                             aaCATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCttt                                                                   <  2:1013965/151‑1 (MQ=255)
                                                                                  gatgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCg                                                               <  1:433157/150‑1 (MQ=255)
                                                                                  gatgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCg                                                               <  1:269371/150‑1 (MQ=255)
                                                                                   atgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTg                                                             >  1:519272/1‑151 (MQ=255)
                                                                                   atgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTg                                                             <  1:1122446/151‑1 (MQ=255)
                                                                                   atgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTg                                                             <  2:8491/151‑1 (MQ=255)
                                                                                   atgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTg                                                             >  2:974844/1‑151 (MQ=255)
                                                                                   atgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGAGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTg                                                             >  1:881090/1‑151 (MQ=255)
                                                                                       cAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTg                                                         >  2:1693567/1‑151 (MQ=255)
                                                                                       cAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTg                                                         <  2:1360686/151‑1 (MQ=255)
                                                                                        aaCACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTgg                                                        <  2:1339178/151‑1 (MQ=255)
                                                                                        aaCACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTgg                                                        <  1:1703526/151‑1 (MQ=255)
                                                                                          cacaAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGt                                                      >  1:367787/1‑151 (MQ=255)
                                                                                           acaAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGtt                                                     <  2:1519903/151‑1 (MQ=255)
                                                                                           acaAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGtt                                                     <  2:1483984/151‑1 (MQ=255)
                                                                                             aaaGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGa                                                   <  2:506030/151‑1 (MQ=255)
                                                                                               aGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAAt                                                 <  2:1196176/151‑1 (MQ=255)
                                                                                               aGGCAACT‑ATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAAtt                                                >  1:56441/1‑151 (MQ=255)
                                                                                                  cAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTc                                              <  2:1648131/151‑1 (MQ=255)
                                                                                                  cAAATAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTc                                              >  1:438753/1‑151 (MQ=255)
                                                                                                   aaCTTATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGGCTCGCTTTGGCGTGAGTGGGTTGAATTTc                                              <  2:269588/150‑1 (MQ=255)
                                                                                                   aaCTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCg                                             <  2:1557774/151‑1 (MQ=255)
                                                                                                    aCTAATAAAGATTCAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTc                                              <  1:250098/149‑1 (MQ=255)
                                                                                                    aCTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGc                                            >  2:1088257/1‑151 (MQ=255)
                                                                                                     cTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGcc                                           >  1:944249/1‑151 (MQ=255)
                                                                                                      taataaAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGcc                                           >  1:1274914/1‑150 (MQ=255)
                                                                                                       aataaAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAg                                         <  2:1028749/151‑1 (MQ=255)
                                                                                                            aGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCAtt                                    >  1:583360/1‑151 (MQ=255)
                                                                                                             gATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTg                                   <  1:86197/151‑1 (MQ=255)
                                                                                                              aTACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGc                                  >  2:1540565/1‑151 (MQ=255)
                                                                                                               tACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCg                                 >  1:640698/1‑151 (MQ=255)
                                                                                                                  aGATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGaa                               >  1:1298644/1‑150 (MQ=255)
                                                                                                                   gATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGaaaa                             <  1:1358311/151‑1 (MQ=255)
                                                                                                                   gATAAAGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGaaaa                             <  1:1284242/151‑1 (MQ=255)
                                                                                                                        aGATCTAAAAACAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGAt                        <  1:1088172/151‑1 (MQ=255)
                                                                                                                              aaaaaCAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGTATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGa                    <  1:1349196/149‑1 (MQ=255)
                                                                                                                              aaaaaCAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGa                    <  1:866833/149‑1 (MQ=255)
                                                                                                                                 agcAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGa               <  2:286709/149‑1 (MQ=255)
                                                                                                                                 aaCAGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGa               <  1:800812/151‑1 (MQ=255)
                                                                                                                                    aGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGc             >  2:537430/1‑150 (MQ=255)
                                                                                                                                    aGATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGc             >  2:928509/1‑150 (MQ=255)
                                                                                                                                     gATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGc             >  1:277159/1‑149 (MQ=255)
                                                                                                                                     gATCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAg           <  2:261037/151‑1 (MQ=255)
                                                                                                                                      aTCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAgg          >  1:1116534/1‑151 (MQ=255)
                                                                                                                                      aTCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAgg          <  1:1360750/151‑1 (MQ=255)
                                                                                                                                       tCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAggg         >  2:1243900/1‑151 (MQ=255)
                                                                                                                                       tCCCCCCCTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAggg         >  1:1513665/1‑151 (MQ=255)
                                                                                                                                        cccccccTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGgc        >  1:694706/1‑151 (MQ=255)
                                                                                                                                             ccTAAATC‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTa    >  1:1471854/1‑150 (MQ=255)
                                                                                                                                              cTAAAT‑‑CCCCCTCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTTGAATGGAGCAGGGCGCTAAc  >  1:1689355/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑‑CCCCC‑CCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAAC  >  NC_000913/1421496‑1421787

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