Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,315,516 Δ1 bp 6.0% coding (255/501 nt) yciF ← putative rubrerythrin/ferritin‑like metal‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,315,5160A. 6.0% 214.2 / 6.8 67coding (255/501 nt)yciFputative rubrerythrin/ferritin‑like metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base . (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

CGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCA  >  NC_000913/1315372‑1315660
                                                                                                                                                |                                                                                                                                                
cGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGACATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGc                                                                                                                                                <  1:370884/147‑1 (MQ=255)
        tACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAAt                                                                                                                                       >  1:321840/1‑148 (MQ=255)
         aCGGTAACCTAATTGTTCAGCCAGCGTCGCCAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATa                                                                                                                                      <  8:122173/148‑1 (MQ=255)
           ggTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGcttctt                                                                                                                                           >  6:264668/1‑141 (MQ=255)
           ggTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGcttctt                                                                                                                                           <  5:264668/141‑1 (MQ=255)
           ggTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGa                                                                                                                                   <  4:343516/149‑1 (MQ=255)
                 cTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCCCGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTc                                                                                                                              >  5:244842/1‑148 (MQ=255)
                    aTTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTg                                                                                                                         >  5:19819/1‑149 (MQ=255)
                           agccagcGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCcaca                                                                                                                    >  5:70856/1‑148 (MQ=255)
                           agccagcGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACGTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCcaca                                                                                                                    >  8:92274/1‑148 (MQ=255)
                                  gTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCAt                                                                                                            <  7:114992/149‑1 (MQ=255)
                                    cGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATgc                                                                                                          >  3:92937/1‑149 (MQ=255)
                                             cccGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAAtttt                                                                                                >  1:373858/1‑149 (MQ=255)
                                             cccGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAAtttt                                                                                                >  1:379333/1‑149 (MQ=255)
                                               cGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTGCTTCAATAAGCCCTTCCATTGCCACACATTTCATGCGCTTAATTTTc                                                                                               <  3:74575/149‑1 (MQ=255)
                                                gTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAAtttt                                                                                                >  5:360078/1‑147 (MQ=255)
                                                       gCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGa                                                                                       >  1:191203/1‑149 (MQ=255)
                                                         aaTCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGAtt                                                                                     <  6:70856/149‑1 (MQ=255)
                                                           tctcATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCg                                                                                   >  7:261548/1‑149 (MQ=255)
                                                           tctcATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCAGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCg                                                                                   >  1:208520/1‑149 (MQ=255)
                                                             tcATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCAATTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTGCTTCAATAAGACGTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGa                                                                                 >  4:18425/1‑149 (MQ=255)
                                                               aTAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCg                                                                                   <  2:321840/145‑1 (MQ=255)
                                                               aTAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACt                                                                               >  8:299542/1‑149 (MQ=255)
                                                                tAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCAGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTc                                                                              <  2:208520/149‑1 (MQ=255)
                                                                      tCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGa                                                                                       >  2:146930/1‑134 (MQ=255)
                                                                      tCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGa                                                                                       <  1:146930/134‑1 (MQ=255)
                                                                      tCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTAGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATATCATGCGCTTAATTTTCAGTTTCGATTCGGACTCCACAAc                                                                        >  8:126476/1‑149 (MQ=255)
                                                                       cGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAAc                                                                        >  5:307136/1‑148 (MQ=255)
                                                                       cGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAAc                                                                        <  6:307136/148‑1 (MQ=255)
                                                                        gACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACtt                                                                      <  8:97461/149‑1 (MQ=255)
                                                                             tCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGAt                                                                   <  7:299542/147‑1 (MQ=255)
                                                                               tgtgCTGTGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCaa                                                                >  6:79555/1‑148 (MQ=255)
                                                                                 tgctgcGGCAATCCGTCCGGCATCACGCCCTTCGTGTTGCGCGGTACTATCGAGGACCTCATACGCTTCTTCAATCAGACCTTCCCTTGCCACACATTTCATCCGCATAATATTCACGCTCGATTCGGACTCCACAACTTGATCAATAc                                                             <  6:301948/149‑1 (MQ=255)
                                                                                   ctgcGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATc                                                                  <  6:45614/142‑1 (MQ=255)
                                                                                   ctgcGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATc                                                                  >  5:45614/1‑142 (MQ=255)
                                                                                       ggCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCATCCATTGCCACACATTTCATGCGATTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCaa                                                       <  7:92274/149‑1 (MQ=255)
                                                                                        gCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCTGACTCCACAACTTGATCAATACGTTCagt                                                      >  7:342632/1‑147 (MQ=255)
                                                                                         cAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATc                                                     >  2:140566/1‑149 (MQ=255)
                                                                                            tCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTTAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACTACTTGATCAATACGTTCAATCTgg                                                  <  1:110994/149‑2 (MQ=255)
                                                                                              aGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTcc                                                <  6:244842/149‑1 (MQ=255)
                                                                                                  cGGCATCACTCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATg                                             <  7:255030/148‑1 (MQ=255)
                                                                                                  cGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGa                                            <  3:241516/149‑1 (MQ=255)
                                                                                                  cGGCATCACGCACTTCGTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGa                                            <  6:19819/148‑1 (MQ=255)
                                                                                                   ggCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCAt                                              >  4:316342/1‑146 (MQ=255)
                                                                                                   ggCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGa                                            <  3:316342/148‑1 (MQ=255)
                                                                                                      aTCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTc                                       <  2:379333/149‑1 (MQ=255)
                                                                                                      aTCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTc                                       <  2:373858/149‑1 (MQ=255)
                                                                                                         ccGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAATACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGAtt                                                                                     <  2:200918/100‑1 (MQ=255)
                                                                                                         aCGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCt                                     <  1:204184/149‑1 (MQ=255)
                                                                                                           gCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTc                                    >  5:279696/1‑148 (MQ=255)
                                                                                                             aCTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAg                                 >  2:259351/1‑149 (MQ=255)
                                                                                                               ttCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGt                               >  1:174709/1‑149 (MQ=255)
                                                                                                               ttCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGt                               >  2:302922/1‑149 (MQ=255)
                                                                                                                  gTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTGCGATTCGCACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTgcg                            >  8:320451/1‑149 (MQ=255)
                                                                                                                     tttCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCa                          <  1:259351/148‑1 (MQ=255)
                                                                                                                         tcGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATTCGCTTAATTTTCAGGTTCGATTCGGCCTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGaa                      >  5:167996/1‑148 (MQ=255)
                                                                                                                         tcGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGaaa                     >  2:92275/1‑149 (MQ=255)
                                                                                                                           ggTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAg                   >  1:253782/1‑149 (MQ=255)
                                                                                                                            gTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGc                  <  1:302922/149‑1 (MQ=255)
                                                                                                                               cTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCt                <  6:360078/148‑1 (MQ=255)
                                                                                                                               cTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCt                <  8:261548/148‑1 (MQ=255)
                                                                                                                                    gATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACtt           >  2:46173/1‑148 (MQ=255)
                                                                                                                                      tGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAAt        <  2:174709/149‑1 (MQ=255)
                                                                                                                                          cTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAAttttt    >  3:288411/1‑149 (MQ=255)
                                                                                                                                          cTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAAttttt    <  4:92937/149‑1 (MQ=255)
                                                                                                                                            cATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCa  >  2:324325/1‑149 (MQ=255)
                                                                                                                                             aTTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTc   <  1:140566/147‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCA  >  NC_000913/1315372‑1315660

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