Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,356,810 T→C 5.8% K651E (AAG→GAG)  infB ← translation initiation factor IF‑2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,356,8100TC 5.8% 202.7 / 2.3 69K651E (AAG→GAG) infBtranslation initiation factor IF‑2
Reads supporting (aligned to +/- strand):  ref base T (34/31);  new base C (2/2);  total (36/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 = 8.92e-01

AGTACTGCATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTTTCTTCTCGTC  >  NC_002947/5356662‑5356954
                                                                                                                                                    |                                                                                                                                                
aGTACTGCATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCtt                                                                                                                                                  <  6:143164/149‑1 (MQ=255)
aGTACTGCATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCtt                                                                                                                                                  <  3:164977/149‑1 (MQ=255)
  tACTGCATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTga                                                                                                                                                <  6:89407/149‑1 (MQ=255)
  tACTGCATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTga                                                                                                                                                <  1:233517/149‑1 (MQ=255)
     tGCATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGaa                                                                                                                                             >  6:59901/1‑149 (MQ=255)
     tGCATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGaa                                                                                                                                             >  2:5320/1‑149 (MQ=255)
     tGCATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGaa                                                                                                                                             >  2:22198/1‑149 (MQ=255)
            gAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGtt                                                                                                                                      >  8:164652/1‑149 (MQ=255)
            gAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGtt                                                                                                                                      >  1:242541/1‑149 (MQ=255)
             aaGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTg                                                                                                                                     <  1:106376/149‑1 (MQ=255)
              aGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGtt                                                                                                                                      >  8:132549/1‑147 (MQ=255)
               gCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGTTCTCAAGACCGATGtgcgt                                                                                                                                   <  8:79271/149‑5 (MQ=255)
               gCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAg                                                                                                                                   <  5:131653/149‑1 (MQ=255)
                  aGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAg                                                                                                                                   <  3:249891/146‑1 (MQ=255)
                  aGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAg                                                                                                                                   >  4:249891/1‑146 (MQ=255)
                   gcgcCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGtt                                                                                                                               >  5:295165/1‑149 (MQ=255)
                    cgcCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACGAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGttt                                                                                                                              >  2:223497/1‑149 (MQ=255)
                      ccAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCTCTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTATTGAGAACGATGTTGAGGGTTttc                                                                                                                            <  6:208376/149‑1 (MQ=255)
                               cATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAgg                                                                                                                                  <  2:50199/134‑1 (MQ=255)
                               cATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAgg                                                                                                                                  >  1:50199/1‑134 (MQ=255)
                               cATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTtcttct                                                                                                                        <  6:137556/144‑1 (MQ=255)
                               cATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTtcttct                                                                                                                        >  5:137556/1‑144 (MQ=255)
                               cATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTcc                                                                                                                   >  3:112358/1‑149 (MQ=255)
                                    ctctCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGAcc                                                                                                              >  1:172332/1‑149 (MQ=255)
                                     tctcGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGAcc                                                                                                              >  3:156558/1‑148 (MQ=255)
                                       tcGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCAt                                                                                                           >  5:118966/1‑149 (MQ=255)
                                         ggTTATACCAACGACGCCACCACCAATCACGCGCACTTGCAATTAGTCGTTGCCCACGCCGCCGACCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATg                                                                                                          <  4:13842/148‑1 (MQ=255)
                                         ggTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATg                                                                                                          <  2:172332/148‑1 (MQ=255)
                                           tGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGGCCCATGGtc                                                                                                       >  2:22235/1‑149 (MQ=255)
                                             aTACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGtct                                                                                                      <  7:132549/148‑1 (MQ=255)
                                             aTACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGtct                                                                                                      <  4:112358/148‑1 (MQ=255)
                                                 caccGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCg                                                                                                    >  7:250961/1‑146 (MQ=255)
                                                   ccGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGGCCCATGGTCTCGAAGAt                                                                                               <  1:22235/149‑1 (MQ=255)
                                                      aCGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGtt                                                                                            <  6:118966/149‑1 (MQ=255)
                                                        gCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCg                                                                                                    >  6:13536/1‑139 (MQ=255)
                                                        gCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCg                                                                                                    <  5:13536/139‑1 (MQ=255)
                                                          caccaccaATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTc                                                                                         >  3:95915/1‑148 (MQ=255)
                                                           accaccaATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTc                                                                                         >  7:259197/1‑147 (MQ=255)
                                                           accaccaATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCa                                                                                       <  5:59901/149‑1 (MQ=255)
                                                            ccaccaATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTttct                                                                                                                           <  4:243316/112‑1 (MQ=255)
                                                            ccaccaATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATg                                                                                                          >  7:231448/1‑129 (MQ=255)
                                                            ccaccaATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATg                                                                                                          <  8:231448/129‑1 (MQ=255)
                                                               ccaATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCtt                                                                                   <  2:242541/149‑1 (MQ=255)
                                                                     aCGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCgg                                                                              >  7:176429/1‑148 (MQ=255)
                                                                         gcacttgcacttCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGa                                                                         <  6:295165/149‑1 (MQ=255)
                                                                           acttgcacttCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAg                                                                        >  5:321316/1‑148 (MQ=255)
                                                                                 acttCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGc                                                                 >  4:73767/1‑149 (MQ=255)
                                                                                   ttCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGACGTGAGCACGGGc                                                                 >  5:131370/1‑147 (MQ=255)
                                                                                     cgtcgtTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGc                                                             >  8:238592/1‑149 (MQ=255)
                                                                                          ttGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGAc                                                        <  1:223497/149‑1 (MQ=255)
                                                                                               cAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTcgc                                                   <  4:95915/149‑1 (MQ=255)
                                                                                                 ggCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGtt                                                                                            <  8:293314/106‑1 (MQ=255)
                                                                                                 ggCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGtt                                                                                            >  7:293314/1‑106 (MQ=255)
                                                                                                          gCGAACCCTGCAGCGCTTCCAGACAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGATTGCCGCCGTGAGCACGGGCCAGCTTGACATCGCGATACTTGCCt                                        >  4:158723/1‑149 (MQ=255)
                                                                                                             aaCCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCATGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGa                                     <  8:176429/149‑1 (MQ=255)
                                                                                                                 cTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGAAgg                                  >  7:294956/1‑148 (MQ=255)
                                                                                                                     aGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAg                             <  7:164652/149‑1 (MQ=255)
                                                                                                                        gcTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGc                          <  1:106651/149‑1 (MQ=255)
                                                                                                                            ccAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATATTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACt                      >  6:103052/1‑149 (MQ=255)
                                                                                                                               gagaACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGcc                                         >  6:127208/1‑127 (MQ=255)
                                                                                                                               gagaACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGcc                                         <  5:127208/127‑1 (MQ=255)
                                                                                                                                agaACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCAc                  >  2:6111/1‑149 (MQ=255)
                                                                                                                                agaACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCAc                  <  6:321316/149‑1 (MQ=255)
                                                                                                                                       tcaacaTCGTTCTCAAGACCGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGa                                     <  2:218788/120‑1 (MQ=255)
                                                                                                                                       tcaacaTCGTTCTCAAGACCGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGa                                     >  1:218788/4‑123 (MQ=255)
                                                                                                                                          acaTCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTttctt        <  1:6111/149‑1 (MQ=255)
                                                                                                                                          acaTCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTttctt        <  7:238592/149‑1 (MQ=255)
                                                                                                                                           caTCGTTCTCAAGACCGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTttctt        >  3:55650/1‑148 (MQ=255)
                                                                                                                                           caTCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTttcttc       >  7:228703/1‑149 (MQ=255)
                                                                                                                                           caTCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTttcttc       <  1:22198/149‑1 (MQ=255)
                                                                                                                                                gTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAAATTCACGGGCTTTCTTCtcgtc  >  6:182978/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AGTACTGCATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTTTCTTCTCGTC  >  NC_002947/5356662‑5356954

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