Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,647,873 C→T 5.2% intergenic (‑23/‑61) dicC ← / → dicA Qin prophage; DNA‑binding transcriptional regulator for DicB/Qin prophage; putative regulator for DicB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,647,8730CT 5.2% 321.4 / 6.1 96intergenic (‑23/‑61)dicC/dicAQin prophage; DNA‑binding transcriptional regulator for DicB/Qin prophage; putative regulator for DicB
Reads supporting (aligned to +/- strand):  ref base C (45/46);  new base T (2/3);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

ACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGC  >  NC_000913/1647725‑1648019
                                                                                                                                                    |                                                                                                                                                  
aCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTc                                                                                                                                                    <  3:269744/149‑1 (MQ=255)
  gCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAt                                                                                                                                                  >  2:693840/1‑149 (MQ=255)
    aTCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAtaa                                                                                                                                                >  4:1039540/1‑149 (MQ=255)
         gcgACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTTAtgactaac                                                                                                                                           <  3:1423400/149‑1 (MQ=255)
               ttCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagata                                                                                                                                     <  1:693840/149‑1 (MQ=255)
                   ggAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATg                                                                                                                                 <  3:1039540/149‑1 (MQ=255)
                     aaCTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTTATGACTAACaat                                                                                                                                        >  1:1164385/1‑139 (MQ=255)
                     aaCTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTTATGACTAACaat                                                                                                                                        <  2:1164385/140‑2 (MQ=255)
                     aaCTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGtt                                                                                                                               >  1:381333/1‑149 (MQ=255)
                       cTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTaa                                                                                                                             <  4:727442/149‑1 (MQ=255)
                          aaTCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTaa                                                                                                                             <  3:523543/146‑1 (MQ=255)
                          aaTCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTaa                                                                                                                             >  4:523543/1‑146 (MQ=255)
                            tCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACaaaa                                                                                                                        >  1:509395/1‑149 (MQ=255)
                             ccccTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACaaaaa                                                                                                                       <  2:381333/149‑1 (MQ=255)
                              cccTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCa                                                                                                                                                   <  3:627690/120‑1 (MQ=255)
                              cccTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCa                                                                                                                                                   >  4:627690/1‑120 (MQ=255)
                                 tttCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTTCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATa                                                                                                                   <  1:574731/149‑1 (MQ=255)
                                 tttCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATa                                                                                                                   <  1:699360/149‑1 (MQ=255)
                                   tCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAtaactaa                                                                                                                                            >  4:960594/1‑122 (MQ=255)
                                   tCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAtaactaa                                                                                                                                            <  3:960594/122‑1 (MQ=255)
                                     cAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAtaa                                                                                                                                                <  1:1216441/116‑1 (MQ=255)
                                     cAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAtaa                                                                                                                                                >  2:1216441/1‑116 (MQ=255)
                                     cAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTATATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCa                                                                                                               >  3:1190919/1‑149 (MQ=255)
                                      aGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTc                                                                                                                <  3:469446/147‑1 (MQ=255)
                                      aGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTc                                                                                                                >  4:469446/1‑147 (MQ=255)
                                      aGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCaa                                                                                                              <  1:1196061/149‑1 (MQ=255)
                                       gCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATa                                                                                                                   >  1:662369/1‑143 (MQ=255)
                                       gCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATa                                                                                                                   <  2:662369/143‑1 (MQ=255)
                                        cTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATa                                                                                                            <  2:509395/149‑1 (MQ=255)
                                           taAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGa                                                                                                         >  2:85787/1‑149 (MQ=255)
                                               gCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTa                                                                                                     >  4:470700/1‑149 (MQ=255)
                                                   aGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCAt                                                                                                 <  1:85787/149‑1 (MQ=255)
                                                       aaaCGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACaa                                                                                                                                         >  1:289529/1‑105 (MQ=255)
                                                       aaaCGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACaa                                                                                                                                         <  2:289529/105‑1 (MQ=255)
                                                       aaaCGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagata                                                                                                                                     <  4:1420670/109‑1 (MQ=255)
                                                       aaaCGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagata                                                                                                                                     >  3:1420670/1‑109 (MQ=255)
                                                       aaaCGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACaaaaa                                                                                                                       >  4:248376/1‑123 (MQ=255)
                                                       aaaCGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACaaaaa                                                                                                                       <  3:248376/123‑1 (MQ=255)
                                                       aaaCGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGc                                                                                             >  1:1326184/1‑149 (MQ=255)
                                                       aaaCGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGc                                                                                             <  2:1326184/149‑1 (MQ=255)
                                                                ccTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTATATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATgg                                                                                    <  4:1190919/149‑1 (MQ=255)
                                                                                      tttttGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACaa                                                                                                                                         >  3:1156084/1‑74 (MQ=255)
                                                                                      tttttGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACaa                                                                                                                                         <  4:1156084/74‑1 (MQ=255)
                                                                                        tttGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagat                                                                                                                                      >  1:953974/1‑75 (MQ=255)
                                                                                        tttGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagat                                                                                                                                      <  4:596651/75‑1 (MQ=255)
                                                                                        tttGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagat                                                                                                                                      <  2:953974/75‑1 (MQ=255)
                                                                                        tttGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagat                                                                                                                                      >  3:596651/1‑75 (MQ=255)
                                                                                         ttGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTAt                                                                                      >  4:1110954/1‑122 (MQ=255)
                                                                                         ttGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTAt                                                                                      <  3:1110954/122‑1 (MQ=255)
                                                                                            aaCCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACg                                                        >  2:315162/1‑149 (MQ=255)
                                                                                                       aaaaGAGCGTCAGTTTTAAGCATTTAAAACATCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCAtt                                                                                                <  4:948188/98‑1 (MQ=255)
                                                                                                       aaaaGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTc                                                                                                                >  2:584165/1‑82 (MQ=255)
                                                                                                       aaaaGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTc                                                                                                                <  1:584165/82‑1 (MQ=255)
                                                                                                       aaaaGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCAtt                                                                                                >  3:948188/1‑98 (MQ=255)
                                                                                                        aaaGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTAtcg                                            >  2:1195018/1‑149 (MQ=255)
                                                                                                           gagCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTAtcgtcg                                         <  3:470700/149‑1 (MQ=255)
                                                                                                                tCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATa                                                                                                                   >  2:182853/1‑70 (MQ=255)
                                                                                                                tCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATa                                                                                                                   <  1:182853/70‑1 (MQ=255)
                                                                                                                        aaGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCAt                                                     <  4:1053680/124‑1 (MQ=255)
                                                                                                                        aaGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCAt                                                     >  3:1053680/1‑124 (MQ=255)
                                                                                                                         aGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGAtt                                                                                                       <  2:1222100/73‑1 (MQ=255)
                                                                                                                         aGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGAtt                                                                                                       >  1:1222100/1‑73 (MQ=255)
                                                                                                                            aTTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCaa                             >  2:479987/1‑144 (MQ=255)
                                                                                                                            aTTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCaa                             <  1:479987/144‑1 (MQ=255)
                                                                                                                             tttAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATa                                                                                                                   >  4:89259/1‑57 (MQ=255)
                                                                                                                             tttAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATa                                                                                                                   <  3:89259/57‑1 (MQ=255)
                                                                                                                                aaaaCACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAAc                                                                                         >  2:257546/1‑80 (MQ=255)
                                                                                                                                aaaaCACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAAc                                                                                         <  1:257546/80‑1 (MQ=255)
                                                                                                                                  aaCACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACaaaaaa                                                                          <  2:140377/93‑1 (MQ=255)
                                                                                                                                  aaCACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACaaaaaa                                                                          >  1:140377/1‑93 (MQ=255)
                                                                                                                                  aaCACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCaaa                  <  1:1195018/149‑1 (MQ=255)
                                                                                                                                   acacCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAg                 <  1:315162/149‑1 (MQ=255)
                                                                                                                                    cacCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTa                                                                                       <  3:272027/78‑1 (MQ=255)
                                                                                                                                    cacCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTa                                                                                       >  4:272027/1‑78 (MQ=255)
                                                                                                                                     acCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGa                                                           <  3:1115024/105‑1 (MQ=255)
                                                                                                                                     acCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGa                                                           >  4:1115024/1‑105 (MQ=255)
                                                                                                                                        tttATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTAt                                                                 <  2:690085/96‑1 (MQ=255)
                                                                                                                                        tttATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCgg                                        >  3:794809/1‑121 (MQ=255)
                                                                                                                                        tttATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCgg                                        <  4:794809/121‑1 (MQ=255)
                                                                                                                                        tttATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAg                 <  1:181406/144‑1 (MQ=255)
                                                                                                                                        tttATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAg                 >  2:181406/1‑144 (MQ=255)
                                                                                                                                        tttATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGtctc            >  2:216627/1‑149 (MQ=255)
                                                                                                                                        tttATTGTTAGTCATAACTAACAAGATAGAGGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTAt                                                                 >  1:690085/1‑96 (MQ=255)
                                                                                                                                           cttGTTAGTTATGACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAacacac                       >  1:706628/2‑135 (MQ=255)
                                                                                                                                           cttGTTAGTTATGACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAacacac                       <  2:706628/134‑1 (MQ=255)
                                                                                                                                           aTTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAAc                                  <  4:343122/124‑1 (MQ=255)
                                                                                                                                           aTTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAAc                                  >  3:343122/1‑124 (MQ=255)
                                                                                                                                              gTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGc        <  1:1297640/147‑1 (MQ=255)
                                                                                                                                              gTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGc        >  2:1297640/1‑147 (MQ=255)
                                                                                                                                               ttAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGa                                                                                   >  4:630056/1‑71 (MQ=255)
                                                                                                                                               ttAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGa                                                                                   <  3:630056/71‑1 (MQ=255)
                                                                                                                                               ttAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTg         <  1:91113/145‑1 (MQ=255)
                                                                                                                                               ttAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTg         >  2:91113/1‑145 (MQ=255)
                                                                                                                                               ctAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTaa     >  4:837563/2‑149 (MQ=255)
                                                                                                                                                  gTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGc  >  4:534648/1‑149 (MQ=255)
                                                                                                                                                   tCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAAc                                  >  1:831625/1‑116 (MQ=255)
                                                                                                                                                   tCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAAc                                  <  2:831625/116‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGC  >  NC_000913/1647725‑1648019

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