Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 68,598 T→G 5.1% P225P (CCA→CCC PP_0058 ← 1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294768,5980TG 5.1% 1145.2 / 34.3 335P225P (CCA→CCCPP_00581‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (164/152);  new base G (9/8);  total (175/160)
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.35e-01

Only 100 of 335 total aligned reads displayed.

ATCCAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCGGTAC  >  NC_002947/68453‑68746
                                                                                                                                                 |                                                                                                                                                    
  ccAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGcccc                                                                                                                                                >  2:6589429/1‑150 (MQ=255)
    aaaaaCTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGt                                                                                                                                              <  2:2029960/150‑1 (MQ=255)
    aaaaaCTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGt                                                                                                                                              <  2:1848904/150‑1 (MQ=255)
    aaaaaCTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGt                                                                                                                                              >  2:2408699/1‑150 (MQ=255)
    aaaaaCTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGt                                                                                                                                              >  1:746541/1‑150 (MQ=255)
           gCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTg                                                                                                                                       <  2:6024518/150‑1 (MQ=255)
           gCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTg                                                                                                                                       <  2:6143534/150‑1 (MQ=255)
            cacccaAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAACGGCACGGGGCCACGTGCCATTgc                                                                                                                                      <  1:4556373/150‑2 (MQ=255)
            cacacaAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTgg                                                                                                                                      <  2:858010/150‑1 (MQ=255)
             acacaAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGc                                                                                                                                     <  1:265337/150‑1 (MQ=255)
             acacaAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGc                                                                                                                                     >  1:3072201/1‑150 (MQ=255)
             acacaAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGCCAACGGCACGGGGCCACGTGCCATTgcc                                                                                                                                     <  2:708371/150‑3 (MQ=255)
               acaAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGt                                                                                                                                   >  2:4936347/1‑150 (MQ=255)
               acaAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTGGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGt                                                                                                                                   >  2:2465110/1‑150 (MQ=255)
                 aaGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTAc                                                                                                                                 >  2:5401504/1‑150 (MQ=255)
                      aTCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCgg                                                                                                                            <  1:5508356/150‑1 (MQ=255)
                       tCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGc                                                                                                                           >  1:4495508/1‑150 (MQ=255)
                          acaGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCa                                                                                                                        <  1:385682/150‑1 (MQ=255)
                          acaGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCa                                                                                                                        >  1:358097/1‑150 (MQ=255)
                              aaGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCaa                                                                                                                    <  1:4626774/150‑1 (MQ=255)
                               aGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAAt                                                                                                                   <  2:4746833/150‑1 (MQ=255)
                               aGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAAt                                                                                                                   <  1:2129130/150‑1 (MQ=255)
                                ggTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAAtt                                                                                                                  <  2:746337/150‑1 (MQ=255)
                                ggTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAAtt                                                                                                                  <  1:4377876/150‑1 (MQ=255)
                                  tCAAGCCATCTGCTCCTGAGGATTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTAc                                                                                                                <  1:241495/150‑1 (MQ=255)
                                    aaGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCa                                                                                                              >  2:2714370/1‑150 (MQ=255)
                                    aaGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCa                                                                                                              >  2:6118735/1‑150 (MQ=255)
                                     aGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCAc                                                                                                             <  2:314140/150‑1 (MQ=255)
                                        cATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACATc                                                                                                          >  1:6062480/1‑150 (MQ=255)
                                          tCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGa                                                                                                        <  2:2107361/150‑1 (MQ=255)
                                           cTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGAt                                                                                                       <  1:919219/150‑1 (MQ=255)
                                            tGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCACTTACCACTTCGAtt                                                                                                      <  2:3868153/150‑1 (MQ=255)
                                            tGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGAtt                                                                                                      <  2:5507926/150‑1 (MQ=255)
                                            tGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGAtt                                                                                                      <  1:3284293/150‑1 (MQ=255)
                                               tCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTa                                                                                                   >  1:4699818/1‑150 (MQ=255)
                                                 cTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTAcc                                                                                                 >  1:5199428/1‑150 (MQ=255)
                                                  tGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCg                                                                                                >  2:5722282/1‑150 (MQ=255)
                                                  tGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCg                                                                                                >  2:3299739/1‑150 (MQ=255)
                                                  tGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCg                                                                                                >  1:1754571/1‑150 (MQ=255)
                                                  tGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACAACGGCTCACCAATTACCACTTCGATTGTAACg                                                                                                >  1:5749737/1‑150 (MQ=255)
                                                   gAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCgg                                                                                               >  1:3716223/1‑150 (MQ=255)
                                                           tCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTg                                                                                       >  2:1412907/1‑150 (MQ=255)
                                                                  cTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAAc                                                                                <  2:789218/150‑1 (MQ=255)
                                                                  cTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAAc                                                                                <  1:3524620/150‑1 (MQ=255)
                                                                   tGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCATTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAAcc                                                                               <  1:3936266/150‑1 (MQ=255)
                                                                   tGGCGGCAGCGAACCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCGATTACCACTTCGATTGTACCGGCACGCTTGCCCCAAcc                                                                               >  1:1588951/1‑150 (MQ=255)
                                                                       ggcAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTc                                                                           >  1:1886485/1‑150 (MQ=255)
                                                                          aGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGTGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTg                                                                        >  1:1052036/1‑150 (MQ=255)
                                                                          aGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTg                                                                        >  2:2315305/1‑150 (MQ=255)
                                                                           gCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTgg                                                                       <  2:3279392/150‑1 (MQ=255)
                                                                           gCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTgg                                                                       <  1:5261987/150‑1 (MQ=255)
                                                                           gCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTgg                                                                       <  2:6215266/150‑1 (MQ=255)
                                                                           gCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTgg                                                                       <  2:4601713/150‑1 (MQ=255)
                                                                            cGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATTGCACGTGGCCCAGTGCCGTTGGCGTACATCGGCTCACCAAATACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTggg                                                                      >  1:2877162/1‑150 (MQ=255)
                                                                            cGATCCCATGGCCCGCTGGGTCGCTTCGGTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTggg                                                                      <  2:5954008/150‑1 (MQ=255)
                                                                             gATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGc                                                                     <  2:6569820/150‑1 (MQ=255)
                                                                                cccATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAAATACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAg                                                                  >  1:5899728/1‑150 (MQ=255)
                                                                                  cATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGaa                                                                >  2:10245/1‑150 (MQ=255)
                                                                                  cATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGaa                                                                <  2:1974076/150‑1 (MQ=255)
                                                                                  cATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGaa                                                                <  1:6458469/150‑1 (MQ=255)
                                                                                  cATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGaa                                                                <  1:1403340/150‑1 (MQ=255)
                                                                                  cATGGCCCGCTGGGTCACTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGaa                                                                <  2:1527916/150‑1 (MQ=255)
                                                                                    tGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACATCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACt                                                              >  1:4052718/1‑150 (MQ=255)
                                                                                     ggCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACtt                                                             <  2:5775049/150‑1 (MQ=255)
                                                                                        ccGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGcc                                                          >  1:3247064/1‑150 (MQ=255)
                                                                                         cGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCg                                                         <  1:6176441/150‑1 (MQ=255)
                                                                                             gggTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGt                                                     >  2:3759539/1‑150 (MQ=255)
                                                                                              ggTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGtt                                                    <  1:1937945/150‑1 (MQ=255)
                                                                                              ggTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGAATGTAACGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGtt                                                    >  1:6544931/1‑150 (MQ=255)
                                                                                               gTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTtg                                                   <  2:3716223/150‑1 (MQ=255)
                                                                                                 cGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTtgtg                                                 <  2:4280422/150‑1 (MQ=255)
                                                                                                  gCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTgg                                                <  2:6482051/150‑1 (MQ=255)
                                                                                                  gCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTgg                                                <  1:120249/150‑1 (MQ=255)
                                                                                                    ttCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGc                                              >  1:6839625/1‑150 (MQ=255)
                                                                                                    ttCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGc                                              >  2:5309005/1‑150 (MQ=255)
                                                                                                    ttCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACATGCCGGCGTTGTGGGc                                              >  1:6383360/1‑150 (MQ=255)
                                                                                                        ttCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATc                                          <  1:4000112/150‑1 (MQ=255)
                                                                                                           cAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGgc                                       >  1:6621982/1‑150 (MQ=255)
                                                                                                           cAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGgc                                       >  1:5655864/1‑150 (MQ=255)
                                                                                                           cAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGgc                                       >  2:2652961/1‑150 (MQ=255)
                                                                                                                tGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCAc                                  <  2:1773958/150‑1 (MQ=255)
                                                                                                                tGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCAc                                  <  2:2495087/150‑1 (MQ=255)
                                                                                                                 gCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCAcc                                 <  1:1458745/150‑1 (MQ=255)
                                                                                                                    gcgcGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTAGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGATGTGGGCAATCGGCAGCACCGGc                              >  2:3476732/1‑150 (MQ=255)
                                                                                                                       cgGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTAACGGCACGCTTTCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGc                           >  1:4177176/1‑150 (MQ=255)
                                                                                                                        ggTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGcc                          >  1:4781056/1‑150 (MQ=255)
                                                                                                                          tCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCg                        >  1:675848/1‑150 (MQ=255)
                                                                                                                             ttGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCa                     <  2:1246485/150‑1 (MQ=255)
                                                                                                                             ttGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCa                     >  2:6734404/1‑150 (MQ=255)
                                                                                                                                  ttCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCAc                <  1:4963116/150‑1 (MQ=255)
                                                                                                                                  ttCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCAc                <  1:6666257/150‑1 (MQ=255)
                                                                                                                                  ttCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCAc                >  2:4677233/1‑150 (MQ=255)
                                                                                                                                     ggCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGc             >  1:2255759/1‑150 (MQ=255)
                                                                                                                                       cAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCa           >  1:1150932/1‑150 (MQ=255)
                                                                                                                                       cAACGGCACGGGGCCACGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCa           >  1:2981642/1‑150 (MQ=255)
                                                                                                                                        aaTGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCgg                               <  2:3975136/129‑1 (MQ=255)
                                                                                                                                          tGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAggg        >  2:6818857/1‑150 (MQ=255)
                                                                                                                                           ggCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGc       <  2:2816642/150‑1 (MQ=255)
                                                                                                                                            gCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCGATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCg      <  2:1588951/150‑1 (MQ=255)
                                                                                                                                            gCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCg      <  2:3064754/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ATCCAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCGGTAC  >  NC_002947/68453‑68746

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