Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,021,409 C→G 5.8% *314S (TGA→TCA)  PP_3547 ← NAD(P)‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,021,4090CG 5.8% 218.6 / 4.5 69*314S (TGA→TCA) PP_3547NAD(P)‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (34/31);  new base G (2/2);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

CTCAAGGTGCTGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTCCGGGTCA  >  NC_002947/4021260‑4021558
                                                                                                                                                     |                                                                                                                                                     
cTCAAGGTGCTGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTc                                                                                                                                                       >  1:60348/1‑150 (MQ=255)
cTCAAGGTGCTGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTc                                                                                                                                                       >  2:1254041/1‑150 (MQ=255)
cTCAAGGTGCTGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTc                                                                                                                                                       >  1:345355/1‑150 (MQ=255)
 tCAAGGTGCTGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCa                                                                                                                                                      <  1:1199032/150‑1 (MQ=255)
   aaGGTGCTGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTGAgc                                                                                                                                                    >  2:128791/1‑150 (MQ=255)
    aGGTGCTGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCATCGACGGCCAACCGCTGGAAGGTACGTTTGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAgcg                                                                                                                                                   <  2:649481/150‑1 (MQ=255)
        gCTGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTGAGcgcccg                                                                                                                                               <  1:128791/150‑1 (MQ=255)
           gTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTGAGcgcccgtca                                                                                                                                            >  1:562645/1‑147 (MQ=255)
            ttGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTGAGcgcccgtcaa                                                                                                                                           <  2:562645/150‑5 (MQ=255)
               aTCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCa                                                                                                                                        >  1:56826/1‑150 (MQ=255)
                 cACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAgc                                                                                                                                                    >  1:625523/1‑136 (MQ=255)
                 cACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAgc                                                                                                                                                    <  2:625523/136‑1 (MQ=255)
                  aCCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGACAGAGCGGCTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCCGCGCAGGCTGCCCatgc                                                                                                                                     >  2:659859/1‑147 (MQ=255)
                   ccGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTc                                                                                                                                    >  2:1262982/1‑150 (MQ=255)
                        tACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGa                                                                                                                               <  1:863204/150‑1 (MQ=255)
                                    gcggcgTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTgg                                                                                                                   <  2:1227911/150‑1 (MQ=255)
                                       gcgTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACTGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCAt                                                                                                                >  1:414294/1‑150 (MQ=255)
                                       gcgTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCAt                                                                                                                >  2:1190556/1‑150 (MQ=255)
                                       gcgTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCAt                                                                                                                <  2:105223/150‑1 (MQ=255)
                                              gCGACGGCCAACCGTTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGc                                                                                                         <  1:659859/150‑1 (MQ=255)
                                                     ccAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTcc                                                                                                  <  2:345355/150‑1 (MQ=255)
                                                     ccAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTcc                                                                                                  >  1:104726/1‑150 (MQ=255)
                                                      cAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCg                                                                                                 >  2:845404/1‑150 (MQ=255)
                                                         ccGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGaaa                                                                                              <  2:56826/150‑1 (MQ=255)
                                                         ccGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGaaa                                                                                              <  2:60348/150‑1 (MQ=255)
                                                         ccGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGaaa                                                                                              >  2:1052711/1‑150 (MQ=255)
                                                         ccGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGaaa                                                                                              <  2:970575/150‑1 (MQ=255)
                                                              gCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAAc                                                                                         <  1:892139/150‑1 (MQ=255)
                                                                     aCGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGACCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGa                                                                                  <  2:1027336/150‑1 (MQ=255)
                                                                       gTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATg                                                                                >  1:1042740/1‑150 (MQ=255)
                                                                       gTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATg                                                                                >  1:1190864/1‑150 (MQ=255)
                                                                        ttGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCttt                                                                                                      <  2:48802/127‑1 (MQ=255)
                                                                        ttGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCttt                                                                                                      >  1:48802/1‑127 (MQ=255)
                                                                        ttGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATgg                                                                               <  1:423460/150‑1 (MQ=255)
                                                                           gtgtTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCt                                                                            <  2:657920/150‑1 (MQ=255)
                                                                               tGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCTTTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGc                                                                        <  1:747833/150‑1 (MQ=255)
                                                                                 aCCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCa                                                                      >  2:1165044/1‑150 (MQ=255)
                                                                                          ttCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCgg                                                             >  2:781420/1‑150 (MQ=255)
                                                                                             gAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGt                                                                                                                     <  1:765908/91‑1 (MQ=255)
                                                                                             gAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGt                                                                                                                     >  2:765908/1‑91 (MQ=255)
                                                                                             gAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGAc                                                          >  2:251233/1‑150 (MQ=255)
                                                                                               aCTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGc                                                        <  1:978003/150‑1 (MQ=255)
                                                                                                cTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGcc                                                       >  2:1381946/1‑150 (MQ=255)
                                                                                                            gCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCg                                                                                                                       <  2:1121041/74‑1 (MQ=255)
                                                                                                            gCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCg                                                                                                                       >  1:1121041/1‑74 (MQ=255)
                                                                                                            gCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATa                                                                                                               >  1:225591/1‑82 (MQ=255)
                                                                                                            gCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATa                                                                                                               <  2:225591/82‑1 (MQ=255)
                                                                                                            gCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCa                                           >  2:551844/1‑150 (MQ=255)
                                                                                                            gCAGAGCGGGTTCGCAAGCTGCTATAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCa                                           >  2:837282/1‑150 (MQ=255)
                                                                                                              agagCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATg                                         >  2:607622/1‑150 (MQ=255)
                                                                                                               gagCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAgg                                                                                                                             >  2:938532/1‑65 (MQ=255)
                                                                                                               gagCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAgg                                                                                                                             <  1:938532/65‑1 (MQ=255)
                                                                                                                   gggTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCaa                                    <  1:1190556/150‑1 (MQ=255)
                                                                                                                     gTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGc                                  >  1:568403/1‑150 (MQ=255)
                                                                                                                      ttCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGcc                                 <  1:781420/150‑1 (MQ=255)
                                                                                                                               tgctgGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAg                        <  1:1052711/150‑1 (MQ=255)
                                                                                                                                gctgGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGc                       >  1:353872/1‑150 (MQ=255)
                                                                                                                                      aCCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCTGGTGGCATAGGAGGCtt                                                                                                       <  2:727479/64‑1 (MQ=255)
                                                                                                                                      aCCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCtt                                                                                                       >  1:727479/1‑64 (MQ=255)
                                                                                                                                       ccTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCCCTTCTTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGa                >  1:1150731/1‑150 (MQ=255)
                                                                                                                                        cTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGAc               <  1:1381946/150‑1 (MQ=255)
                                                                                                                                        cTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGAc               >  1:826351/1‑150 (MQ=255)
                                                                                                                                        cTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGAc               >  1:637140/1‑150 (MQ=255)
                                                                                                                                        cTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCAGCCAATGGGCAAGCCAGGTTCGAGCTTGCTGAc               >  2:227094/1‑150 (MQ=255)
                                                                                                                                          tGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGAcgc             <  1:551844/150‑1 (MQ=255)
                                                                                                                                               ggCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTcc        <  1:1262982/150‑1 (MQ=255)
                                                                                                                                                  gcTCAGCGCCGGCCGCAGCGGCCGATGAGGCGGCCGGGGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTCCggg     <  2:1150731/150‑1 (MQ=255)
                                                                                                                                                    tCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATg                                                                                >  2:1047811/1‑73 (MQ=255)
                                                                                                                                                    tCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATg                                                                                <  1:1047811/73‑1 (MQ=255)
                                                                                                                                                    tCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTCCGGGTc   >  2:861915/1‑150 (MQ=255)
                                                                                                                                                    tCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTCCGGGTc   >  1:78008/1‑150 (MQ=255)
                                                                                                                                                    tCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGACAAGCCAGGTTCGAGCTTGCTGACGCGGTCCGGGTc   >  2:583332/1‑150 (MQ=255)
                                                                                                                                                     cAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACCCGGACCGCGTCa  <  1:837282/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCAAGGTGCTGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTCCGGGTCA  >  NC_002947/4021260‑4021558

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