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.

AAAGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAG  >  NC_000913/1648195‑1648491
                                                                                                                                                    |                                                                                                                                                     
aaaGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAAca                                                                                                                                                       >  4:251094/1‑149 (MQ=255)
                ttGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATc                                                                                                                                      >  1:1152913/1‑149 (MQ=255)
                     cTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTa                                                                                                                                 <  3:477769/149‑1 (MQ=255)
                        tttAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCttt                                                                                                                                   >  2:1592862/1‑144 (MQ=255)
                                  tGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTa                                                                                                                                 <  1:1351324/136‑1 (MQ=255)
                                  tGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTa                                                                                                                                 <  3:1478141/136‑1 (MQ=255)
                                  tGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTa                                                                                                                                 >  4:1478141/1‑136 (MQ=255)
                                  tGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTa                                                                                                                                 >  2:1351324/1‑136 (MQ=255)
                                          tCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCTGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTa                                                                                                                                 >  4:436292/1‑128 (MQ=255)
                                          tCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTa                                                                                                                                 >  4:280436/1‑128 (MQ=255)
                                          tCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTa                                                                                                                                 <  3:436292/128‑1 (MQ=255)
                                          tCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTAAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTa                                                                                                                                 <  3:280436/128‑1 (MQ=255)
                                             gAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTa                                                                                                                <  2:660144/142‑1 (MQ=255)
                                                  ggATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACGTCATCAATGA‑GTTATCTTTTACCACATCAATTa                                                                                                                     <  2:1188260/132‑1 (MQ=255)
                                                   gATACCCAGCTCGCCGAAATGCGAGCTCGAGTAATAGACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGAGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGtt                                                                                                                 <  4:170874/135‑1 (MQ=255)
                                                   gATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGtt                                                                                                                 >  3:170874/1‑135 (MQ=255)
                                                   gATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGtt                                                                                                                 <  3:405270/135‑1 (MQ=255)
                                                   gATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGtt                                                                                                                 >  4:405270/1‑135 (MQ=255)
                                                        ccAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATc                                                                                                                                      <  2:1636265/109‑1 (MQ=255)
                                                        ccAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATc                                                                                                                                      >  1:1636265/1‑109 (MQ=255)
                                                        ccAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATc                                                                                              >  2:1195436/1‑149 (MQ=255)
                                                         cAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAAC‑TCATTGATGATGTTATCt                                                                                                                                     >  4:820118/1‑109 (MQ=38)
                                                         cAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAAC‑TCATTGATGATGTTATCt                                                                                                                                     <  3:820118/109‑1 (MQ=38)
                                                            cTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACtt                                                                                          <  2:1152913/149‑1 (MQ=255)
                                                             tCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTg                                                                                         >  4:1344391/1‑149 (MQ=255)
                                                                  gAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCCTCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAAt                                                                                               <  4:1192814/138‑1 (MQ=255)
                                                                  gAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAAt                                                                                               >  3:1192814/1‑138 (MQ=255)
                                                                     aTGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGatat                                                                                 >  1:718324/1‑149 (MQ=255)
                                                                         gAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                             <  3:1873702/149‑1 (MQ=255)
                                                                           gCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCtttt                                                                                                                                  >  4:830101/1‑94 (MQ=255)
                                                                           gCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCtttt                                                                                                                                  <  3:830101/94‑1 (MQ=255)
                                                                               gAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTATCACATCAATTATGTTAGCTATCGCATACAAAATCACTTGACCGATATGTTAGTCATg                                                                       >  1:1717023/1‑149 (MQ=255)
                                                                               gAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                             >  4:942276/1‑143 (MQ=255)
                                                                               gAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                             <  3:942276/143‑1 (MQ=255)
                                                                               gAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATg                                                                       <  2:718324/149‑1 (MQ=255)
                                                                                    aaaaaCTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACGAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCtata                                                                                                          >  1:907755/1‑109 (MQ=255)
                                                                                    aaaaaCTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACGAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCtata                                                                                                          <  2:907755/109‑1 (MQ=255)
                                                                                    aaaaaCTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACGAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCtata                                                                                                          <  4:1683311/109‑1 (MQ=255)
                                                                                    aaaaaCTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACGAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCtata                                                                                                          >  3:1683311/1‑109 (MQ=255)
                                                                                       aaCTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCt                                                               >  1:400421/1‑149 (MQ=255)
                                                                                          ttCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAACTAAAAGATAACATCATCAATGA‑GTTATCTTTTACcaca                                                                                                                            >  1:282480/1‑85 (MQ=255)
                                                                                          ttCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACcaca                                                                                                                            <  2:282480/85‑1 (MQ=255)
                                                                                          ttCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGt                                                            >  3:655718/1‑149 (MQ=255)
                                                                                           tCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGtt                                                           <  2:773781/149‑1 (MQ=255)
                                                                                           tCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGtt                                                           <  2:1717023/149‑1 (MQ=255)
                                                                                           tCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGtt                                                           <  3:1344391/149‑1 (MQ=255)
                                                                                             aaTAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGtt                                                           >  3:754516/1‑147 (MQ=255)
                                                                                             aaTAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGtt                                                           <  4:754516/147‑1 (MQ=255)
                                                                                              aTAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGata                                                                                  <  2:106432/123‑1 (MQ=255)
                                                                                              aTAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGata                                                                                  >  1:106432/1‑123 (MQ=255)
                                                                                               tAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑TTTATCTTTTACCACATCAATTATGTTACCtata                                                                                                          <  1:537926/98‑1 (MQ=255)
                                                                                               tAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCtata                                                                                                          >  2:537926/1‑98 (MQ=255)
                                                                                               tAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCa                                                       >  4:1454339/1‑149 (MQ=255)
                                                                                                 aaCTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATc                                                     <  4:655718/149‑1 (MQ=255)
                                                                                                    tctTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGctgagaa                                                  >  2:685947/1‑143 (MQ=255)
                                                                                                    tctTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGctgagaa                                                  <  1:685947/149‑7 (MQ=255)
                                                                                                    tctTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAac                                                  >  1:1405534/1‑149 (MQ=255)
                                                                                                      tttGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCAGACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAacac                                                >  3:1632547/1‑149 (MQ=255)
                                                                                                        tgatgaATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAacacac                                              <  4:1632547/149‑1 (MQ=255)
                                                                                                         gaagaaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                             >  3:960738/1‑117 (MQ=255)
                                                                                                         gaagaaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                             <  4:960738/117‑1 (MQ=255)
                                                                                                           agaaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCAc                                                                                            >  2:1851653/1‑100 (MQ=255)
                                                                                                           agaaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGc                                           >  3:1137829/1‑149 (MQ=255)
                                                                                                                tACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGgt                                      <  2:400421/149‑1 (MQ=255)
                                                                                                                   tAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                   >  1:540327/1‑149 (MQ=255)
                                                                                                                    aaaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAAt                                                                 >  2:1371849/1‑118 (MQ=255)
                                                                                                                    aaaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAAt                                                                 <  1:1371849/118‑1 (MQ=255)
                                                                                                                    aaaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTtc                                  <  1:1195436/149‑1 (MQ=255)
                                                                                                                    aaaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTtc                                  <  2:540327/149‑1 (MQ=255)
                                                                                                                     aaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCAc                                                                                            >  3:1267373/1‑90 (MQ=255)
                                                                                                                     aaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCAc                                                                                            <  4:1267373/90‑1 (MQ=255)
                                                                                                                     aaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTtct                                 <  1:1819790/149‑1 (MQ=255)
                                                                                                                       aGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCa                               <  4:1137829/149‑1 (MQ=255)
                                                                                                                                gCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAg                                                                                                         <  3:972571/66‑1 (MQ=255)
                                                                                                                                  gtaCAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTc                    >  1:885844/3‑149 (MQ=255)
                                                                                                                                   gACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCaa                                                   >  4:1625867/1‑117 (MQ=255)
                                                                                                                                   gACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCa                               >  2:1302346/1‑137 (MQ=255)
                                                                                                                                   gACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTcc                   >  4:592151/1‑149 (MQ=255)
                                                                                                                                      aaaTAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCt                <  3:592151/149‑1 (MQ=255)
                                                                                                                                      aaaTAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCt                <  1:1031569/149‑1 (MQ=255)
                                                                                                                                         tAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                             <  4:1264464/85‑1 (MQ=255)
                                                                                                                                         tAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATg                                                                       <  3:360362/91‑1 (MQ=255)
                                                                                                                                         tAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAAt                                                                 <  1:550378/97‑1 (MQ=255)
                                                                                                                                         tAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGc                                                        >  1:1239508/1‑106 (MQ=255)
                                                                                                                                         tAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAcc             >  3:1767451/1‑149 (MQ=255)
                                                                                                                                         tAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAcc             >  3:1205427/1‑149 (MQ=255)
                                                                                                                                          aaaaGATAAC‑TCATTGATGATGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAccc            >  4:498640/1‑149 (MQ=255)
                                                                                                                                             aGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCtt                                                              >  2:1210617/1‑96 (MQ=255)
                                                                                                                                             aGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCtt                                                              <  1:1210617/96‑1 (MQ=255)
                                                                                                                                             aGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATACCAGACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAg         >  2:1001315/1‑149 (MQ=255)
                                                                                                                                             aGATAAC‑TCATTGATGATGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATa                                                                                                     >  1:769826/1‑57 (MQ=17)
                                                                                                                                             aGATAAC‑TCATTGATGATGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAg         >  3:1475653/1‑149 (MQ=255)
                                                                                                                                              gATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGc        <  3:1521405/149‑1 (MQ=255)
                                                                                                                                                   catcatCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTc                    <  4:366063/132‑1 (MQ=255)
                                                                                                                                                   catcatCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTc                    >  3:366063/1‑132 (MQ=255)
                                                                                                                                                    atcatcAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                   <  4:1826118/116‑1 (MQ=255)
                                                                                                                                                    atcatcAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                   >  3:1826118/1‑116 (MQ=255)
                                                                                                                                                    atcatcAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAg  <  2:885844/149‑1 (MQ=255)
                                                                                                                                                    atcatcAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAg  >  4:1846274/1‑149 (MQ=255)
                                                                                                                                                    atcatcAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAg  >  4:1502561/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAAGAACTAGAGCTCCTTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAG  >  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: ‑