Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,685,384 G→T 7.8% intergenic (+23/‑36) PP_5593 → / → PP_5594 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,685,3840GT 7.8% 195.5 / 5.8 64intergenic (+23/‑36)PP_5593/PP_5594hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/34);  new base T (2/3);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.24e-01

AAGAACCCGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTAAAGATTCCGACCC  >  NC_002947/3685236‑3685517
                                                                                                                                                    |                                                                                                                                     
aaGAACCCGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGc                                                                                                                                      >  2:1514491/1‑150 (MQ=255)
aaGAACCCGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGc                                                                                                                                      <  1:470096/150‑1 (MQ=255)
aaGAACCCGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCAAGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGc                                                                                                                                      <  2:1532498/150‑1 (MQ=255)
    aCCCGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGAcc                                                                                                                                  <  1:835517/150‑1 (MQ=255)
     cccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGc                                                                                                                                      <  2:985200/145‑1 (MQ=255)
     cccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGc                                                                                                                                      >  1:985200/1‑145 (MQ=255)
     cccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCt                                                                                                                                 >  2:1794331/1‑150 (MQ=255)
     cccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCt                                                                                                                                 <  1:1836110/150‑1 (MQ=255)
      ccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTg                                                                                                                                >  1:1457553/1‑150 (MQ=255)
      ccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTg                                                                                                                                <  1:659947/150‑1 (MQ=255)
      ccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTg                                                                                                                                <  2:1688412/150‑1 (MQ=255)
      ccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGAAGTTACGGAGCGAccat                                                                                                                                >  2:1564557/1‑148 (MQ=255)
       cGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTgg                                                                                                                               <  1:1737561/150‑1 (MQ=255)
            aaTTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCAGGTCGCCCCGGCACTg                                                                                                                          <  2:44124/150‑1 (MQ=255)
                    cAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATc                                                                                                                  >  2:725702/1‑150 (MQ=255)
                    cAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCAGGTCGCCCCGGCACTGcgc                                                                                                                       >  1:344882/1‑143 (MQ=255)
                    cAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCAGGTCGCCCCGGCACTGcgc                                                                                                                       <  2:344882/145‑3 (MQ=255)
                        gACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATcc                                                                                                                 <  2:444025/147‑1 (MQ=255)
                        gACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATcc                                                                                                                 >  1:444025/1‑147 (MQ=255)
                         aCCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACTGAGGCGCAGTGCCGGGGCGACATGGCACTGCAATTATCCAAgg                                                                                                             <  2:275201/150‑1 (MQ=255)
                           cGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAgggg                                                                                                           <  2:1925457/150‑1 (MQ=255)
                           cGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAgggg                                                                                                           >  1:1520962/1‑150 (MQ=255)
                             tGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTgg                                                                                                                               >  2:1633268/1‑128 (MQ=255)
                             tGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTgg                                                                                                                               <  1:1633268/128‑1 (MQ=255)
                                        ccGACGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCAACAAGGACGGCCGTTCGTAACTGACCGAGTCTTAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATgcgc                                                                                              <  1:725702/150‑1 (MQ=255)
                                                 gTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCaa                                                                                                                       >  2:886213/1‑116 (MQ=255)
                                                 gTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCaa                                                                                                                       <  1:886213/116‑1 (MQ=255)
                                                  tCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCAc                                                                                    <  1:1564557/150‑1 (MQ=255)
                                                   cGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACt                                                                                   <  2:104050/150‑1 (MQ=255)
                                                        aGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGa                                                                                 <  2:247053/147‑1 (MQ=255)
                                                        aGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGa                                                                                 >  1:247053/1‑147 (MQ=255)
                                                         ggcggcCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcat                                                                             <  1:431262/150‑1 (MQ=255)
                                                          gcggcCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            <  1:357792/150‑1 (MQ=255)
                                                          gcggcCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            <  2:1457553/150‑1 (MQ=255)
                                                              ccAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            >  2:254098/1‑146 (MQ=255)
                                                              ccAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            <  1:254098/146‑1 (MQ=255)
                                                              ccAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCAGGTCGCCCCGGCACTGcg                                                                                                                        >  1:709126/1‑101 (MQ=38)
                                                              ccAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCAGGTCGCCCCGGCACTGcg                                                                                                                        <  2:709126/102‑2 (MQ=38)
                                                                           ggcAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTgcg                                                           >  2:1698902/1‑148 (MQ=255)
                                                                               aCATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTg                                                       <  1:1794010/150‑1 (MQ=255)
                                                                                  tGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCt                                                                                            <  2:1095197/110‑1 (MQ=255)
                                                                                  tGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCt                                                                                            >  1:1095197/1‑110 (MQ=255)
                                                                                       cACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTcggcgg                                               >  2:1584896/1‑150 (MQ=255)
                                                                                             cgcgAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATGATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGttt                                         >  1:1105209/1‑150 (MQ=255)
                                                                                               cgAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTcggcg                                                <  2:1465636/141‑1 (MQ=255)
                                                                                               cgAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTcggcg                                                >  1:1465636/1‑141 (MQ=255)
                                                                                                   tCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACcg                                   >  2:862278/1‑150 (MQ=255)
                                                                                                     gggCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGa                                                            >  2:224630/1‑123 (MQ=255)
                                                                                                     gggCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGa                                                            <  1:224630/123‑1 (MQ=255)
                                                                                                      ggCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGa                                >  2:404611/1‑150 (MQ=255)
                                                                                                         cacaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            <  2:1345835/103‑1 (MQ=255)
                                                                                                         cacaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            >  1:1345835/1‑103 (MQ=255)
                                                                                                          acaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTgg                                                                                                                               <  2:531596/51‑1 (MQ=255)
                                                                                                          acaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTgg                                                                                                                               >  1:531596/1‑51 (MQ=255)
                                                                                                          acaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCCTCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAAt                            <  2:1520962/150‑1 (MQ=255)
                                                                                                          acaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAAt                            <  2:1426930/150‑1 (MQ=255)
                                                                                                            aaGGGCGGCCGTTCGTAACTGACCGAGGAGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGa                          >  2:1677846/1‑150 (MQ=255)
                                                                                                             aGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGac                         <  1:1794331/150‑1 (MQ=255)
                                                                                                                       ttCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCg                                       >  2:1273805/1‑126 (MQ=255)
                                                                                                                       ttCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCg                                       <  1:1273805/126‑1 (MQ=255)
                                                                                                                        tCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTaa              >  1:1931327/1‑150 (MQ=255)
                                                                                                                              cTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTAAAGATTc        <  1:712385/150‑1 (MQ=255)
                                                                                                                                  ccGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTAAAGATTCCGAc    <  1:862278/150‑1 (MQ=255)
                                                                                                                                    gAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTAAAGATTCCGAccc  <  2:1931327/150‑1 (MQ=255)
                                                                                                                                                  gggCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAg                                                           <  1:892836/79‑1 (MQ=255)
                                                                                                                                                  gggCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAg                                                           <  1:893037/79‑1 (MQ=255)
                                                                                                                                                  gggCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAg                                                           >  2:892836/1‑79 (MQ=255)
                                                                                                                                                  gggCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAg                                                           >  2:893037/1‑79 (MQ=255)
                                                                                                                                                    |                                                                                                                                     
AAGAACCCGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTAAAGATTCCGACCC  >  NC_002947/3685236‑3685517

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