Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,648,343 Δ1 bp 5.6% intergenic (+2/‑165) 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,3430A. 5.6% 451.0 / 17.9 124intergenic (+2/‑165)dicA/ydfAQin prophage; putative regulator for DicB/Qin prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (61/55);  new base . (5/2);  total (66/58)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

Only 100 of 136 total aligned reads displayed.

AAAGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAG  >  NC_000913/1648195‑1648491
                                                                                                                                                    |                                                                                                                                                     
aaaGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAAca                                                                                                                                                       >  4:251094/1‑149 (MQ=255)
                ttGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATc                                                                                                                                      >  1:1152913/1‑149 (MQ=255)
                     cTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                 <  3:477769/149‑1 (MQ=255)
                        tttAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCttt                                                                                                                                   >  2:1592862/1‑144 (MQ=255)
                                  tGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                 <  1:1351324/136‑1 (MQ=255)
                                  tGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                 <  3:1478141/136‑1 (MQ=255)
                                  tGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                 >  4:1478141/1‑136 (MQ=255)
                                  tGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                 >  2:1351324/1‑136 (MQ=255)
                                          tCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCTGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                 >  4:436292/1‑128 (MQ=255)
                                          tCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                 >  4:280436/1‑128 (MQ=255)
                                          tCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                 <  3:436292/128‑1 (MQ=255)
                                          tCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTAAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                 <  3:280436/128‑1 (MQ=255)
                                             gAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTa                                                                                                                <  2:660144/142‑1 (MQ=255)
                                                  ggATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACGTCATCAATGAGTTATCTTTTACCACATCAATTa                                                                                                                     <  2:1188260/132‑1 (MQ=255)
                                                   gATACCCAGCTCGCCGAAATGCGAGCTCGAGTAATAGACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGAGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGtt                                                                                                                 <  4:170874/135‑1 (MQ=255)
                                                   gATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGtt                                                                                                                 >  3:170874/1‑135 (MQ=255)
                                                   gATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGtt                                                                                                                 <  3:405270/135‑1 (MQ=255)
                                                   gATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGtt                                                                                                                 >  4:405270/1‑135 (MQ=255)
                                                        ccAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATc                                                                                                                                      <  2:1636265/109‑1 (MQ=255)
                                                        ccAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATc                                                                                                                                      >  1:1636265/1‑109 (MQ=255)
                                                        ccAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATc                                                                                              >  2:1195436/1‑149 (MQ=255)
                                                         cAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACTCATTGATGATGTTATCt                                                                                                                                     >  4:820118/1‑109 (MQ=38)
                                                         cAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACTCATTGATGATGTTATCt                                                                                                                                     <  3:820118/109‑1 (MQ=38)
                                                            cTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACtt                                                                                          <  2:1152913/149‑1 (MQ=255)
                                                             tCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTg                                                                                         >  4:1344391/1‑149 (MQ=255)
                                                                  gAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCCTCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAAt                                                                                               <  4:1192814/138‑1 (MQ=255)
                                                                  gAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAAt                                                                                               >  3:1192814/1‑138 (MQ=255)
                                                                     aTGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGatat                                                                                 >  1:718324/1‑149 (MQ=255)
                                                                         gAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                             <  3:1873702/149‑1 (MQ=255)
                                                                           gCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCtttt                                                                                                                                  >  4:830101/1‑94 (MQ=255)
                                                                           gCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCtttt                                                                                                                                  <  3:830101/94‑1 (MQ=255)
                                                                               gAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTATCACATCAATTATGTTAGCTATCGCATACAAAATCACTTGACCGATATGTTAGTCATg                                                                       >  1:1717023/1‑149 (MQ=255)
                                                                               gAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                             >  4:942276/1‑143 (MQ=255)
                                                                               gAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                             <  3:942276/143‑1 (MQ=255)
                                                                               gAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATg                                                                       <  2:718324/149‑1 (MQ=255)
                                                                                    aaaaaCTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACGAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCtata                                                                                                          >  1:907755/1‑109 (MQ=255)
                                                                                    aaaaaCTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACGAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCtata                                                                                                          <  2:907755/109‑1 (MQ=255)
                                                                                    aaaaaCTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACGAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCtata                                                                                                          <  4:1683311/109‑1 (MQ=255)
                                                                                    aaaaaCTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACGAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCtata                                                                                                          >  3:1683311/1‑109 (MQ=255)
                                                                                       aaCTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCt                                                               >  1:400421/1‑149 (MQ=255)
                                                                                          ttCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAACTAAAAGATAACATCATCAATGAGTTATCTTTTACcaca                                                                                                                            >  1:282480/1‑85 (MQ=255)
                                                                                          ttCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACcaca                                                                                                                            <  2:282480/85‑1 (MQ=255)
                                                                                          ttCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGt                                                            >  3:655718/1‑149 (MQ=255)
                                                                                           tCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGtt                                                           <  2:773781/149‑1 (MQ=255)
                                                                                           tCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGtt                                                           <  2:1717023/149‑1 (MQ=255)
                                                                                           tCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGtt                                                           <  3:1344391/149‑1 (MQ=255)
                                                                                             aaTAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGtt                                                           >  3:754516/1‑147 (MQ=255)
                                                                                             aaTAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGtt                                                           <  4:754516/147‑1 (MQ=255)
                                                                                              aTAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGata                                                                                  <  2:106432/123‑1 (MQ=255)
                                                                                              aTAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGata                                                                                  >  1:106432/1‑123 (MQ=255)
                                                                                               tAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGATTTATCTTTTACCACATCAATTATGTTACCtata                                                                                                          <  1:537926/98‑1 (MQ=255)
                                                                                               tAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCtata                                                                                                          >  2:537926/1‑98 (MQ=255)
                                                                                               tAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCa                                                       >  4:1454339/1‑149 (MQ=255)
                                                                                                 aaCTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATc                                                     <  4:655718/149‑1 (MQ=255)
                                                                                                    tctTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGctgagaa                                                  >  2:685947/1‑143 (MQ=255)
                                                                                                    tctTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGctgagaa                                                  <  1:685947/149‑7 (MQ=255)
                                                                                                    tctTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAac                                                  >  1:1405534/1‑149 (MQ=255)
                                                                                                      tttGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCAGACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAacac                                                >  3:1632547/1‑149 (MQ=255)
                                                                                                        tgatgaATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAacacac                                              <  4:1632547/149‑1 (MQ=255)
                                                                                                         gaagaaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                             >  3:960738/1‑117 (MQ=255)
                                                                                                         gaagaaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                             <  4:960738/117‑1 (MQ=255)
                                                                                                           agaaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCAc                                                                                            >  2:1851653/1‑100 (MQ=255)
                                                                                                           agaaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGc                                           >  3:1137829/1‑149 (MQ=255)
                                                                                                                tACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGgt                                      <  2:400421/149‑1 (MQ=255)
                                                                                                                   tAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                   >  1:540327/1‑149 (MQ=255)
                                                                                                                    aaaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAAt                                                                 >  2:1371849/1‑118 (MQ=255)
                                                                                                                    aaaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAAt                                                                 <  1:1371849/118‑1 (MQ=255)
                                                                                                                    aaaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTtc                                  <  1:1195436/149‑1 (MQ=255)
                                                                                                                    aaaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTtc                                  <  2:540327/149‑1 (MQ=255)
                                                                                                                     aaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCAc                                                                                            >  3:1267373/1‑90 (MQ=255)
                                                                                                                     aaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCAc                                                                                            <  4:1267373/90‑1 (MQ=255)
                                                                                                                     aaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTtct                                 <  1:1819790/149‑1 (MQ=255)
                                                                                                                       aGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCa                               <  4:1137829/149‑1 (MQ=255)
                                                                                                                                gCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAg                                                                                                         <  3:972571/66‑1 (MQ=255)
                                                                                                                                  gtaCAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTc                    >  1:885844/3‑149 (MQ=255)
                                                                                                                                   gACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCaa                                                   >  4:1625867/1‑117 (MQ=255)
                                                                                                                                   gACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCa                               >  2:1302346/1‑137 (MQ=255)
                                                                                                                                   gACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTcc                   >  4:592151/1‑149 (MQ=255)
                                                                                                                                      aaaTAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCt                <  3:592151/149‑1 (MQ=255)
                                                                                                                                      aaaTAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCt                <  1:1031569/149‑1 (MQ=255)
                                                                                                                                         tAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                             <  4:1264464/85‑1 (MQ=255)
                                                                                                                                         tAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATg                                                                       <  3:360362/91‑1 (MQ=255)
                                                                                                                                         tAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAAt                                                                 <  1:550378/97‑1 (MQ=255)
                                                                                                                                         tAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGc                                                        >  1:1239508/1‑106 (MQ=255)
                                                                                                                                         tAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAcc             >  3:1767451/1‑149 (MQ=255)
                                                                                                                                         tAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAcc             >  3:1205427/1‑149 (MQ=255)
                                                                                                                                          aaaaGATAACTCATTGATGATGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAccc            >  4:498640/1‑149 (MQ=255)
                                                                                                                                             aGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCtt                                                              >  2:1210617/1‑96 (MQ=255)
                                                                                                                                             aGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCtt                                                              <  1:1210617/96‑1 (MQ=255)
                                                                                                                                             aGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATACCAGACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAg         >  2:1001315/1‑149 (MQ=255)
                                                                                                                                             aGATAACTCATTGATGATGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATa                                                                                                     >  1:769826/1‑57 (MQ=17)
                                                                                                                                             aGATAACTCATTGATGATGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAg         >  3:1475653/1‑149 (MQ=255)
                                                                                                                                              gATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGc        <  3:1521405/149‑1 (MQ=255)
                                                                                                                                                   catcatCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTc                    <  4:366063/132‑1 (MQ=255)
                                                                                                                                                   catcatCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTc                    >  3:366063/1‑132 (MQ=255)
                                                                                                                                                    atcatcAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                   <  4:1826118/116‑1 (MQ=255)
                                                                                                                                                    atcatcAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                   >  3:1826118/1‑116 (MQ=255)
                                                                                                                                                    atcatcAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAg  <  2:885844/149‑1 (MQ=255)
                                                                                                                                                    atcatcAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAg  >  4:1846274/1‑149 (MQ=255)
                                                                                                                                                    atcatcAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAg  >  4:1502561/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAAGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAG  >  NC_000913/1648195‑1648491

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: