Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,013,761 G→T 6.1% P208T (CCT→ACT)  potF‑I ← putrescine‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,013,7610GT 6.1% 208.5 / 3.1 66P208T (CCT→ACT) potF‑Iputrescine‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base T (2/2);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

AATAGCCCACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCAGCT  >  NC_002947/1013617‑1013901
                                                                                                                                                |                                                                                                                                            
aaTAGCCCACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCg                                                                                                                                           >  6:41954/1‑148 (MQ=255)
     cccACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACttt                                                                                                                                     >  5:253280/1‑149 (MQ=255)
      ccACCTTCACGCCATTGTTCGCCTCGTCTGCGCGGTTCTTCGCCTGCAGCACATCCCCTTAGAACCCCACTACCACGCAGATATCGCCGTTGGCCATGTAGCCCGTTTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTtc                                                                                                                                    <  4:4226/149‑1 (MQ=255)
      ccACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACttt                                                                                                                                     >  1:131470/1‑148 (MQ=255)
      ccACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACttt                                                                                                                                     >  7:241456/1‑148 (MQ=255)
      ccACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTtc                                                                                                                                    <  3:198438/149‑1 (MQ=255)
           ttCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATcag                                                                                                                               >  8:14744/1‑149 (MQ=255)
           ttCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATcag                                                                                                                               >  5:232575/1‑149 (MQ=255)
           ttCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATcag                                                                                                                               >  3:242989/1‑149 (MQ=255)
             cACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGAAAGTGCGGCCTTACATcagc                                                                                                                              <  7:25707/148‑1 (MQ=255)
                gCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTtcatc                                                                                                                                 >  4:87434/1‑142 (MQ=255)
                gCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTtcatc                                                                                                                                 <  3:87434/142‑1 (MQ=255)
                 ccATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGAAAGTGCGGCCTTACATcagctattt                                                                                                                         <  8:279889/149‑6 (MQ=38)
                    ttGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCtt                                                                                                                      <  5:12806/149‑1 (MQ=255)
                      gTTCGGCTCGTCGGCGCGTTTCTTCGCCTTCAGCACATCCCCTGAGAACACCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTTTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTg                                                                                                                     <  2:131470/148‑1 (MQ=255)
                      gTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTg                                                                                                                     <  8:241456/148‑1 (MQ=255)
                          gCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCt                                                                                                                 >  6:187297/1‑148 (MQ=255)
                            cTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTtct                                                                                                              >  2:146506/1‑149 (MQ=255)
                             tcgtcgGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTtctc                                                                                                             >  7:66421/1‑149 (MQ=255)
                             tagtcgGCGCGGATGGTCGCCTGCCGCAAATCCCCTGAGCAGCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGGGTAGTTGGACGAGTGGATATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGGCCTTGGCCTtctc                                                                                                             >  6:188651/3‑149 (MQ=255)
                              cggcgGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCGGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGTCCTTCTCg                                                                                                            <  5:188651/146‑1 (MQ=255)
                               gtcgGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCg                                                                                                            <  5:41954/148‑1 (MQ=255)
                                tcgGCTCGTTTCTTCGCCTTCAGCACATCCCCTGAGAACCCCCCTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGGGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTa                                                                                                          <  6:327740/149‑1 (MQ=255)
                                       ggTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTg                                                                                                   >  5:256626/1‑149 (MQ=255)
                                         ttcttcGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTggg                                                                                                 <  5:238216/149‑1 (MQ=255)
                                          tcttcGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCGTTGGCCTTCTCGTAGTCCTTgtgg                                                                                                >  5:286722/1‑149 (MQ=255)
                                            ttcGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTg                                                                                                   <  3:188894/144‑1 (MQ=255)
                                            ttcGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTg                                                                                                   >  4:188894/1‑144 (MQ=255)
                                                   gcagcaCATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCtg                                                                                       >  6:27628/1‑149 (MQ=255)
                                                    cagcaCATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCtgt                                                                                      <  5:187297/149‑1 (MQ=255)
                                                     agcaCATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCtgtg                                                                                     <  4:242989/149‑1 (MQ=255)
                                                     agcaCATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCtgta                                                                                     >  3:180872/1‑148 (MQ=255)
                                                     agcaCATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCtg                                                                                       <  4:180872/147‑1 (MQ=255)
                                                         caTCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATgg                                                                                 <  6:256626/149‑1 (MQ=255)
                                                         caTCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGATGGTCCTTGGCCTTAACGTAGTCCTTGGGGTCTTGGATGTTATgg                                                                                 >  3:234411/1‑149 (MQ=255)
                                                          aTCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATggg                                                                                <  1:146506/149‑1 (MQ=255)
                                                            ccccTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATgggg                                                                               >  8:45341/1‑148 (MQ=255)
                                                             cccTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCTCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTGTTCGCTGTGATGGGGc                                                                              >  4:209925/1‑148 (MQ=255)
                                                             accTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGc                                                                              <  5:318127/147‑1 (MQ=255)
                                                               cTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTgg                                                                                                  >  4:220391/1‑126 (MQ=255)
                                                               cTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTgg                                                                                                  <  3:220391/126‑1 (MQ=255)
                                                                   gaACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGcc                                                                         >  2:241428/1‑147 (MQ=255)
                                                                          ccTAACACGCAGATTTCGGCGTTGGCCAGGGCGCCCGTGGACTTGGACGAGTGGAAATATCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTg                                                                <  4:234411/148‑1 (MQ=255)
                                                                          aCTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTg                                                                <  6:253280/149‑1 (MQ=255)
                                                                                    aGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATc                                                      <  7:45341/149‑1 (MQ=255)
                                                                                      atatCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATCgg                                                    <  8:66421/149‑1 (MQ=255)
                                                                                        atCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATCgg                                                    >  1:321393/1‑147 (MQ=255)
                                                                                           gCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCACGTAGTGCAGGGCGATCGgcagca                                               >  6:246181/1‑149 (MQ=255)
                                                                                                      ggTCCCCCGGGAACTGGTACCCGTGTAAATACCTGATGTACGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGGCTTGGCTGTTATGGGGCAAGCCCAGTTAGTGCAGGGCTATCGCCAGCAGTTCGGGGCCa                                    <  2:86498/149‑1 (MQ=255)
                                                                                                       gTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATCGGCAGCAGTTCGGGGCCAt                                   <  6:232575/149‑1 (MQ=255)
                                                                                                          gCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATCGGCAGCAGTTCgg                                         <  3:41155/140‑1 (MQ=255)
                                                                                                          gCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATCGGCAGCAGTTCgg                                         >  4:41155/1‑140 (MQ=255)
                                                                                                                       gACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGAc                   >  3:62622/1‑149 (MQ=255)
                                                                                                                             tGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCca               >  2:311402/1‑147 (MQ=255)
                                                                                                                               gAAATTGCTTATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCa                             <  2:289807/131‑1 (MQ=255)
                                                                                                                               gAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGc           <  7:14744/149‑1 (MQ=255)
                                                                                                                               gAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGc           >  5:268714/1‑149 (MQ=255)
                                                                                                                                 aaTAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCtt         >  5:14802/1‑149 (MQ=255)
                                                                                                                                  cctgCTGATGAAAGTGCGGCCTTACATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCttt        >  5:74829/4‑149 (MQ=255)
                                                                                                                                     gCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCtttt       >  6:153667/1‑147 (MQ=255)
                                                                                                                                     gCTGATGAAAGTGCGGCCTTACATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTCTCCAGCACAGCGACGCCACAGCtttt       >  1:125985/1‑147 (MQ=255)
                                                                                                                                      cTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGcc                                                                         <  6:141597/80‑1 (MQ=255)
                                                                                                                                      cTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGcc                                                                         <  6:146358/80‑1 (MQ=255)
                                                                                                                                      cTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGcc                                                                         >  5:146358/1‑80 (MQ=255)
                                                                                                                                      cTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGcc                                                                         >  5:141597/1‑80 (MQ=255)
                                                                                                                                        gATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCAGCt  >  2:124936/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                            
AATAGCCCACCTTCACGCCATTGTTCGCCTCGTCGGCGCGGTTCTTCGCCTGCAGCACATCCCCTGAGAACCCCACTACCACGCAGATATCGCCGTTGGCCAGGTCGCCCGTGTACTTGGACGAGTGGAAATAGCTGATGTAAGGCCGCACTTTCATCAGCAGGTCCTTGGCCTTCTCGTAGTCCTTGGGGTCTTGGCTGTGATGGGGCAAGCCCAGGTAGTGCAGGGCGATCGGCAGCAGTTCGGGGCCATTGTCCAGCACCGCGACGCCACAGCTTTTCAGCT  >  NC_002947/1013617‑1013901

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