Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,071,345 T→A 5.7% noncoding (29/195 nt) isrC → novel sRNA, CP4‑44; putative prophage remnant

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,071,3450TA 5.7% 218.2 / 2.3 82noncoding (29/195 nt)isrCnovel sRNA, CP4‑44; putative prophage remnant
Reads supporting (aligned to +/- strand):  ref base T (43/34);  new base A (2/3);  total (45/37)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTT  >  NC_000913/2071199‑2071493
                                                                                                                                                  |                                                                                                                                                    
gATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTtata                                                                                                                                                     <  3:274997/148‑1 (MQ=255)
   tttCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTGGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATatcga                                                                                                                                                 >  8:148915/1‑149 (MQ=255)
     tcttctGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATatcgat                                                                                                                                                <  7:183916/148‑1 (MQ=255)
     tcttctGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATatcgat                                                                                                                                                <  7:389301/148‑1 (MQ=255)
           gAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAg                                                                                                                                         >  3:354788/1‑149 (MQ=255)
           gAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAg                                                                                                                                         >  2:184687/1‑149 (MQ=255)
            acccTGTCGTGAATGCTGCACTCCCTGACTCTGCTTCTGCTCACAAAAGCAATGTTTTCGTTACTGTCTCTATTGTCCGTGCAATAGCTCAATAATAGTATAAAACGATCAATATATATTTTATCGATCGATATAAACGATCGATaaaa                                                                                                                                        <  3:94981/147‑3 (MQ=255)
            aaCCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGATATAAACGATCGATa                                                                                                                                           <  1:191439/146‑1 (MQ=255)
               cTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATATGCtaa                                                                                                                                     >  2:251626/1‑149 (MQ=255)
               cTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaa                                                                                                                                     <  4:151833/149‑1 (MQ=255)
               cTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaa                                                                                                                                     >  6:217557/1‑149 (MQ=255)
                    gTGACTGATGCCCTCCCTGCCTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaataata                                                                                                                                <  1:269005/149‑1 (MQ=255)
                         tGATGCCCTCCCTGACTCTGAGTCTGCTCACAACAGCAGTGTGTTCGTTACTGTCTCTCTTGTCCCTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCCATCGTTTATATCGATCGATAAGCTAATAATAACCtt                                                                                                                           >  6:191602/1‑149 (MQ=255)
                         tGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCt                                                                                                                            <  8:151312/148‑1 (MQ=255)
                           aTGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCttt                                                                                                                          <  5:367108/148‑1 (MQ=255)
                            tGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGt                                                                                                                        <  1:251626/149‑1 (MQ=255)
                              ccctccctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCttt                                                                                                                          >  4:336700/1‑145 (MQ=255)
                              ccctccctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCttt                                                                                                                          <  3:336700/145‑1 (MQ=255)
                                 tccctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGt                                                                                                                    >  4:141509/1‑148 (MQ=255)
                                 tccctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGt                                                                                                                    >  2:278541/1‑148 (MQ=255)
                                 tccctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCACAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGt                                                                                                                    >  8:364284/1‑148 (MQ=255)
                                   cctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTaa                                                                                                                  <  1:278541/148‑1 (MQ=255)
                                       aCTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATg                                                                                                              <  4:354788/148‑1 (MQ=255)
                                           tGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGcac                                                                                                           <  1:335429/147‑1 (MQ=255)
                                            gAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGca                                                                                                            >  4:199285/1‑145 (MQ=255)
                                             aGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGa                                                                                                        <  1:184687/148‑1 (MQ=255)
                                                 tGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACg                                                                                                    >  7:376292/1‑148 (MQ=255)
                                                 tGCTCACAAAAGCACTGTTTTCGTTACTGTATCACTTGACCGTGAAATAGAACAATAATAGAATAAAACGATCAATATCTATTTTATCTATCGTTTATATCGATAGATAACCTAATAATAACAATTGTAATTAACATTCACAGATACg                                                                                                    >  5:114566/1‑148 (MQ=255)
                                                          aaGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATATCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGa                                                                                          <  7:372241/149‑1 (MQ=255)
                                                             cACTGTTTTCGTTACTGTCACTCTTGTCCGTGCAATAGCTAAATAATAGAATAAAACGATAAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAAACTTAGTCAGTAAAATGCACAGATACGTACAGAAAGACAt                                                                                       >  1:286442/1‑149 (MQ=255)
                                                              aCTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACAtt                                                                                      <  6:43791/149‑1 (MQ=255)
                                                                            tGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAg                                                                         >  6:344020/1‑148 (MQ=255)
                                                                              tctctgTTCTCCGTGCAAAAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTAAATCGATAGATAAGCTAATAATAAACTTTGTCACAAACATGCACAGATAAGTAAAGAAAGAAATTCAGGGAACAACAGAAc                                                                      >  3:129124/1‑149 (MQ=255)
                                                                              tctctcTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAAc                                                                      >  7:60809/1‑149 (MQ=255)
                                                                              tctctcTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAAc                                                                      >  3:180227/1‑149 (MQ=255)
                                                                              tctctcTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGATATAAAcgatcgat                                                                                                                                            >  6:79041/1‑79 (MQ=37)
                                                                              tctctcTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGATATAAAcgatcgat                                                                                                                                            <  5:79041/79‑1 (MQ=37)
                                                                                tctctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATAtcgatcga                                                                                                                                             >  5:140660/1‑76 (MQ=255)
                                                                                tctctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATAtcgatcga                                                                                                                                             <  6:140660/76‑1 (MQ=255)
                                                                                tctctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCa                                                                                    >  4:365205/1‑133 (MQ=255)
                                                                                tctctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCa                                                                                    <  3:365205/133‑1 (MQ=255)
                                                                                tctctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACca                                                                    >  5:201361/1‑149 (MQ=255)
                                                                                tctctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACca                                                                    >  6:164349/1‑149 (MQ=255)
                                                                                    ttGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATATGTACAGAAAGACATTCAGGGAACAACAGAACCACAAt                                                                >  3:80047/1‑149 (MQ=255)
                                                                                    ttGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACcaca                                                                  <  5:254287/147‑1 (MQ=255)
                                                                                    ttGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACaa                                                                 >  5:60499/1‑148 (MQ=255)
                                                                                      gTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAAtt                                                               >  3:232562/1‑148 (MQ=255)
                                                                                      gTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAAtt                                                               >  7:232068/1‑148 (MQ=255)
                                                                                       tCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTc                                                              <  5:217557/148‑1 (MQ=255)
                                                                                         cGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGa                                                                                          >  5:242715/1‑118 (MQ=255)
                                                                                         cGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGa                                                                                          <  6:242715/118‑1 (MQ=255)
                                                                                         cGTGCAATAGCTCAATAATAGAATAAAAAGATCATTATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGa                                                           <  2:286442/149‑1 (MQ=255)
                                                                                                  tcTCAATAATAGAATAAAACGATGAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATATTAACCTTTGTCAGTAACATGCTCAGATACGTACAGAAAGACATTCAGGGATCAACAGAACCACAATTCAGAAACTccga                                                   <  6:114566/147‑3 (MQ=255)
                                                                                                  gCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCca                                                   <  8:60809/148‑1 (MQ=255)
                                                                                                  gCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCca                                                   <  4:180227/148‑1 (MQ=255)
                                                                                                  gCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCca                                                   <  5:164349/148‑1 (MQ=255)
                                                                                                    tCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAg                                                >  8:398454/1‑149 (MQ=255)
                                                                                                     cAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGc                                               >  5:379660/1‑149 (MQ=255)
                                                                                                      aataatAGAATAAAACGATCAATCTCTATTTTATCGATCGTTTATAAAGAACGATAACCTAATAATAACCATTCTGAGTAACATGCACAGATACGTACAGAAAGACAATCACCGAACAACAGAACAAAAATTCAGAAACACCCACAGcc                                              >  3:18905/1‑149 (MQ=255)
                                                                                                      aataatAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAAtt                                                               <  8:336370/132‑1 (MQ=255)
                                                                                                      aataatAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAAtt                                                               >  7:336370/1‑132 (MQ=255)
                                                                                                        taataGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCtt                                                                                                                           >  1:349232/1‑70 (MQ=255)
                                                                                                        taataGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCtt                                                                                                                           <  2:349232/70‑1 (MQ=255)
                                                                                                        taataGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCgg                                            >  6:246469/1‑149 (MQ=255)
                                                                                                         aataGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTAc                                                                                                 >  5:213267/1‑95 (MQ=255)
                                                                                                         aataGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTAc                                                                                                 <  6:213267/95‑1 (MQ=255)
                                                                                                           taGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGAcc                                         >  8:381674/1‑149 (MQ=255)
                                                                                                             gAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTc                                       >  5:1936/1‑149 (MQ=255)
                                                                                                                tAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAAt                                                                >  7:329618/1‑121 (MQ=255)
                                                                                                                tAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAAt                                                                <  8:329618/121‑1 (MQ=255)
                                                                                                                   aaCGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCAc                                 >  4:147854/1‑149 (MQ=255)
                                                                                                                    aCGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACt                                <  6:201361/149‑1 (MQ=255)
                                                                                                                    aCGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACt                                <  8:232068/149‑1 (MQ=255)
                                                                                                                          aaTATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAAcc                          <  6:60499/149‑1 (MQ=255)
                                                                                                                           atatCTTTTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAAcc                          >  8:18805/1‑148 (MQ=255)
                                                                                                                                tATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTAc                    >  7:118791/1‑149 (MQ=255)
                                                                                                                                  ttttATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCttt                      <  7:364284/145‑1 (MQ=255)
                                                                                                                                   tttATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTg                 >  4:354347/1‑149 (MQ=255)
                                                                                                                                   tttATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTg                 >  5:195459/1‑149 (MQ=255)
                                                                                                                                     tATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGcc               >  7:159010/1‑149 (MQ=255)
                                                                                                                                         gatcgTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTa           >  6:334164/1‑149 (MQ=255)
                                                                                                                                         gatcgGTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTCACAAGCACAGATACGTACACAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGCCCTCCGGCACTGTAACCCTTTACCTGCCGGTa           <  6:116002/149‑1 (MQ=255)
                                                                                                                                          atcgTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTAt          <  7:18805/149‑1 (MQ=255)
                                                                                                                                                  tatCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGttt  >  7:155338/1‑149 (MQ=255)
                                                                                                                                                  tatCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGttt  <  5:191602/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GATTTTCTTCTGAACCTGTCGTGACTGATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTT  >  NC_000913/2071199‑2071493

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: