Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,648,348 C→T 5.5% intergenic (+7/‑160) 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,3480CT 5.5% 429.1 / 9.9 126intergenic (+7/‑160)dicA/ydfAQin prophage; putative regulator for DicB/Qin prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (62/57);  new base T (5/2);  total (67/59)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 137 total aligned reads displayed.

TTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGG  >  NC_000913/1648211‑1648494
                                                                                                                                         |                                                                                                                                                   
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                                                                                                                                    <  3:472776/136‑1 (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)
                    ccAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGt                                                                                                                     <  3:251094/149‑1 (MQ=255)
                          tCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCTGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTa                                                                                                                                    >  4:436292/1‑128 (MQ=255)
                          tCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTa                                                                                                                                    >  4:280436/1‑128 (MQ=255)
                          tCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTAAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTa                                                                                                                                    <  3:280436/128‑1 (MQ=255)
                             gAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTa                                                                                                                   >  1:660144/1‑142 (MQ=255)
                                  ggATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACGTCATCAATGA‑GTTATCTTTTACCACATCAATTa                                                                                                                        <  2:1188260/132‑1 (MQ=255)
                                   gATACCCAGCTCGCCGAAATGCGAGCTCGAGTAATAGACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGAGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGtt                                                                                                                    <  4:170874/135‑1 (MQ=255)
                                   gATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGtt                                                                                                                    <  3:405270/135‑1 (MQ=255)
                                   gATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGtt                                                                                                                    >  3:170874/1‑135 (MQ=255)
                                   gATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGtt                                                                                                                    >  4:405270/1‑135 (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)
                                             tCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTg                                                                                            >  4:1344391/1‑149 (MQ=255)
                                               gCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGAc                                                                                          >  4:1873702/1‑149 (MQ=255)
                                                  gAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCCTCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAAt                                                                                                  <  4:1192814/138‑1 (MQ=255)
                                                  gAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAAt                                                                                                  >  3:1192814/1‑138 (MQ=255)
                                                               gAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                                <  3:942276/143‑1 (MQ=255)
                                                               gAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                                >  4:942276/1‑143 (MQ=255)
                                                               gAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATg                                                                          <  2:718324/149‑1 (MQ=255)
                                                               gAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATg                                                                          >  1:773781/1‑149 (MQ=255)
                                                                    aaaaaCTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACGAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCtata                                                                                                             >  3:1683311/1‑109 (MQ=255)
                                                                          ttCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAACTAAAAGATAACATCATCAATGA‑GTTATCTTTTACcaca                                                                                                                               >  1:282480/1‑85 (MQ=255)
                                                                          ttCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACcaca                                                                                                                               <  2:282480/85‑1 (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                                                              <  4:754516/147‑1 (MQ=255)
                                                                             aaTAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGtt                                                              >  3:754516/1‑147 (MQ=255)
                                                                              aTAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGata                                                                                     >  1:106432/1‑123 (MQ=255)
                                                                              aTAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGata                                                                                     <  2:106432/123‑1 (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)
                                                                                  aCTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                                >  4:1012841/1‑124 (MQ=255)
                                                                                  aCTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                                <  3:1012841/124‑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                                                                                               <  1:1851653/100‑1 (MQ=255)
                                                                                           agaaTTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGc                                              >  3:1137829/1‑149 (MQ=255)
                                                                                                tACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGgt                                         <  2:1405534/149‑1 (MQ=255)
                                                                                                tACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGgt                                         <  2:400421/149‑1 (MQ=255)
                                                                                                   tAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCaa                                                      >  1:487166/1‑133 (MQ=255)
                                                                                                   tAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                      >  2:1819790/1‑149 (MQ=255)
                                                                                                   tAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                      >  4:1684286/1‑149 (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)
                                                                                                             tCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCa                               <  3:704247/146‑1 (MQ=255)
                                                                                                             tCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCa                               >  4:704247/1‑146 (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)
                                                                                                                   gACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTcc                      >  2:1031569/1‑149 (MQ=255)
                                                                                                                      aaaTAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCt                   <  1:1031569/149‑1 (MQ=255)
                                                                                                                      aaaTAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCt                   <  3:592151/149‑1 (MQ=255)
                                                                                                                        gtAAAAGATAAC‑TCATTGATGATGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAc                 >  3:394956/2‑149 (MQ=255)
                                                                                                                         tAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                                <  4:1264464/85‑1 (MQ=255)
                                                                                                                         tAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTa                                                                                >  3:1264464/1‑85 (MQ=255)
                                                                                                                         tAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATg                                                                          <  3:360362/91‑1 (MQ=255)
                                                                                                                         tAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATg                                                                          >  4:360362/1‑91 (MQ=255)
                                                                                                                         tAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAAt                                                                    >  2:550378/1‑97 (MQ=255)
                                                                                                                         tAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGc                                                           <  2:1239508/106‑1 (MQ=255)
                                                                                                                         tAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGc                                                           >  1:1239508/1‑106 (MQ=255)
                                                                                                                         tAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAcc                >  2:1183454/1‑149 (MQ=255)
                                                                                                                         tAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAcc                >  3:1767451/1‑149 (MQ=255)
                                                                                                                          aaaaGATAAC‑TCATTGATGATGTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTAccc               >  4:498640/1‑149 (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)
                                                                                                                                    atcatcAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                      >  3:1826118/1‑116 (MQ=255)
                                                                                                                                    atcatcAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGtt                                      <  4:1826118/116‑1 (MQ=255)
                                                                                                                                    atcatcAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAg     <  2:885844/149‑1 (MQ=255)
                                                                                                                                    atcatcAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAg     >  4:1502561/1‑149 (MQ=255)
                                                                                                                                     tcatcaATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCa                                  <  1:1646044/119‑1 (MQ=255)
                                                                                                                                     tcatcaATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCa                                  >  2:1646044/1‑119 (MQ=255)
                                                                                                                                       atcaATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAgggg  >  4:525586/1‑149 (MQ=255)
                                                                                                                                       atcaATGA‑GTTATCTTTTACCAAATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAgggg  >  2:713855/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                   
TTGAGCTGTTTAATGCACTGCCAGAATCAGAACAGGATACCCAGCTCGCCGAAATGCGAGCTCGAGTAAAAAACTTCAATAAACTCTTTGAAGAATTACTAAAAGCCCGTCAGCGGACAAATAAAAGATAACATCATCAATGA‑GTTATCTTTTACCACATCAATTATGTTAGCTATAGCATACAAAATCACTTGACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGG  >  NC_000913/1648211‑1648494

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: ‑