Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 68,603 C→A 7.0% G224W (GGG→TGG)  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,6030CA 7.0% 467.9 / 18.6 142G224W (GGG→TGG) PP_00581‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (72/60);  new base A (4/6);  total (76/66)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 146 total aligned reads displayed.

TCCAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCGGTACCGCCGC  >  NC_002947/68454‑68752
                                                                                                                                                     |                                                                                                                                                     
  cAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCg                                                                                                                                                     >  1:3716013/1‑150 (MQ=255)
  cAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCg                                                                                                                                                     >  1:494720/1‑150 (MQ=255)
   aaaaaCTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGt                                                                                                                                                    >  1:4702272/1‑150 (MQ=255)
   aaaaaCTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGt                                                                                                                                                    >  2:440582/1‑150 (MQ=255)
     aaaCTGCACACAAGGTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCCTGGCCCGCTGGGTCGCTTCGGTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGGGGCCACGTGc                                                                                                                                                  <  2:3163893/150‑1 (MQ=255)
        cTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGt                                                                                                                                               <  2:2471574/150‑1 (MQ=255)
         tGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGtt                                                                                                                                              >  1:1977983/1‑150 (MQ=255)
          gCACACAAGTTATCCCCCGACGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCACGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTg                                                                                                                                             <  2:4349279/150‑1 (MQ=255)
          gCACACAAGTTATCCACCGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTg                                                                                                                                             <  2:4349336/150‑1 (MQ=255)
          gCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTg                                                                                                                                             <  2:2006251/150‑1 (MQ=255)
            acacaAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGc                                                                                                                                           <  1:4541257/150‑1 (MQ=255)
             cacaAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCg                                                                                                                                          >  1:401189/1‑150 (MQ=255)
              acaAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGt                                                                                                                                         >  1:848377/1‑150 (MQ=255)
               caAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGCCAACGGCACGGGGCCACGTGCCATTGCCgaa                                                                                                                                        <  2:2150243/150‑3 (MQ=255)
                 aGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACa                                                                                                                                      >  1:2679548/1‑150 (MQ=255)
                     aTCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCCTGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCgg                                                                                                                                  <  1:3982404/150‑1 (MQ=255)
                     aTCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCGGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGCCAACGGCACGGGGCCACGTGCCATTGCCGAActcaac                                                                                                                                  <  1:4707505/150‑6 (MQ=38)
                     aTCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGCCAACGGCACGGGGCCACGTGCCATTGCCGAActcaac                                                                                                                                  <  2:1645372/150‑6 (MQ=38)
                      tCCACAGCAGGTCAAGCCATCTGCTGCTGAGGTTTTTCCGGAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCACGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGc                                                                                                                                 <  2:4152181/150‑1 (MQ=255)
                      tCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGc                                                                                                                                 <  2:1362046/150‑1 (MQ=255)
                         acaGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCa                                                                                                                              >  2:484696/1‑150 (MQ=255)
                         acaGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCa                                                                                                                              >  1:4222340/1‑150 (MQ=255)
                         acaGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCa                                                                                                                              >  2:2896880/1‑150 (MQ=255)
                              aGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAAt                                                                                                                         <  2:3159572/150‑1 (MQ=255)
                               ggTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAAtt                                                                                                                        <  2:203339/150‑1 (MQ=255)
                               ggTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAAtt                                                                                                                        <  1:4510766/150‑1 (MQ=255)
                               ggTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCGCACCAAtt                                                                                                                        >  1:4917295/1‑150 (MQ=255)
                                                cTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTAcc                                                                                                       >  2:2636880/1‑150 (MQ=255)
                                                cTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTAcc                                                                                                       >  2:4645988/1‑150 (MQ=255)
                                                 tGAGGTTTTTCCGCAGCTGGTGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCg                                                                                                      >  1:1586744/1‑150 (MQ=255)
                                                 tGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCg                                                                                                      <  1:98559/150‑1 (MQ=255)
                                                 tGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCg                                                                                                      >  1:652056/1‑150 (MQ=255)
                                                 tGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCg                                                                                                      <  1:549436/150‑1 (MQ=255)
                                                 ngAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCg                                                                                                      >  1:745597/2‑150 (MQ=255)
                                                  gAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCgg                                                                                                     >  2:512185/1‑150 (MQ=255)
                                                   aGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGc                                                                                                    <  1:3329590/150‑1 (MQ=255)
                                                     gTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCAc                                                                                                  <  1:497100/150‑1 (MQ=255)
                                                         ttCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCtt                                                                                              <  1:736116/150‑1 (MQ=255)
                                                          tCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTg                                                                                             <  2:4796154/150‑1 (MQ=255)
                                                           ccGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGc                                                                                            <  1:3951908/150‑1 (MQ=255)
                                                             gCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGccc                                                                                          >  2:3428196/1‑150 (MQ=255)
                                                             gCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGccc                                                                                          >  1:2481288/1‑150 (MQ=255)
                                                             gCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGccc                                                                                          >  1:4637752/1‑150 (MQ=255)
                                                             gCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGccc                                                                                          >  2:1059557/1‑150 (MQ=255)
                                                              cAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGAGTGTACCGGCACGCTTGcccc                                                                                         >  1:4673480/1‑150 (MQ=255)
                                                               aGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCa                                                                                        >  1:3581306/1‑150 (MQ=255)
                                                                 cTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCTTTGGCGTACATCGGCTCACCAATTACCACTTCGTTTGTACCGGCACGCTTGCCCCAAc                                                                                      <  1:690493/150‑1 (MQ=255)
                                                                         aGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTg                                                                              <  1:4463645/150‑1 (MQ=255)
                                                                         aGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTg                                                                              <  2:3275850/150‑1 (MQ=255)
                                                                         aGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTg                                                                              <  2:2359514/150‑1 (MQ=255)
                                                                         aGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTg                                                                              >  1:4973837/1‑150 (MQ=255)
                                                                         aGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTAACGGCACGCTTGCCCCAACCGGTCTTg                                                                              >  2:1286107/1‑150 (MQ=255)
                                                                         aGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATAGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTg                                                                              >  1:1417075/1‑150 (MQ=255)
                                                                            gATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACATCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGc                                                                           >  1:867205/1‑150 (MQ=255)
                                                                             aTCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGcc                                                                          >  2:3308218/1‑150 (MQ=255)
                                                                              tCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCa                                                                         >  2:4474929/1‑150 (MQ=255)
                                                                                ccATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGa                                                                       >  2:3051072/1‑150 (MQ=255)
                                                                                 cATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGaa                                                                      >  1:3836747/1‑150 (MQ=255)
                                                                                 cATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGaa                                                                      >  2:34540/1‑150 (MQ=255)
                                                                                  aTGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAAATACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCGTGGGCCAGAAc                                                                     >  1:3073590/1‑150 (MQ=255)
                                                                                    ggCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACtt                                                                   <  2:2674375/150‑1 (MQ=255)
                                                                                      cccGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGc                                                                 <  2:4658564/150‑1 (MQ=255)
                                                                                        cGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCg                                                               <  1:4822949/150‑1 (MQ=255)
                                                                                            gggTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGt                                                           <  2:3718894/150‑1 (MQ=255)
                                                                                              gTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTtg                                                         <  1:2154109/150‑1 (MQ=255)
                                                                                               tCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTtgt                                                        <  1:1266196/150‑1 (MQ=255)
                                                                                                 gCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTgg                                                      <  1:1783769/150‑1 (MQ=255)
                                                                                                 gCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTgg                                                      >  1:4349805/1‑150 (MQ=255)
                                                                                                 gCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTgg                                                      <  2:2679548/150‑1 (MQ=255)
                                                                                                   ttCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGc                                                    >  1:1215500/1‑150 (MQ=255)
                                                                                                   ttCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGc                                                    >  2:4547277/1‑150 (MQ=255)
                                                                                                       ttCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATc                                                <  1:3872245/150‑1 (MQ=255)
                                                                                                         ccAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCgg                                              >  1:3905198/1‑150 (MQ=255)
                                                                                                         ccAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCgg                                              >  2:1428851/1‑150 (MQ=255)
                                                                                                           aaGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGgca                                            <  2:2696471/150‑1 (MQ=255)
                                                                                                           aaGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGgca                                            >  1:1850190/1‑150 (MQ=255)
                                                                                                                         tCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCg                              <  2:133076/150‑1 (MQ=255)
                                                                                                                           gTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCgg                                                                                                     >  2:4919121/1‑77 (MQ=255)
                                                                                                                           gTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCgg                                                                                                     <  1:4919121/77‑1 (MQ=255)
                                                                                                                           gTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGc                            >  2:4033156/1‑150 (MQ=255)
                                                                                                                            ttGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCa                           >  2:4236445/1‑150 (MQ=255)
                                                                                                                            ttGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCa                           <  1:2681257/150‑1 (MQ=255)
                                                                                                                                 ttCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCAc                      >  1:3442973/1‑150 (MQ=255)
                                                                                                                                   cGGCAATGGCACGTTGCCCCGTTCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACgg                    <  2:2772660/150‑1 (MQ=255)
                                                                                                                                    ggCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGc                   <  2:2627343/150‑1 (MQ=255)
                                                                                                                                     ccAACGGCACGGGGCCACGTGCCATTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGcc                  >  1:4595761/2‑150 (MQ=255)
                                                                                                                                      cAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCa                 >  1:3590297/1‑150 (MQ=255)
                                                                                                                                      cAACGGCACGGGGCCACGTGCCATTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCa                 >  1:3286223/1‑150 (MQ=255)
                                                                                                                                       aacggCACGGGGCCACGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAg                >  1:1716646/4‑150 (MQ=255)
                                                                                                                                         tGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAggg              <  1:1167002/150‑1 (MQ=255)
                                                                                                                                          ggCACGGGGCCACGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGc             >  1:3143837/1‑150 (MQ=255)
                                                                                                                                            cACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCgg           >  1:249822/1‑150 (MQ=255)
                                                                                                                                             aCGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCGGt          >  1:4624858/1‑150 (MQ=255)
                                                                                                                                             aCGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCGGt          >  1:4624903/1‑150 (MQ=255)
                                                                                                                                             aCGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCGGt          >  2:5002142/1‑150 (MQ=255)
                                                                                                                                               gTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCGGTAc        <  1:4097250/150‑1 (MQ=255)
                                                                                                                                                 ggCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCGGTAccg      <  1:1795603/150‑1 (MQ=255)
                                                                                                                                                  gCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCGGTAccgc     >  1:4493850/1‑150 (MQ=255)
                                                                                                                                                     ccGTGCCGTTGGCGTACATCGGCTCCCCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCGGTACcgccgc  <  2:5096935/150‑1 (MQ=255)
                                                                                                                                                     ccGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCGGTACcgccgc  >  1:2560144/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCAAAAACTGCACACAAGTTATCCACAGAAGGTCAAGCCATCTGCTCCTGAGGTTTTTCCGCAGCTGGCGGCAGCGATCCCATGGCCCGCTGGGTCGCTTCGTTCCAAGCTGCAGCGCGGTCGTTGAGTTCGGCAATGGCACGTGGCCCCGTGCCGTTGGCGTACATCGGCTCACCAATTACCACTTCGATTGTACCGGCACGCTTGCCCCAACCGGTCTTGGGCCAGAACTTGCCGGCGTTGTGGGCAATCGGCAGCACCGGCAGCCCGGCATTCACGGCCAGGGCGGTACCGCCGC  >  NC_002947/68454‑68752

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