Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,784,316 C→T 6.4% Y242Y (TAC→TAT PP_2437 → acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,784,3160CT 6.4% 180.9 / 2.5 63Y242Y (TAC→TATPP_2437acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (34/24);  new base T (2/2);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

GGGCCACGCCACCAGCAGCATTCCATGGTGCTGGTGCCGACCGACACGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCAAGCGTTCGGTGGAGCGTACTGCTTT  >  NC_002947/2784170‑2784465
                                                                                                                                                  |                                                                                                                                                     
gggCCACGCCACCAGCAGCATTCCATGGTGCTGGTGCCGACCGACACGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGaa                                                                                                                                                    >  1:153522/1‑150 (MQ=255)
     acgccacCAGCAGCATTCCATGGTGCTGGTGCCGACCGACACGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAAtgt                                                                                                                                               >  1:36488/1‑150 (MQ=255)
     acgccacCAGCAGCATTCCATGGTGCTGGTGCCGACCGACACGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAAtgt                                                                                                                                               >  1:414259/1‑150 (MQ=255)
         caccaGCAGCATTCCATGGTGCTGGTGCCGACCGACACGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATc                                                                                                                                           >  2:1127010/1‑150 (MQ=255)
           ccaGCAGCATTCCATGGTGCTGGTGCCGACCGACACGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAACGTGCGGGTGCCATACGAAAATGTGATCCt                                                                                                                                         >  1:861955/1‑150 (MQ=255)
             agcagcATTCCATGGTGCTGGTGCCGACCGACACGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTg                                                                                                                                       <  2:240232/150‑1 (MQ=255)
                  caTTCCATGGTGCTGGTGCCGACCGACACGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTaaa                                                                                                                                  <  1:1281787/150‑1 (MQ=255)
                  caTTCCATGGTGCTGGTGCCGACCGACACGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGaa                                                                                                                                  <  2:411872/150‑1 (MQ=255)
                    ttCCATGGTGCTGGTGCCGACCGACACGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGGGGGACGAAAATGTGATCCTTGGTGAAgg                                                                                                                                >  1:965470/1‑150 (MQ=255)
                    ttCCATGGTGCTGGTGCCGACCGACACGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAgg                                                                                                                                >  2:770479/1‑150 (MQ=255)
                          ggtgctgttgcCGACCGACACGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCgg                                                                                                                          >  2:1282631/1‑150 (MQ=255)
                                  tgcCGACCGACACGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGaaa                                                                                                                  <  2:388103/150‑1 (MQ=255)
                                      gaccgacACGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGc                                                                                                              >  1:854445/1‑150 (MQ=255)
                                       accgacACGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACACGGCTTCGAAATTGcc                                                                                                             >  1:873351/1‑150 (MQ=255)
                                          gacACGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAg                                                                                                          >  2:1206403/1‑150 (MQ=255)
                                           acacGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATATGAAAATGTGATc                                                                                                                                           <  1:382462/116‑1 (MQ=255)
                                           acacGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATATGAAAATGTGATc                                                                                                                                           <  1:382899/116‑1 (MQ=255)
                                           acacGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATATGAAAATGTGATc                                                                                                                                           >  2:382462/1‑116 (MQ=255)
                                           acacGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATATGAAAATGTGATc                                                                                                                                           >  2:382899/1‑116 (MQ=255)
                                           acacGCCGGGCGTGAAGAGCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAgg                                                                                                         >  2:399831/1‑150 (MQ=255)
                                                                     ccgcTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATccac                                                                               <  1:1282631/150‑1 (MQ=255)
                                                                         tGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACt                                                                           <  2:153522/150‑1 (MQ=255)
                                                                                  ccGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCGGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCt                                                                  >  2:879467/2‑150 (MQ=255)
                                                                                    ggTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCAGTATCCACTACGCCATTCGGTCg                                                                >  2:997377/1‑150 (MQ=255)
                                                                                      ttATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCCTGCGCTCGAt                                                              >  2:153308/1‑150 (MQ=255)
                                                                                      ttATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGAt                                                              <  1:879467/150‑1 (MQ=255)
                                                                                        aTGACGATGCGCCCCATGGTCATGCCGAAGTCATTATCGAGACAGAACGAGGGCCAAAAGATAATGAGATCCCGGTTGAACTACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCg                                                            >  2:1109966/1‑150 (MQ=255)
                                                                                          gACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTgg                                                                                                <  1:867972/112‑1 (MQ=255)
                                                                                          gACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTgg                                                                                                >  2:867972/1‑112 (MQ=255)
                                                                                             gATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGc                                                                         >  2:604716/1‑132 (MQ=255)
                                                                                             gATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGc                                                                         <  1:604716/132‑1 (MQ=255)
                                                                                             gATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATg                                                       <  1:38419/150‑1 (MQ=255)
                                                                                              aTGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATgg                                                      <  2:854445/150‑1 (MQ=255)
                                                                                                     ccATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAgcg                                               >  1:195454/1‑150 (MQ=255)
                                                                                                     ccATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAgc                                                >  2:1153170/1‑149 (MQ=255)
                                                                                                     ccATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAgc                                                <  1:1153170/149‑1 (MQ=255)
                                                                                                            cATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGACATCTCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTc                                        <  1:997377/150‑1 (MQ=255)
                                                                                                              tGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGa                                      <  1:770479/150‑1 (MQ=255)
                                                                                                                    aGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGAt                                >  1:585061/1‑150 (MQ=255)
                                                                                                                    aGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGAt                                <  2:195454/150‑1 (MQ=255)
                                                                                                                      tCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGAtgt                              <  2:429013/150‑1 (MQ=255)
                                                                                                                       ccgtttCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAg                                     >  1:1284020/4‑142 (MQ=255)
                                                                                                                       ccTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAg                                     <  2:1284020/142‑1 (MQ=255)
                                                                                                                        cTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGGAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGc                            >  2:548745/1‑150 (MQ=255)
                                                                                                                        cTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGc                            >  1:839950/1‑150 (MQ=255)
                                                                                                                        cTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGc                            <  1:1127010/150‑1 (MQ=255)
                                                                                                                         ttttCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGAtgt                              <  2:202980/147‑1 (MQ=255)
                                                                                                                         ttttCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGAtgt                              >  1:202980/1‑147 (MQ=255)
                                                                                                                         ttttCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCa                           >  1:327123/1‑150 (MQ=255)
                                                                                                                          tttCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCaa                          <  1:1206403/150‑1 (MQ=255)
                                                                                                                              gagaATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCAAGCGt                      >  2:538857/1‑150 (MQ=255)
                                                                                                                              gagaATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCAAGCGt                      >  1:627097/1‑150 (MQ=255)
                                                                                                                                  aTGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCAAGCGTTCgg                  <  2:627097/150‑1 (MQ=255)
                                                                                                                                    gtgCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCAAGCGTTCGGTg                <  1:153308/150‑1 (MQ=255)
                                                                                                                                      gCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCAAGCGTTCGGTGGa              >  2:1106889/1‑150 (MQ=255)
                                                                                                                                          gTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCAAGCGTTCGGTGGAGCGt          >  2:293262/1‑150 (MQ=255)
                                                                                                                                           tGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCAAGCGTTCGGTGGAGCGTa         <  1:538857/150‑1 (MQ=255)
                                                                                                                                              cATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCAAGCGTTCGGTGGAGCGTACTg      >  2:581480/1‑150 (MQ=255)
                                                                                                                                              cATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCAAGCGTTCGGTGGAGCGTACTg      >  1:454671/1‑150 (MQ=255)
                                                                                                                                                tACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCAAGCGTTCGGTGGAGCGTACTGCt    >  2:5134/1‑150 (MQ=255)
                                                                                                                                                tACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCAAGCGTTCGGTGGAGCGTACTGCt    >  1:255577/1‑150 (MQ=255)
                                                                                                                                                tACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCAAGCGTTCGGTGGAGCGTACTGCt    >  2:805502/1‑150 (MQ=255)
                                                                                                                                                 aCGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCAAGCGTTCGGTGGAGCGTACTGCtt   <  2:36488/150‑1 (MQ=255)
                                                                                                                                                  cGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCAAGCGTTCGGTGGAGCGTACTGCttt  >  2:376555/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGGCCACGCCACCAGCAGCATTCCATGGTGCTGGTGCCGACCGACACGCCGGGCGTGAAGATCGTCCGCCCGCTGCCGGTGTTCGGTTATGACGATGCGCCCCATGGTCATGCCGAAGTCCTTTTCGAGAATGTGCGGGTGCCATACGAAAATGTGATCCTTGGTGAAGGACGCGGCTTCGAAATTGCCCAGGGCCGCCTGGGCCCAGGCCGTATCCACCACTGCATGCGCTCGATCGGCATGGCCGAGCGCGCACTCGAGTTGATGTGCAAGCGTTCGGTGGAGCGTACTGCTTT  >  NC_002947/2784170‑2784465

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