Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,146,831 A→T 6.8% L411H (CTC→CAC)  PP_1006 ← heme receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,146,8310AT 6.8% 215.8 / 5.1 74L411H (CTC→CAC) PP_1006heme receptor
Reads supporting (aligned to +/- strand):  ref base A (30/39);  new base T (2/3);  total (32/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.94e-01

AATACCCACTGGCGATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGTCACGGTCATCC  >  NC_002947/1146684‑1146971
                                                                                                                                                    |                                                                                                                                            
aagaCCCACTGGCGATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAg                                                                                                                                             <  5:188851/146‑1 (MQ=255)
 aTACCCACTGGCGATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAgg                                                                                                                                            >  2:43405/1‑149 (MQ=255)
  tACCCACTGGCGATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGt                                                                                                                                           >  7:70098/1‑149 (MQ=255)
      cACTGGCTATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCTCGATAGCGTGGTCgcg                                                                                                                                       <  1:38588/149‑3 (MQ=255)
       aCTGGCGATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAAt                                                                                                                                      <  6:84866/149‑1 (MQ=255)
         tGGCGATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAg                                                                                                                                             >  5:206812/1‑140 (MQ=255)
         tGGCGATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAg                                                                                                                                             <  6:206812/140‑1 (MQ=255)
             gATCCTTGTAGGTAGTCTGCCGTTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGa                                                                                                                                <  8:248260/149‑1 (MQ=255)
             gATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGa                                                                                                                                >  2:150339/1‑149 (MQ=255)
             gATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGa                                                                                                                                <  3:251017/149‑1 (MQ=255)
               tCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTTGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGAtt                                                                                                                              <  8:6838/149‑1 (MQ=255)
                 cTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGa                                                                                                                            <  3:143984/149‑1 (MQ=255)
                 cTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCACGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTGATGAtt                                                                                                                              >  1:284114/1‑147 (MQ=255)
                  ttGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGTGACCACGCTATCGAGGTCGAATTCATGAtt                                                                                                                              >  6:36658/1‑146 (MQ=255)
                  ttGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAg                                                                                                                           >  6:368441/1‑149 (MQ=255)
                    gtaggtagTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGTGACCACGCTATCGAGGTCGAATTCATGATTGa                                                                                                                            <  5:36658/146‑1 (MQ=255)
                        gtagTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGaa                                                                                                                     <  1:43405/149‑1 (MQ=255)
                            tCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAg                                                                                                                           <  6:274761/139‑1 (MQ=255)
                            tCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAg                                                                                                                           >  5:274761/1‑139 (MQ=255)
                            tCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGt                                                                                                                 >  8:322312/1‑149 (MQ=255)
                              tGCCGGTCACGAAACACTTGGCGTCCCTAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGtt                                                                                                                >  7:106422/1‑148 (MQ=255)
                               gCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGt                                                                                                                 <  7:116831/146‑1 (MQ=255)
                               gCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTAATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGt                                                                                                                 >  8:116831/1‑146 (MQ=255)
                                   gTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACgg                                                                                                           >  5:189480/1‑148 (MQ=255)
                                   gTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACgg                                                                                                           <  6:189480/148‑1 (MQ=255)
                                     cccGAACCACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTTGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTg                                                                                                        <  2:284114/149‑1 (MQ=255)
                                     cccGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGt                                                                                                         >  6:355795/1‑148 (MQ=255)
                                        gAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATc                                                                                                     >  4:314137/1‑149 (MQ=255)
                                        gAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATc                                                                                                     <  6:178176/149‑1 (MQ=255)
                                        gAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GCGGGGCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATc                                                                                                     <  1:302798/149‑1 (MQ=255)
                                         aaaaACATGGCGTCACGAGGCTACAAACACCTCAAACGTAC‑CTTCGTCCGTCTACCCGATCAATTAGTTCAAACACCATTTGACGTGATCCGAGACCACGCTATCGAGGTCGAATTCAAGATTCACTCACAAGCTTTCACgggagctc                                                                                                     >  4:286669/1‑143 (MQ=255)
                                          aaCACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCgt                                                                                                   <  3:236961/149‑1 (MQ=255)
                                           acacTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCgtg                                                                                                  <  5:368441/149‑1 (MQ=255)
                                               ttGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGAtt                                                                                                                              >  3:47826/1‑117 (MQ=255)
                                               ttGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGAtt                                                                                                                              <  4:47826/117‑1 (MQ=255)
                                                  gCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGc                                                                                           <  1:150339/149‑1 (MQ=255)
                                                     tCCCGAGGCGACATACAGCTCATCGGTACGGTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCgtgt                                                                                                 >  1:215893/1‑141 (MQ=255)
                                                     tCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCgtgt                                                                                                 <  2:215893/140‑1 (MQ=255)
                                                     tCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCtt                                                                                        <  2:231290/149‑1 (MQ=255)
                                                       ccGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTgg                                                                                      >  3:201245/1‑149 (MQ=255)
                                                            gCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCTCGATAGCGTGGTCgcgg                                                                                                                                      >  7:234792/1‑93 (MQ=255)
                                                            gCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCTCGATAGCGTGGTCgcgg                                                                                                                                      <  8:234792/96‑4 (MQ=255)
                                                                  tACAGCTCATCGGTAA‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCTAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTa                                                                           <  5:219921/149‑1 (MQ=255)
                                                                    cAGCTCATCGGTAC‑GTTGGTCGGTCTTCGAGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTAc                                                                          <  3:314137/148‑1 (MQ=255)
                                                                       cTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTg                                                                                       <  2:150035/132‑1 (MQ=255)
                                                                       cTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTg                                                                                       >  1:150035/1‑132 (MQ=255)
                                                                          aTCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCa                                                                   <  7:322312/149‑1 (MQ=255)
                                                                           tCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCAt                                                                  <  8:70098/149‑1 (MQ=255)
                                                                            cGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCAt                                                                  <  8:106422/148‑1 (MQ=255)
                                                                             ggTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGc                                                                <  5:355795/149‑1 (MQ=255)
                                                                              gTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCTCGATAGCGTGGTCgc                                                                                                                                        >  5:11986/1‑75 (MQ=38)
                                                                              gTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCTCGATAGCGTGGTCgc                                                                                                                                        <  6:11986/76‑2 (MQ=38)
                                                                                       gtcggtcTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGc                                                       >  1:83369/1‑149 (MQ=255)
                                                                                        tcggtcTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGcc                                                      >  4:64506/1‑149 (MQ=255)
                                                                                          ggtcTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCg                                                    <  2:83369/149‑1 (MQ=255)
                                                                                            tcTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCgg                                                   >  5:215729/1‑148 (MQ=255)
                                                                                              ttCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATg                                                                                >  8:329557/1‑117 (MQ=255)
                                                                                              ttCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATg                                                                                <  7:329557/117‑1 (MQ=255)
                                                                                              ttCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGa                                                  >  4:151799/1‑147 (MQ=255)
                                                                                              ttCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGa                                                  <  3:151799/147‑1 (MQ=255)
                                                                                                 gcgATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGaa                                             <  2:180237/149‑1 (MQ=255)
                                                                                                        ggTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGGCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGccca                                      >  2:317522/1‑149 (MQ=255)
                                                                                                         gTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGccca                                      <  3:64506/148‑1 (MQ=255)
                                                                                                                   cTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCt                           >  5:155151/1‑149 (MQ=255)
                                                                                                                            aCGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGtt                  <  3:286669/149‑1 (MQ=255)
                                                                                                                                 aTCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGt             >  4:49756/1‑149 (MQ=255)
                                                                                                                                   ccGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGTCa           >  3:359135/1‑149 (MQ=255)
                                                                                                                                    gcggatcACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTc                 <  5:299591/136‑1 (MQ=255)
                                                                                                                                    gcggatcACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTc                 >  6:299591/7‑142 (MQ=255)
                                                                                                                                     gcgACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGTCACg         >  2:149985/1‑149 (MQ=255)
                                                                                                                                     gcgACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGTCACg         <  4:359135/149‑1 (MQ=255)
                                                                                                                                       gACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGGCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGTCACGGt       <  1:317522/149‑1 (MQ=255)
                                                                                                                                           aCGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGTCACGGTCATc   <  1:149985/149‑1 (MQ=255)
                                                                                                                                            cGCTATCGAGGGCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGTCACGGTCATcc  <  4:201245/149‑1 (MQ=255)
                                                                                                                                                    aGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTc                 <  3:266226/126‑1 (MQ=255)
                                                                                                                                                    aGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTc                 >  4:266226/1‑126 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
AATACCCACTGGCGATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGTCACGGTCATCC  >  NC_002947/1146684‑1146971

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