Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,648,450 C→T 5.4% intergenic (+109/‑58) dicA → / → ydfA Qin prophage; putative regulator for DicB/Qin prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,648,4500CT 5.4% 422.0 / 9.0 129intergenic (+109/‑58)dicA/ydfAQin prophage; putative regulator for DicB/Qin prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (62/60);  new base T (2/5);  total (64/65)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 131 total aligned reads displayed.

AGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCA  >  NC_000913/1648302‑1648598
                                                                                                                                                    |                                                                                                                                                    
  aaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCAc                                                                                                                                                    <  4:329367/149‑1 (MQ=255)
  aaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCAc                                                                                                                                                    >  4:165541/1‑149 (MQ=255)
  aaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCAc                                                                                                                                                    <  2:151916/149‑1 (MQ=255)
    ttACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCCATGAGTTATCTTTTACCAAATCAATTATGTTAACAATAGCATGCAAAATCACTTGACCGATATGTTAGTCATGCCTAATCTTGTTTGCATCAACACACCGCACgg                                                                                                                                                  >  3:667519/1‑149 (MQ=255)
    ttACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACgg                                                                                                                                                  <  4:481920/149‑1 (MQ=255)
        tAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGTGCGgt                                                                                                                                                 <  1:228693/146‑1 (MQ=255)
        tAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGTGCGgt                                                                                                                                                 >  2:228693/1‑146 (MQ=255)
        tAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                                                                                                                              >  4:952244/1‑149 (MQ=255)
        tAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                                                                                                                              <  1:915123/149‑1 (MQ=255)
        tAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGATGtt                                                                                                                                              <  4:885694/149‑1 (MQ=255)
         aaaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTtc                                                                                                                                             <  2:578487/149‑1 (MQ=255)
         aaaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTtc                                                                                                                                             <  2:55669/149‑1 (MQ=255)
          aaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGTGCGgtgtgtt                                                                                                                                            <  2:112914/149‑4 (MQ=255)
            aGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGTGCGgtgtgttga                                                                                                                                          <  3:92393/149‑6 (MQ=255)
            aGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCa                                                                                                                                          <  2:977538/149‑1 (MQ=255)
            aGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCa                                                                                                                                          <  3:165541/149‑1 (MQ=255)
             gCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGTGCGgtgtgttgat                                                                                                                                         <  3:473580/149‑7 (MQ=255)
                  tCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAAc                                                                                                                                    <  1:601073/149‑1 (MQ=255)
                  tCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAAc                                                                                                                                    >  2:284639/1‑149 (MQ=255)
                  tCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAAc                                                                                                                                    >  2:88292/1‑149 (MQ=255)
                        gACAAATAAAAGATAACTTCATCAATGAGTTATCTTTTACCCCATCAATCATGTTAGCTATAGCATACAGAATCACTTGACCGATATGTTAGGCATGGCTAATCTTGTTTGCAAAAACGCACCGCACGGTGTTCTCACCAAACAGTTcc                                                                                                                              >  4:667983/1‑149 (MQ=255)
                          cAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGc                                                                                                                            <  2:846470/149‑1 (MQ=255)
                           aaaTAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCt                                                                                                                           >  4:743164/1‑149 (MQ=255)
                           aaaTAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCt                                                                                                                           >  2:38878/1‑149 (MQ=255)
                           aaaTAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCt                                                                                                                           >  1:1002294/1‑149 (MQ=255)
                                   gATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGc                                                                                                                   <  3:169761/149‑1 (MQ=255)
                                     tAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGt                                                                                                                 >  4:501797/1‑149 (MQ=255)
                                      aaCATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGtt                                                                                                                <  2:1039128/149‑1 (MQ=255)
                                        catcatCAATGAGTTATATTTTACCACATCAATTATGTTAGCTATAGCAAACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGCTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTaa                                                                                                              <  4:911864/149‑1 (MQ=255)
                                         atcatcAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGtt                                                                                                                                >  3:133221/1‑130 (MQ=255)
                                         atcatcAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGtt                                                                                                                                <  4:133221/130‑1 (MQ=255)
                                          tcatcaATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAcc                                                                                                                        >  1:984635/1‑137 (MQ=255)
                                          tcatcaATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAgg                                                                                                            <  2:1002294/149‑1 (MQ=255)
                                                aTGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAAt                                                                                                      <  2:458641/149‑1 (MQ=255)
                                                   aGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAg                                                                                                   <  4:44515/149‑1 (MQ=255)
                                                    gTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGCAAATGAgg                                                                                                  <  4:334453/149‑1 (MQ=255)
                                                          ttttACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCa                                                                                            >  4:715575/1‑149 (MQ=255)
                                                          ttttACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCa                                                                                            <  1:38878/149‑1 (MQ=255)
                                                           tttACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCAt                                                                                           <  1:253979/149‑1 (MQ=255)
                                                                cacaTCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGgt                                                                                                                                                 <  2:68362/90‑1 (MQ=255)
                                                                cacaTCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGgt                                                                                                                                                 >  1:68362/1‑90 (MQ=255)
                                                                cacaTCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATa                                                                                      >  1:1019762/1‑149 (MQ=255)
                                                                   aTCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGTGCGgtgtgtt                                                                                                                                            <  4:1064679/92‑4 (MQ=255)
                                                                   aTCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGTGCGgtgtgtt                                                                                                                                            >  3:1064679/1‑89 (MQ=255)
                                                                   aTCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTa                                                                                   >  3:658368/1‑149 (MQ=255)
                                                                     cAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATc                                                                                 <  3:952244/149‑1 (MQ=255)
                                                                      aaTTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCg                                                                                >  1:30635/1‑149 (MQ=255)
                                                                       aTTATGTTAGCTATAGCATATAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGa                                                                               >  1:878340/1‑149 (MQ=255)
                                                                       aTTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGa                                                                               <  3:970229/149‑1 (MQ=255)
                                                                            gTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTg                                                                          <  2:98292/149‑1 (MQ=255)
                                                                             ttAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTgg                                                                         <  2:165673/149‑1 (MQ=255)
                                                                                         tACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAg                                                                                                   >  1:442450/1‑111 (MQ=255)
                                                                                         tACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAg                                                                                                   <  2:442450/111‑1 (MQ=255)
                                                                                          aCAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATg                                                                                          >  3:560758/1‑119 (MQ=255)
                                                                                          aCAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATg                                                                                          <  4:560758/119‑1 (MQ=255)
                                                                                              aaTCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTgg                                                        <  1:1049693/149‑1 (MQ=255)
                                                                                              aaTCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTgg                                                        <  2:193932/149‑1 (MQ=255)
                                                                                                    ttGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACg                                                  <  3:743164/149‑1 (MQ=255)
                                                                                                      gACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGc                                                <  4:644526/149‑1 (MQ=255)
                                                                                                       aCCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCg                                               >  1:587117/1‑149 (MQ=255)
                                                                                                        ccGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgt                                              <  1:589988/149‑1 (MQ=255)
                                                                                                        ccGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgt                                              <  2:857663/149‑1 (MQ=255)
                                                                                                         cGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgtg                                             >  3:553170/1‑149 (MQ=255)
                                                                                                         cGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgtg                                             >  2:817673/1‑149 (MQ=255)
                                                                                                         cGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgtg                                             <  2:518701/149‑1 (MQ=255)
                                                                                                          gATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACTGTGTTCTCAGCAAACAGTTACGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgtgt                                            >  4:59617/1‑149 (MQ=255)
                                                                                                          gATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgtgt                                            <  3:17165/149‑1 (MQ=255)
                                                                                                           atatGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGtt                                           >  4:799052/1‑149 (MQ=255)
                                                                                                           atatGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGtt                                           <  4:654824/149‑1 (MQ=255)
                                                                                                           atatGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGtt                                           >  1:795998/1‑149 (MQ=255)
                                                                                                           atatGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGtt                                           <  2:765734/149‑1 (MQ=255)
                                                                                                               gtaaGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCa                                       >  4:723905/4‑149 (MQ=255)
                                                                                                               gTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCa                                       >  2:30789/1‑149 (MQ=255)
                                                                                                                ttAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAAt                                                                                                      >  1:504322/1‑85 (MQ=255)
                                                                                                                ttAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCaa                                      >  2:816600/1‑149 (MQ=255)
                                                                                                                ttAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCaa                                      >  3:615488/1‑149 (MQ=255)
                                                                                                                ttAGTCATCGCTAAACTTGTTTGCATCAACATACCGCACGGGGTTCTCAGCAAACAGTTCCGCTACCCCAACGTTAAGGGGAAATGGGGTCACCATGGATACTACTGATCTTCGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCaa                                      <  3:667983/149‑1 (MQ=255)
                                                                                                                   gTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAAcc                                   <  3:723905/149‑1 (MQ=255)
                                                                                                                       tGGCTAGTCTTGTTTGCATCAACACACCGCACGCTGTTCTCAGCAAACATTTCCGCTACCCCAGCGATAAGGGGGTGTGAGGTCCGCATGGTTAATATCGCTCTTGGCAACAGGGAAACTCGGGTGTACGGGGTGTTTCCAAACCTggc                               >  2:640362/1‑148 (MQ=255)
                                                                                                                       tGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGa                               >  4:171891/1‑149 (MQ=255)
                                                                                                                           tAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGc                           >  3:1066431/1‑149 (MQ=255)
                                                                                                                             aTCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAc                         >  2:348097/1‑149 (MQ=255)
                                                                                                                                ttGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTa                                                                                   >  1:966264/1‑88 (MQ=255)
                                                                                                                                ttGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTa                                                                                   <  2:966264/88‑1 (MQ=255)
                                                                                                                                   tttGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAc                   >  2:51906/1‑149 (MQ=255)
                                                                                                                                   tttGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAc                   >  1:117779/1‑149 (MQ=255)
                                                                                                                                      gCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGCTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                >  2:192216/1‑149 (MQ=255)
                                                                                                                                        aTCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa              >  4:61847/1‑149 (MQ=255)
                                                                                                                                          cAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATg            >  2:548690/1‑149 (MQ=255)
                                                                                                                                            acacacCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc          >  1:255613/1‑149 (MQ=255)
                                                                                                                                            acacacCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc          >  4:777985/1‑149 (MQ=255)
                                                                                                                                              acacCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGt        >  1:295639/1‑149 (MQ=255)
                                                                                                                                                acCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAAt                                                                                                      >  3:738148/1‑53 (MQ=255)
                                                                                                                                                acCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAAt                                                                                                      <  4:738148/53‑1 (MQ=255)
                                                                                                                                                acCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc          >  4:780977/1‑145 (MQ=255)
                                                                                                                                                acCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc          <  3:780977/145‑1 (MQ=255)
                                                                                                                                                 ccGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGtata     <  4:161278/149‑1 (MQ=255)
                                                                                                                                                 ccGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGtata     <  4:323103/149‑1 (MQ=255)
                                                                                                                                                   gCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATAcc   <  1:284639/149‑1 (MQ=255)
                                                                                                                                                    cACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCa  <  3:862906/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCA  >  NC_000913/1648302‑1648598

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