Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,648,349 A→G 5.3% intergenic (+8/‑159) 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,3490AG 5.3% 464.9 / 9.4 132intergenic (+8/‑159)dicA/ydfAQin prophage; putative regulator for DicB/Qin prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (66/59);  new base G (5/2);  total (71/61)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

TTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAA  >  NC_000913/1648211‑1648497
                                                                                                                                          |                                                                                                                                                     
ttGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATc                                                                                                                                            >  1:1152913/1‑149 (MQ=255)
     cTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                       <  3:477769/149‑1 (MQ=255)
        tttAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCttt                                                                                                                                         >  2:1592862/1‑144 (MQ=255)
        tttAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCttt                                                                                                                                         <  1:1592862/144‑1 (MQ=255)
                  tGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                       >  4:472776/1‑136 (MQ=255)
                  tGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                       <  1:1351324/136‑1 (MQ=255)
                  tGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                       >  2:1351324/1‑136 (MQ=255)
                  tGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                       >  4:1478141/1‑136 (MQ=255)
                          tCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCTGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                       >  4:436292/1‑128 (MQ=255)
                          tCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                       <  3:436292/128‑1 (MQ=255)
                          tCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                       >  4:280436/1‑128 (MQ=255)
                          tCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTAAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTa                                                                                                                                       <  3:280436/128‑1 (MQ=255)
                             gAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTa                                                                                                                      >  1:660144/1‑142 (MQ=255)
                                  ggATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACGTCATCAATGAGTTATCTTTTACCACATCAATTa                                                                                                                           <  2:1188260/132‑1 (MQ=255)
                                   gATACCCAGCTCGCCGAAATGCGAGCTCGAGTAATAGACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGAGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGtt                                                                                                                       <  4:170874/135‑1 (MQ=255)
                                   gATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGtt                                                                                                                       <  3:405270/135‑1 (MQ=255)
                                   gATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGtt                                                                                                                       >  3:170874/1‑135 (MQ=255)
                                   gATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGtt                                                                                                                       >  4:405270/1‑135 (MQ=255)
                                        ccAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATc                                                                                                                                            <  2:1636265/109‑1 (MQ=255)
                                         cAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACTCATTGATGATGTTATCt                                                                                                                                           >  4:820118/1‑109 (MQ=38)
                                            cTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACtt                                                                                                <  2:1152913/149‑1 (MQ=255)
                                             tCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTg                                                                                               >  4:1344391/1‑149 (MQ=255)
                                               gCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGAc                                                                                             >  4:1873702/1‑149 (MQ=255)
                                                  gAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCCTCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAAt                                                                                                     <  4:1192814/138‑1 (MQ=255)
                                                     aTGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGatat                                                                                       >  1:718324/1‑149 (MQ=255)
                                                         gAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                                   <  3:1873702/149‑1 (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                                                                             >  1:773781/1‑149 (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)
                                                                           tCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGtt                                                                 <  2:1717023/149‑1 (MQ=255)
                                                                           tCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGtt                                                                 <  3:1344391/149‑1 (MQ=255)
                                                                           tCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGtt                                                                 <  2:773781/149‑1 (MQ=255)
                                                                             aaTAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGtt                                                                 >  3:754516/1‑147 (MQ=255)
                                                                             aaTAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGtt                                                                 <  4:754516/147‑1 (MQ=255)
                                                                              aTAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGata                                                                                        >  1:106432/1‑123 (MQ=255)
                                                                               tAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGATTTATCTTTTACCACATCAATTATGTTACCtata                                                                                                                <  1:537926/98‑1 (MQ=255)
                                                                                 aaCTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATc                                                           <  4:655718/149‑1 (MQ=255)
                                                                                  aCTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                                   >  4:1012841/1‑124 (MQ=255)
                                                                                  aCTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                                   <  3:1012841/124‑1 (MQ=255)
                                                                                    tctTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGctgagaa                                                        >  2:685947/1‑143 (MQ=255)
                                                                                    tctTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGctgagaa                                                        <  1:685947/149‑7 (MQ=255)
                                                                                        tgatgaATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAacacac                                                    <  4:1632547/149‑1 (MQ=255)
                                                                                         gaagaaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                                   <  4:960738/117‑1 (MQ=255)
                                                                                           agaaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCAc                                                                                                  <  1:1851653/100‑1 (MQ=255)
                                                                                           agaaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGc                                                 >  3:1137829/1‑149 (MQ=255)
                                                                                                tACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGgt                                            <  2:400421/149‑1 (MQ=255)
                                                                                                   tAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCaa                                                         <  2:487166/133‑1 (MQ=255)
                                                                                                   tAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                         >  2:1819790/1‑149 (MQ=255)
                                                                                                   tAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                         >  1:540327/1‑149 (MQ=255)
                                                                                                    aaaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAAt                                                                       <  1:1371849/118‑1 (MQ=255)
                                                                                                    aaaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAAt                                                                       >  2:1371849/1‑118 (MQ=255)
                                                                                                    aaaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTtc                                        <  1:1195436/149‑1 (MQ=255)
                                                                                                    aaaaGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTtc                                        <  2:540327/149‑1 (MQ=255)
                                                                                                             tCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCa                                  <  3:704247/146‑1 (MQ=255)
                                                                                                             tCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCa                                  >  4:704247/1‑146 (MQ=255)
                                                                                                                gCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAg                                                                                                               <  3:972571/66‑1 (MQ=255)
                                                                                                                  gtaCAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTc                          >  1:885844/3‑149 (MQ=255)
                                                                                                                  ggACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTc                          <  3:1454339/149‑1 (MQ=255)
                                                                                                                   gACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCaa                                                         >  4:1625867/1‑117 (MQ=255)
                                                                                                                   gACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCaa                                                         <  3:1625867/117‑1 (MQ=255)
                                                                                                                   gACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCa                                     <  1:1302346/137‑1 (MQ=255)
                                                                                                                   gACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTcc                         >  4:592151/1‑149 (MQ=255)
                                                                                                                   gACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTcc                         >  2:1031569/1‑149 (MQ=255)
                                                                                                                      aaaTAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCt                      <  1:1031569/149‑1 (MQ=255)
                                                                                                                      aaaTAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCt                      <  3:592151/149‑1 (MQ=255)
                                                                                                                        gtAAAAGATAACTCATTGATGATGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAc                    >  3:394956/2‑149 (MQ=255)
                                                                                                                         tAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                                   <  4:1264464/85‑1 (MQ=255)
                                                                                                                         tAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                                   >  3:1264464/1‑85 (MQ=255)
                                                                                                                         tAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATg                                                                             >  4:360362/1‑91 (MQ=255)
                                                                                                                         tAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATg                                                                             <  3:360362/91‑1 (MQ=255)
                                                                                                                         tAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAAt                                                                       >  2:550378/1‑97 (MQ=255)
                                                                                                                         tAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGc                                                              <  2:1239508/106‑1 (MQ=255)
                                                                                                                         tAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAcc                   >  2:1183454/1‑149 (MQ=255)
                                                                                                                         tAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAcc                   >  3:1205427/1‑149 (MQ=255)
                                                                                                                          aaaaGATAACTCATTGATGATGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAccc                  >  4:498640/1‑149 (MQ=255)
                                                                                                                           aaaGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAcccc                 >  4:124086/1‑149 (MQ=255)
                                                                                                                             aGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCtt                                                                    <  1:1210617/96‑1 (MQ=255)
                                                                                                                             aGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCtt                                                                    >  2:1210617/1‑96 (MQ=255)
                                                                                                                             aGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATACCAGACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAg               >  2:1001315/1‑149 (MQ=255)
                                                                                                                             aGATAACTCATTGATGATGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATa                                                                                                           <  2:769826/57‑1 (MQ=17)
                                                                                                                             aGATAACTCATTGATGATGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAg               >  3:1475653/1‑149 (MQ=255)
                                                                                                                              gATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGc              <  3:1521405/149‑1 (MQ=255)
                                                                                                                                    atcatcAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                         <  4:1826118/116‑1 (MQ=255)
                                                                                                                                    atcatcAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                         >  3:1826118/1‑116 (MQ=255)
                                                                                                                                    atcatcAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAg        >  4:1502561/1‑149 (MQ=255)
                                                                                                                                    atcatcAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAg        >  4:1372830/1‑149 (MQ=255)
                                                                                                                                    atcatcAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAg        <  2:885844/149‑1 (MQ=255)
                                                                                                                                     tcatcaATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCa                                     <  1:1646044/119‑1 (MQ=255)
                                                                                                                                      catcaATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAggg      <  3:1684286/149‑1 (MQ=255)
                                                                                                                                       atcaATGAGTTATCTTTTACCAAATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAgggg     >  2:713855/1‑149 (MQ=255)
                                                                                                                                          aaTGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTtctc                                      >  1:1524340/1‑113 (MQ=255)
                                                                                                                                          aaTGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTtctc                                      <  2:1524340/113‑1 (MQ=255)
                                                                                                                                          aaTGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTtctc                                      >  3:734938/1‑113 (MQ=255)
                                                                                                                                          |                                                                                                                                                     
TTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGAGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAA  >  NC_000913/1648211‑1648497

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: