Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,648,451 A→G 5.4% intergenic (+110/‑57) 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,4510AG 5.4% 428.8 / 9.0 130intergenic (+110/‑57)dicA/ydfAQin prophage; putative regulator for DicB/Qin prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (62/61);  new base G (2/5);  total (64/66)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

Only 100 of 130 total aligned reads displayed.

AATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCA  >  NC_000913/1648304‑1648598
                                                                                                                                                   |                                                                                                                                                   
aaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCAc                                                                                                                                                    <  2:151916/149‑1 (MQ=255)
aaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCAc                                                                                                                                                    >  4:165541/1‑149 (MQ=255)
aaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCAc                                                                                                                                                    <  4:329367/149‑1 (MQ=255)
  ttACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCCATGAGTTATCTTTTACCAAATCAATTATGTTAACAATAGCATGCAAAATCACTTGACCGATATGTTAGTCATGCCTAATCTTGTTTGCATCAACACACCGCACgg                                                                                                                                                  >  3:667519/1‑149 (MQ=255)
  ttACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACgg                                                                                                                                                  >  3:93416/1‑149 (MQ=255)
  ttACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACgg                                                                                                                                                  <  4:481920/149‑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:55669/149‑1 (MQ=255)
          aGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCa                                                                                                                                          <  3:165541/149‑1 (MQ=255)
          aGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGACTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                                                                                                                              <  3:140586/145‑1 (MQ=255)
          aGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGACTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                                                                                                                              >  4:140586/1‑145 (MQ=255)
           gCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGTGCGgtgtgttgat                                                                                                                                         <  3:473580/149‑7 (MQ=255)
               gTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCaaa                                                                                                                                     >  3:770224/1‑149 (MQ=255)
               gTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCaaa                                                                                                                                     >  2:819689/1‑149 (MQ=255)
                tCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAAc                                                                                                                                    >  2:88292/1‑149 (MQ=255)
                tCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAAc                                                                                                                                    <  1:601073/149‑1 (MQ=255)
                tCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAAc                                                                                                                                    <  2:164261/149‑1 (MQ=255)
                      gACAAATAAAAGATAACTTCATCAATGAGTTATCTTTTACCCCATCAATCATGTTAGCTATAGCATACAGAATCACTTGACCGATATGTTAGGCATGGCTAATCTTGTTTGCAAAAACGCACCGCACGGTGTTCTCACCAAACAGTTcc                                                                                                                              >  4:667983/1‑149 (MQ=255)
                        cAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGc                                                                                                                            <  2:846470/149‑1 (MQ=255)
                         aaaTAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCt                                                                                                                           >  1:1002294/1‑149 (MQ=255)
                         aaaTAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCt                                                                                                                           >  2:38878/1‑149 (MQ=255)
                         aaaTAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCt                                                                                                                           >  4:743164/1‑149 (MQ=255)
                           aTAAAAGATAACATCATCAATGAGTTATCATTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAc                                                                                                                         >  3:644526/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)
                                        tcatcaATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAcc                                                                                                                        >  1:984635/1‑137 (MQ=255)
                                        tcatcaATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAcc                                                                                                                        <  2:984635/137‑1 (MQ=255)
                                        tcatcaATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAgg                                                                                                            <  2:1002294/149‑1 (MQ=255)
                                          atcaATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATAAGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAgggg                                                                                                          >  1:35222/1‑149 (MQ=255)
                                              aTGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAAt                                                                                                      <  2:458641/149‑1 (MQ=255)
                                                 aGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAg                                                                                                   <  4:198616/149‑1 (MQ=255)
                                                 aGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAg                                                                                                   <  4:44515/149‑1 (MQ=255)
                                                         tttACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCAt                                                                                           <  1:253979/149‑1 (MQ=255)
                                                          ttACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATg                                                                                          <  1:271106/149‑1 (MQ=255)
                                                              cacaTCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGgt                                                                                                                                                 >  1:68362/1‑90 (MQ=255)
                                                              cacaTCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGgt                                                                                                                                                 <  2:68362/90‑1 (MQ=255)
                                                              cacaTCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATa                                                                                      >  1:1019762/1‑149 (MQ=255)
                                                                 aTCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGTGCGgtgtgtt                                                                                                                                            >  3:1064679/1‑89 (MQ=255)
                                                                 aTCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGTGCGgtgtgtt                                                                                                                                            <  4:1064679/92‑4 (MQ=255)
                                                                 aTCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTa                                                                                   >  3:658368/1‑149 (MQ=255)
                                                                   cAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATc                                                                                 <  3:952244/149‑1 (MQ=255)
                                                                    aaTTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCg                                                                                >  1:30635/1‑149 (MQ=255)
                                                                    aaTTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCg                                                                                >  2:551854/1‑149 (MQ=255)
                                                                     aTTATGTTAGCTATAGCATATAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGa                                                                               >  1:878340/1‑149 (MQ=255)
                                                                     aTTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGa                                                                               <  3:970229/149‑1 (MQ=255)
                                                                         tGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCtt                                                                           <  3:531473/149‑1 (MQ=255)
                                                                          gTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTg                                                                          <  2:98292/149‑1 (MQ=255)
                                                                           ttAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTgg                                                                         <  2:165673/149‑1 (MQ=255)
                                                                                 ataGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGcaacaa                                                                   <  2:262930/149‑1 (MQ=255)
                                                                                  taGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCaacaac                                                                  <  4:645287/149‑1 (MQ=255)
                                                                                      aTACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAAt                                                              >  1:783550/1‑149 (MQ=255)
                                                                                      aTACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAAt                                                              <  3:68917/149‑1 (MQ=255)
                                                                                       tACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAg                                                                                                   >  1:442450/1‑111 (MQ=255)
                                                                                        aCAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATg                                                                                          <  4:560758/119‑1 (MQ=255)
                                                                                          aaaaTCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATctct                                                          >  4:716283/1‑149 (MQ=255)
                                                                                            aaTCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTgg                                                        <  2:193932/149‑1 (MQ=255)
                                                                                            aaTCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTgg                                                        <  1:1049693/149‑1 (MQ=255)
                                                                                                     aCCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCg                                               >  1:587117/1‑149 (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)
                                                                                                         atatGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGtt                                           >  4:799052/1‑149 (MQ=255)
                                                                                                         atatGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGtt                                           <  4:654824/149‑1 (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)
                                                                                                              ttAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAAt                                                                                                      <  2:504322/85‑1 (MQ=255)
                                                                                                              ttAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCaa                                      >  3:615488/1‑149 (MQ=255)
                                                                                                               tAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCaaa                                     <  4:441328/149‑1 (MQ=255)
                                                                                                                 gTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAAcc                                   <  3:723905/149‑1 (MQ=255)
                                                                                                                  tCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACATCCACACTCTGGTCTTCCGCGTCGTTCCTACCCa                                  <  2:364525/149‑1 (MQ=255)
                                                                                                                     tGGCTAGTCTTGTTTGCATCAACACACCGCACGCTGTTCTCAGCAAACATTTCCGCTACCCCAGCGATAAGGGGGTGTGAGGTCCGCATGGTTAATATCGCTCTTGGCAACAGGGAAACTCGGGTGTACGGGGTGTTTCCAAACCTggc                               >  2:640362/1‑148 (MQ=255)
                                                                                                                         tAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGc                           >  3:1066431/1‑149 (MQ=255)
                                                                                                                           aTCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAc                         >  2:348097/1‑149 (MQ=255)
                                                                                                                            tCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCGTgg                                                                                         <  1:950016/84‑1 (MQ=255)
                                                                                                                            tCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATgg                                                                                         >  2:950016/1‑84 (MQ=255)
                                                                                                                              ttGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTa                                                                                   <  2:966264/88‑1 (MQ=255)
                                                                                                                                 tttGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAc                   >  1:117779/1‑149 (MQ=255)
                                                                                                                                    gCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                >  2:66164/1‑149 (MQ=255)
                                                                                                                                      aTCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa              >  4:61847/1‑149 (MQ=255)
                                                                                                                                        cAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATg            >  2:548690/1‑149 (MQ=255)
                                                                                                                                          acacacCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc          >  2:1038088/1‑149 (MQ=255)
                                                                                                                                          acacacCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc          >  4:777985/1‑149 (MQ=255)
                                                                                                                                          acacacCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc          >  1:255613/1‑149 (MQ=255)
                                                                                                                                            acacCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGt        >  1:295639/1‑149 (MQ=255)
                                                                                                                                              acCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAAt                                                                                                      <  4:738148/53‑1 (MQ=255)
                                                                                                                                              acCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAAt                                                                                                      >  3:738148/1‑53 (MQ=255)
                                                                                                                                              acCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc          >  4:780977/1‑145 (MQ=255)
                                                                                                                                              acCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc          <  3:780977/145‑1 (MQ=255)
                                                                                                                                               ccGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGtata     <  4:323103/149‑1 (MQ=255)
                                                                                                                                               ccGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGtata     <  4:161278/149‑1 (MQ=255)
                                                                                                                                                 gCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATAcc   >  1:509491/1‑149 (MQ=255)
                                                                                                                                                  cACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCa  <  3:862906/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCA  >  NC_000913/1648304‑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: