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‑‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAAC  >  NC_000913/1421496‑1421787
                                                                                                                                                    |                                                                                                                                                  
cAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑‑cc                                                                                                                                                 <  2:223586/150‑1 (MQ=255)
   aaaGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑‑cccccc                                                                                                                                            <  1:602124/151‑1 (MQ=255)
    aaGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑‑ccccccc                                                                                                                                           <  1:806800/151‑1 (MQ=255)
         aGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑‑CCCCCCCcgggg                                                                                                                                      <  1:1030876/151‑1 (MQ=255)
         aGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑‑CCCCCCCcgggg                                                                                                                                      >  1:112021/1‑151 (MQ=255)
           cAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑‑CCCCCCCcggggg                                                                                                                                     >  2:1227169/1‑149 (MQ=255)
                gTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCg                                                                                                                                 <  1:1041543/150‑1 (MQ=255)
                  tCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAg                                                                                                                               <  1:40388/150‑1 (MQ=255)
                    aaacaaaacaAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGgtg                                                                                                                            >  2:1284335/1‑151 (MQ=255)
                    aaacaaaacaAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGgtg                                                                                                                            <  2:1396895/151‑1 (MQ=255)
                     aacaaaacaaGAGAGAAATTTAACAGGTGATCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGgtgt                                                                                                                           >  2:1684298/1‑151 (MQ=255)
                      acaaaacaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTc                                                                                                                          <  2:484251/151‑1 (MQ=255)
                       caaaacaaGATAAAAAACGAACAAGTGGTCAAACATATGTTGAATTAGGGGTTGAACATGATGTAAACACAAATTCAAAGAAAAAAGACACAGATAAAGATCTAAAAACAGACCCCCGCCGAAAA‑‑CCCCCCATGCGGAATCGGGGTGTc                                                                                                                          <  1:778949/149‑1 (MQ=255)
                       caaaacaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCa                                                                                                                         <  1:1110871/151‑1 (MQ=255)
                       caaaacaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATACCCCCCTAAATCCCCCCCCGGGGGAATCGAGgtgt                                                                                                                           >  1:1277213/1‑149 (MQ=255)
                        aaaacaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCa                                                                                                                         >  2:377085/1‑150 (MQ=255)
                        aaaacaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCa                                                                                                                         >  2:433192/1‑150 (MQ=255)
                           acaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCaaaa                                                                                                                      >  2:1013552/1‑150 (MQ=255)
                           acaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTCTCaaaaa                                                                                                                     <  1:1340796/151‑1 (MQ=255)
                              agagagAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGtt                                                                                                                  >  1:1658421/1‑151 (MQ=255)
                               gagagaAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGttt                                                                                                                 >  1:449770/1‑151 (MQ=255)
                                 gagaAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTg                                                                                                                >  1:1263015/1‑150 (MQ=255)
                                   gaAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGAcc                                                                                                             <  2:88398/151‑1 (MQ=255)
                                    aaaTTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGAccc                                                                                                            >  1:1066265/1‑151 (MQ=255)
                                        ttAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCctct                                                                                                         <  1:102392/150‑1 (MQ=255)
                                          aaCAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTg                                                                                                        >  1:1557679/1‑149 (MQ=255)
                                             aGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGatat                                                                                                   <  2:26986/151‑1 (MQ=255)
                                               gtgtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTa                                                                                                 <  1:1402710/151‑1 (MQ=255)
                                                   tCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACtt                                                                                              <  2:1312984/150‑1 (MQ=255)
                                                    cAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTg                                                                                            >  2:48635/1‑151 (MQ=255)
                                                    cAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTg                                                                                            <  1:730393/151‑1 (MQ=255)
                                                     aaaCAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGc                                                                                           >  1:1078664/1‑151 (MQ=255)
                                                     aaaCAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGc                                                                                           >  1:1505755/1‑151 (MQ=255)
                                                       aCAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCg                                                                                         <  1:1079846/151‑1 (MQ=255)
                                                         aGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGaa                                                                                       <  2:1107752/151‑1 (MQ=255)
                                                                gAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGAtt                                                                                <  2:991587/151‑1 (MQ=255)
                                                                gAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGAtt                                                                                <  1:1435369/151‑1 (MQ=255)
                                                                 aaGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGAttt                                                                               <  2:1435828/151‑1 (MQ=255)
                                                                  aGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTc                                                                              <  2:1409572/151‑1 (MQ=255)
                                                                    tggtgTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCt                                                                             >  1:1485001/1‑151 (MQ=255)
                                                                    tggtgTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTg                                                                            >  1:76854/1‑151 (MQ=255)
                                                                    tggtgTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAAATGGATTTCt                                                                             >  1:1689487/1‑150 (MQ=255)
                                                                     ggtgTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGt                                                                           >  2:1358403/1‑151 (MQ=255)
                                                                          ttGAACATTATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGtctc                                                                        <  1:1351655/149‑1 (MQ=255)
                                                                             aaCATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCttt                                                                   >  1:129716/1‑151 (MQ=255)
                                                                             aaCATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCttt                                                                   <  2:1013965/151‑1 (MQ=255)
                                                                                  gatgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCg                                                               <  1:433157/150‑1 (MQ=255)
                                                                                  gatgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCg                                                               <  1:269371/150‑1 (MQ=255)
                                                                                   atgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTg                                                             >  1:519272/1‑151 (MQ=255)
                                                                                   atgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTg                                                             <  1:1122446/151‑1 (MQ=255)
                                                                                   atgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTg                                                             <  2:8491/151‑1 (MQ=255)
                                                                                   atgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTg                                                             >  2:974844/1‑151 (MQ=255)
                                                                                   atgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGAGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTg                                                             >  1:881090/1‑151 (MQ=255)
                                                                                       cAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTg                                                         >  2:1693567/1‑151 (MQ=255)
                                                                                       cAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTg                                                         <  2:1360686/151‑1 (MQ=255)
                                                                                        aaCACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTgg                                                        <  2:1339178/151‑1 (MQ=255)
                                                                                        aaCACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTgg                                                        <  1:1703526/151‑1 (MQ=255)
                                                                                          cacaAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGt                                                      >  1:367787/1‑151 (MQ=255)
                                                                                           acaAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGtt                                                     <  2:1519903/151‑1 (MQ=255)
                                                                                           acaAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGtt                                                     <  2:1483984/151‑1 (MQ=255)
                                                                                             aaaGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGa                                                   <  2:506030/151‑1 (MQ=255)
                                                                                               aGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAAt                                                 <  2:1196176/151‑1 (MQ=255)
                                                                                               aGGCAACTATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAAtt                                                >  1:56441/1‑151 (MQ=255)
                                                                                                  cAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTc                                              <  2:1648131/151‑1 (MQ=255)
                                                                                                  cAAATAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTc                                              >  1:438753/1‑151 (MQ=255)
                                                                                                   aaCTTATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGGCTCGCTTTGGCGTGAGTGGGTTGAATTTc                                              <  2:269588/150‑1 (MQ=255)
                                                                                                   aaCTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCg                                             <  2:1557774/151‑1 (MQ=255)
                                                                                                    aCTAATAAAGATTCAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTc                                              <  1:250098/149‑1 (MQ=255)
                                                                                                    aCTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGc                                            >  2:1088257/1‑151 (MQ=255)
                                                                                                     cTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGcc                                           >  1:944249/1‑151 (MQ=255)
                                                                                                      taataaAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGcc                                           >  1:1274914/1‑150 (MQ=255)
                                                                                                       aataaAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAg                                         <  2:1028749/151‑1 (MQ=255)
                                                                                                            aGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCAtt                                    >  1:583360/1‑151 (MQ=255)
                                                                                                             gATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTg                                   <  1:86197/151‑1 (MQ=255)
                                                                                                              aTACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGc                                  >  2:1540565/1‑151 (MQ=255)
                                                                                                               tACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCg                                 >  1:640698/1‑151 (MQ=255)
                                                                                                                  aGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGaa                               >  1:1298644/1‑150 (MQ=255)
                                                                                                                   gATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGaaaa                             <  1:1358311/151‑1 (MQ=255)
                                                                                                                   gATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGaaaa                             <  1:1284242/151‑1 (MQ=255)
                                                                                                                        aGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGAt                        <  1:1088172/151‑1 (MQ=255)
                                                                                                                              aaaaaCAGATCCCCCCCTAAATCCCCCCCCGGGGGTATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGa                    <  1:1349196/149‑1 (MQ=255)
                                                                                                                              aaaaaCAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGa                    <  1:866833/149‑1 (MQ=255)
                                                                                                                                 agcAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGa               <  2:286709/149‑1 (MQ=255)
                                                                                                                                 aaCAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGa               <  1:800812/151‑1 (MQ=255)
                                                                                                                                    aGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGc             >  2:537430/1‑150 (MQ=255)
                                                                                                                                    aGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGc             >  2:928509/1‑150 (MQ=255)
                                                                                                                                     gATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGc             >  1:277159/1‑149 (MQ=255)
                                                                                                                                     gATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAg           <  2:261037/151‑1 (MQ=255)
                                                                                                                                      aTCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAgg          >  1:1116534/1‑151 (MQ=255)
                                                                                                                                      aTCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAgg          <  1:1360750/151‑1 (MQ=255)
                                                                                                                                       tCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAggg         >  2:1243900/1‑151 (MQ=255)
                                                                                                                                       tCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAggg         >  1:1513665/1‑151 (MQ=255)
                                                                                                                                        cccccccTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGgc        >  1:694706/1‑151 (MQ=255)
                                                                                                                                             ccTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTa    >  1:1471854/1‑150 (MQ=255)
                                                                                                                                              cTAAAT‑‑CCCCCTCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTTGAATGGAGCAGGGCGCTAAc  >  1:1689355/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAAC  >  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: