Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,757,030 G→T 7.1% L29M (CTG→ATG)  PP_3320 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,757,0300GT 7.1% 192.9 / 4.2 65L29M (CTG→ATG) PP_3320hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/33);  new base T (3/2);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-01

TAGGCGTCAGCGTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCA  >  NC_002947/3756883‑3757167
                                                                                                                                                   |                                                                                                                                         
taGGCGTCAGCGTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAg                                                                                                                                           <  2:154705/148‑1 (MQ=255)
 aGGCGTCAGCGTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGCGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGa                                                                                                                                         >  4:292232/1‑149 (MQ=255)
  ggCGTCAGCGTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGa                                                                                                                                         <  1:171870/148‑1 (MQ=255)
  ggCGTCAGCGTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGa                                                                                                                                         <  1:350994/148‑1 (MQ=255)
  ggCGTCAGCGTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGa                                                                                                                                         <  7:348857/148‑1 (MQ=255)
     gTCAGCGTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTc                                                                                                                                     >  3:122220/1‑149 (MQ=255)
         gCGTTGGAAGACAGCACGCTATGCTGACATTTCAGTCTTCGGCACATCCCGTTTGTCATAATGTAACAAGCGACCAGAGTTGCTGAGCAGCAGGAATGTACTGAAGCAGCGCCTTATCGTTCATCCAGGCCAGCCTGAGGACTTCGCCg                                                                                                                                 <  8:328339/149‑1 (MQ=255)
           gTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGAt                                                                                                                               <  5:233113/149‑1 (MQ=255)
               gATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCa                                                                                                                           <  5:36262/149‑1 (MQ=255)
               gATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCACCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTc                                                                                                                            <  1:267379/148‑1 (MQ=255)
                 tGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCAtg                                                                                                                         <  3:292232/149‑1 (MQ=255)
                 tGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCAtg                                                                                                                         <  7:199076/149‑1 (MQ=255)
                        aCGATATGCTGACATTTCCCTCTTCTCCACATCGCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTgg                                                                                                                  <  4:209005/149‑1 (MQ=255)
                         cGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTggg                                                                                                                 >  1:275856/1‑149 (MQ=255)
                              tGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGt                                                                                                                >  7:172276/1‑145 (MQ=255)
                              tGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGt                                                                                                                <  8:172276/145‑1 (MQ=255)
                              tGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGAt                                                                                                             <  5:144692/148‑1 (MQ=255)
                               gCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATgg                                                                                                           >  7:322981/1‑149 (MQ=255)
                                 tGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCg                                                                                                         >  7:277551/1‑149 (MQ=255)
                                   aCATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTgg                                                                                                                  <  2:340726/138‑1 (MQ=255)
                                   aCATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTgg                                                                                                                  >  1:340726/1‑138 (MQ=255)
                                       ttCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGAt                                                                                                    >  2:298957/1‑148 (MQ=255)
                                                 tcCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGc                                                                                                          >  6:365278/1‑132 (MQ=255)
                                                 tcCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGc                                                                                                          <  5:365278/132‑1 (MQ=255)
                                                   cacaTCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGctgctg                                                                                                                   <  2:324140/121‑1 (MQ=255)
                                                   cacaTCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGAGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAggggg                                                                                       >  8:283591/1‑149 (MQ=255)
                                                         ccGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAg                                                                                 <  3:153569/149‑1 (MQ=255)
                                                                   tAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGc                                                                       >  8:36846/1‑149 (MQ=255)
                                                                   tAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGc                                                                       >  2:164698/1‑149 (MQ=255)
                                                                       gTAACAAGCGCCCAGGGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTCGGTGATGGCGGGGATCAGCCCCATGGGGTTACAGGGTATGTGGCGAAg                                                                   >  4:246428/1‑149 (MQ=255)
                                                                       gTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAg                                                                   <  8:192057/149‑1 (MQ=255)
                                                                          aCAAGCGCCCAGAGTTGCTGAGCACCATCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGcac                                                                <  8:322981/149‑1 (MQ=255)
                                                                          aCAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGGGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGcac                                                                <  3:243179/149‑1 (MQ=255)
                                                                            aaGCGCCCAGAGTTGCTGAGCACCAGAAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGaaca                                                               >  4:46490/1‑148 (MQ=255)
                                                                              gcgcCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGACGATGTCATGCTGCTGGGTGATGGCGGGGAGCAGCCCCCGGGGGTTCCAGGGTATGTGGAGAAGCACAGCg                                                            >  2:100596/1‑149 (MQ=255)
                                                                                   cAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTc                                                         >  6:142075/1‑147 (MQ=255)
                                                                                                   cAGCAATGTACTGAAGCAGCGCATTTTCGTTCATCCAGGCAAGAATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCACAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTccc                                       >  6:88813/1‑149 (MQ=255)
                                                                                                     gCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGc                                     <  4:122220/149‑1 (MQ=255)
                                                                                                         tGTACTGAAGGAGCTCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCtgt                                 <  4:160913/149‑1 (MQ=255)
                                                                                                            aCTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGcc                              >  3:109851/1‑149 (MQ=255)
                                                                                                             cTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCt                             <  2:275856/149‑1 (MQ=255)
                                                                                                              tGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGAt                                                                                                                               >  5:30332/1‑50 (MQ=255)
                                                                                                              tGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGAt                                                                                                                               <  6:30332/50‑1 (MQ=255)
                                                                                                              tGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCCTGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTg                            >  5:147212/1‑149 (MQ=255)
                                                                                                              tGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTg                            >  2:77814/1‑149 (MQ=255)
                                                                                                              tGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTAGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCCGTCAAGGCGAAGCGGCTTCCCGCCTTTGCCTg                            >  1:150676/1‑149 (MQ=255)
                                                                                                                aaGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCt                          <  5:142075/149‑1 (MQ=255)
                                                                                                                     gcgcCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAgg                     <  2:5671/149‑1 (MQ=255)
                                                                                                                     gcgcCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAgg                     <  8:277551/149‑1 (MQ=255)
                                                                                                                       gcCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGt                   <  7:36846/149‑1 (MQ=255)
                                                                                                                          ttttCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGAGAACCGGCTTCCCGCCTGTGCCTGAACAAGGGt                   >  7:234526/1‑146 (MQ=255)
                                                                                                                                ttCATCGAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCTGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACAt          <  3:246428/149‑1 (MQ=255)
                                                                                                                                ttCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACAt          >  4:246140/1‑149 (MQ=255)
                                                                                                                                 tCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACAt          >  4:246767/1‑148 (MQ=255)
                                                                                                                                   aTCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAt       >  6:213896/1‑149 (MQ=255)
                                                                                                                                   aTCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAt       >  2:294822/1‑149 (MQ=255)
                                                                                                                                    tCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt      <  4:109851/149‑1 (MQ=255)
                                                                                                                                    tCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt      <  7:283591/149‑1 (MQ=255)
                                                                                                                                    tCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt      <  1:164698/149‑1 (MQ=255)
                                                                                                                                    tCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATCTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGAtt      <  6:152972/149‑1 (MQ=255)
                                                                                                                                      tcacccAGCAGCATGACATCGGCAAAGTAATGATGATGGGTGATGGCGGGGATCATCCCCAAGGGGTTCCAGGGTATGTGGCGAAGCACATCGTTCAAGGCGAAGCGGCTTACAGCATGTTCCTGCTCAAGGGTAAGGAACATGAttt     >  7:56964/5‑148 (MQ=255)
                                                                                                                                        ggCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCa  <  1:298957/149‑1 (MQ=255)
                                                                                                                                         cccAGCAGCATGACATCGGCGAAGTCCTGATGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGc   >  7:275881/2‑147 (MQ=38)
                                                                                                                                         cccAGCAGCATGACATCGGCGAAGTCCTGATGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGc   <  8:275881/146‑1 (MQ=38)
                                                                                                                                             gcagcaTGACATCGGCGAAGTCCTGATGCT‑GGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCt                             >  1:76482/1‑116 (MQ=37)
                                                                                                                                             gcagcaTGACATCGGCGAAGTCCTGATGCT‑GGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCt                             <  2:76482/116‑1 (MQ=37)
                                                                                                                                                   |                                                                                                                                         
TAGGCGTCAGCGTTGGATGCCAGCACGATATGCTGACATTTCCCTCTTCTCCACATCCCGTTTGTCATAACGTAACAAGCGCCCAGAGTTGCTGAGCACCAGCAATGTACTGAAGCAGCGCCTTTTCGTTCATCCAGGCCAGCATCAGGACTTCGCCGATGTCATGCTGCTGGGTGATGGCGGGGATCAGCCCCAGGGGGTTCCAGGGTATGTGGCGAAGCACATCGGTCAAGGCGAAGCGGCTTCCCGCCTGTGCCTGCTCAAGGGTAAGGAACATGATTTGCA  >  NC_002947/3756883‑3757167

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