Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,151,116 A→C 6.4% F117C (TTC→TGC)  PP_4535 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,151,1160AC 6.4% 192.1 / 4.3 63F117C (TTC→TGC) PP_4535hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base C (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

ATGAAGGATGGGAAACTGCGAACAGAGGGCGTTGCACTCCATTTGCCAGTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTGCCCGCCCTACAACAACCGCCATGGTCGCAAAGCGCAGCACGGTG  >  NC_002947/5150972‑5151262
                                                                                                                                                |                                                                                                                                                  
aTGAAGGATGGGAAACTGCGAACAGAGGGCGTTGCACTCCATTTGCCAGTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAAtt                                                                                                                                                <  7:73919/149‑1 (MQ=255)
   aaGGATGGGAAACTGCGAACAGAGGGCGTTGCACTCCATTTGCCAGTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGa                                                                                                                                             <  1:286754/149‑1 (MQ=255)
    aGGATGGGAAACTGCGAACAGAGGGCGTTGCACTCCATTTGCCAGTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGa                                                                                                                                             <  2:134055/148‑1 (MQ=255)
           gAAACTGCGAACAGAGGGCGTTGCACTCCATTTGCCAGTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCaaaa                                                                                                                                      >  8:265321/1‑148 (MQ=255)
            aaaCTGCGAACAGAGGGCGTTGCACTCCATTTGCCAGTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCaaaaaa                                                                                                                                    <  6:92064/149‑1 (MQ=255)
                tGCGAACAGAGGGCGTTGCACTCCATTTGCCAGTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGgtg                                                                                                                                >  3:253841/1‑149 (MQ=255)
                   gAACAGAGGGCGTTGCACTCCATTTGCCAGTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTc                                                                                                                              >  2:58625/1‑148 (MQ=255)
                          gggCGTTGCACTCCATTTGCCAGTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCt                                                                                                                      <  6:269254/149‑1 (MQ=255)
                             cGTTGCACTCCATTTGCCAGTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCaa                                                                                                                   >  8:2120/1‑149 (MQ=255)
                               ttGCACTCCATTTGCCAGTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATc                                                                                                                 <  3:92877/149‑1 (MQ=255)
                                  cACTCCATTTGCCAGTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCaga                                                                                                              <  3:186741/149‑1 (MQ=255)
                                  cACTCCATTTGCCAGTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCaga                                                                                                              <  7:2120/149‑1 (MQ=255)
                                    cTCCATTTGCCAGTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAgg                                                                                                            >  5:118874/1‑149 (MQ=255)
                                             ccAGTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGtt                                                                                                   <  1:58625/149‑1 (MQ=255)
                                                gTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTc                                                                                                <  7:323274/149‑1 (MQ=255)
                                                gTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTc                                                                                                <  7:265321/149‑1 (MQ=255)
                                                   aTCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAAtt                                                                                                                                                >  4:52079/1‑98 (MQ=255)
                                                              cAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGgc                                                                                                           >  1:357631/1‑124 (MQ=255)
                                                              cAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGgc                                                                                                           <  2:357631/124‑1 (MQ=255)
                                                               agagGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCgg                                                                                  >  5:350067/1‑148 (MQ=255)
                                                                   gATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCCCCGCAGGCCCGGAAAcc                                                                             >  6:47305/1‑149 (MQ=255)
                                                                        gtgtCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACg                                                                        >  1:136149/1‑149 (MQ=255)
                                                                          gtCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCt                                                                                                                      <  6:363964/101‑1 (MQ=255)
                                                                          gtCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCt                                                                                                                      >  5:363964/1‑101 (MQ=255)
                                                                           tCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGGTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCAc                                                                         <  2:161622/145‑1 (MQ=255)
                                                                           tCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCACAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACgag                                                                      >  1:161622/1‑146 (MQ=255)
                                                                             gATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAgg                                                                    >  5:59867/1‑148 (MQ=255)
                                                                              aTGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAgg                                                                    >  1:73273/1‑147 (MQ=255)
                                                                               tGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAg                                                                 >  4:141150/1‑149 (MQ=255)
                                                                                   ctcCATGTTGACACCAAAGTAACGACCCAGGATCTTCAGAAATTCGGCCTTGTGCTGGAAGATATTCGAGGCAAAAAATGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGATACCTCACGTAGGAAGCAc                                                              <  6:132085/148‑1 (MQ=255)
                                                                                      cAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGAt                                                          >  8:28810/1‑149 (MQ=255)
                                                                                          gTGACACCAAAGTAACGACCGAGGCTCTTCAGAAAATCGGCATTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCt                                                      <  4:253841/149‑1 (MQ=255)
                                                                                              cacCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGcc                                                  >  1:213461/1‑149 (MQ=255)
                                                                                              cacCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGc                                                   >  5:371464/1‑148 (MQ=255)
                                                                                              cacCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGc                                                   <  6:371464/148‑1 (MQ=255)
                                                                                               acCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGTGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCg                                                 >  1:41935/1‑149 (MQ=255)
                                                                                               acCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGCTTGCTGGCCg                                                 >  4:338070/1‑149 (MQ=255)
                                                                                               acCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCg                                                 >  5:22851/1‑149 (MQ=255)
                                                                                               acCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCg                                                 <  1:6765/149‑1 (MQ=255)
                                                                                                 cAATGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGtt                                                                                                   <  3:300277/97‑1 (MQ=255)
                                                                                                 cAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGtt                                                                                                   >  4:300277/1‑97 (MQ=255)
                                                                                                      tAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAATAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCAAGTAGGAAGCACCGATTGCTGGCCGCTgcgcg                                          >  6:49157/1‑149 (MQ=255)
                                                                                                       aaCGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATc                                                                                                                 <  8:226806/77‑1 (MQ=255)
                                                                                                       aaCGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATc                                                                                                                 >  7:226806/1‑77 (MQ=255)
                                                                                                       aaCGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGgcagc                                                                                                        <  4:155029/86‑1 (MQ=255)
                                                                                                            cccAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTGCCCGCCCTa                                     >  5:308386/1‑148 (MQ=255)
                                                                                                             ccAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTgcccgccc                                       >  6:58525/1‑145 (MQ=255)
                                                                                                             ccAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTgcccgccc                                       <  5:58525/145‑1 (MQ=255)
                                                                                                                   tctTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTgccc                                           <  8:315322/135‑1 (MQ=255)
                                                                                                                   tctTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTgccc                                           >  7:315322/1‑135 (MQ=255)
                                                                                                                      tCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCAGCTGCCCGCCCTACAACAACCGc                           >  3:281095/1‑148 (MQ=255)
                                                                                                                        aGAAAATCGGCCTTGGGCTGGAAGCAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCaga                                                                                                              <  2:331308/63‑1 (MQ=255)
                                                                                                                        aGAAAATCGGCCTTGGGCTGGAAGCAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCaga                                                                                                              <  2:325942/63‑1 (MQ=255)
                                                                                                                        aGAAAATCGGCCTTGGGCTGGAAGCAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCaga                                                                                                              >  1:325942/1‑63 (MQ=255)
                                                                                                                        aGAAAATCGGCCTTGGGCTGGAAGCAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCaga                                                                                                              >  1:331308/1‑63 (MQ=255)
                                                                                                                        aGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTGCCCGCCCTacaa                                  <  3:266392/139‑1 (MQ=255)
                                                                                                                        aGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTGCCCGCCCTacaa                                  >  4:266392/1‑139 (MQ=255)
                                                                                                                         gAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTGCCCGCCCTACAACAACCGCCATg                       >  6:281532/1‑149 (MQ=255)
                                                                                                                           aaaTCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTGCCCGCCCTACAACAACCGCCATgg                      <  6:59867/148‑1 (MQ=255)
                                                                                                                            aaTCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTGCCCGCCCTACAACAACCGCCATGGTc                    <  2:73273/149‑1 (MQ=255)
                                                                                                                                   cTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTGCCCGCCCTACAACAACCGCCATGGTCGCAAAgc             <  2:213461/149‑1 (MQ=255)
                                                                                                                                   cTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTGCCCGCCCTACAACAACCGCCATGGTCGCAAAgc             <  6:22851/149‑1 (MQ=255)
                                                                                                                                          tGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTGCCCGCCCTACAACAACCGCCATGGTCGCAAAGCGCAGCAc      >  3:214041/1‑149 (MQ=255)
                                                                                                                                            gaagaaATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTGCCCGCCCTACAACAACCGCCATGGTCGCAAAGCGCAGCACgg    >  2:170955/1‑149 (MQ=255)
                                                                                                                                             aagaaATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTGCCCGCCCTACAACAACCGCCATGGTCGCAAAGCGCAGCACGGt   <  7:28810/149‑1 (MQ=255)
                                                                                                                                               gaaATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTGCCCGCCCTACAACAACCGCCATGGTCGCAAAGCGCAGCACGGTg  >  1:329515/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
ATGAAGGATGGGAAACTGCGAACAGAGGGCGTTGCACTCCATTTGCCAGTCATCTCCTGCTGCAGAGGATCGGTGTCGATGTTCTCCAGGGTGACACCAAAGTAACGACCCAGGCTCTTCAGAAAATCGGCCTTGGGCTGGAAGAAATTCGAGGCAAAAAAGGTGTCAACATCCTCAATCAGAGGCAGCACGTTGTCCACCGCAGGCCCGGAAACCTCACGTAGGAAGCACCGATTGCTGGCCGCTGCCCGCCCTACAACAACCGCCATGGTCGCAAAGCGCAGCACGGTG  >  NC_002947/5150972‑5151262

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