Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,977,148 C→A 5.3% A271E (GCG→GAG)  PP_1770 → 3‑phosphoshikimate 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,977,1480CA 5.3% 279.9 / 3.8 85A271E (GCG→GAG) PP_17703‑phosphoshikimate 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base C (42/38);  new base A (2/3);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

TTCCGGTACGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCAACGATCTGATTTTCCAGGCCCAAC  >  NC_002947/1977000‑1977288
                                                                                                                                                        |                                                                                                                                             
ttCCGGTACGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGcgcgacgc                                                                                                                                               >  3:228267/1‑149 (MQ=255)
     gTACGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCg                                                                                                                                          <  7:122918/149‑1 (MQ=255)
      tACGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGc                                                                                                                                         <  4:262699/149‑1 (MQ=255)
        cGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCt                                                                                                                                        <  2:213151/148‑1 (MQ=255)
        cGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCGCGTCGCGCAGGGCGTCGAGATCACt                                                                                                                                        <  4:290456/148‑1 (MQ=255)
         gCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGCGATCGCGTCGCGCAGGGCGTCGAGATCACTgc                                                                                                                                      <  3:174599/149‑2 (MQ=38)
          cTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTg                                                                                                                                       >  6:346594/1‑147 (MQ=255)
             cGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGTGCGACGCGATCGCTGAAg                                                                                                                                    <  2:75311/147‑1 (MQ=255)
             cGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAggg                                                                                                                                  <  5:346594/149‑1 (MQ=255)
             cGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGCGATCGCGTCGCGCAGGGCGTCGAGATCACTgcgaaa                                                                                                                                  <  2:388730/149‑6 (MQ=37)
               ggtgggtTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCCGTCCACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAggg                                                                                                                                  <  1:330684/149‑1 (MQ=255)
                        cgcgATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGAt                                                                                                                                            <  3:187872/128‑1 (MQ=255)
                        cgcgATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGAt                                                                                                                                            >  4:187872/1‑128 (MQ=255)
                        cgcgATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCa                                                                                                                       <  5:411648/149‑1 (MQ=255)
                         gcgATTTCACAAGAATCGCCGGCAGCGATCCTACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCg                                                                                                                                          <  8:67124/129‑1 (MQ=255)
                         gcgATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCg                                                                                                                                          >  7:67124/1‑129 (MQ=255)
                         gcgATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGAACCAg                                                                                                                     <  8:165414/149‑1 (MQ=255)
                           gATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAgct                                                                                                                   >  2:37879/1‑149 (MQ=255)
                                 acaAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGGGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTgg                                                                                                              <  7:50336/148‑1 (MQ=255)
                                 acaAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTgg                                                                                                              <  2:401949/148‑1 (MQ=255)
                                  caAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGgt                                                                                                            <  2:183615/149‑1 (MQ=255)
                                   aaGAATCGCCGGCAGCTATCCGACCATGTGGCATGACATCTTCGTCGCCTACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTTATCTCGACGCCCTGCGCGACGCTATCGCTGAAGGGGATGGGCACCAGCTGCtgtgtg                                                                                                           <  4:228267/149‑1 (MQ=255)
                                   aaGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGgtg                                                                                                           >  8:259627/1‑149 (MQ=255)
                                        tCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCTCTGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTAt                                                                                                        <  1:325221/149‑1 (MQ=255)
                                        tCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTc                                                                                                      <  8:112014/149‑1 (MQ=255)
                                           ccGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCAcc                                                                                                   >  7:132881/1‑149 (MQ=255)
                                           ccGGCAGCGATCCGACAATGTGGCATGACATCTTCCTAGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATATCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACAAGCTGCTGGGTGTATTCAcc                                                                                                   >  2:175206/1‑149 (MQ=255)
                                              gCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCcgc                                                                                                >  5:318061/1‑149 (MQ=255)
                                               cAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGAATACAcccagc                                                                                               <  7:105814/149‑4 (MQ=255)
                                                 gCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCa                                                                                             >  7:229584/1‑149 (MQ=255)
                                                      ccGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGt                                                                                        >  3:55864/1‑149 (MQ=255)
                                                        gACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTg                                                                                      >  2:314335/1‑149 (MQ=255)
                                                         aCCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTTTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGc                                                                                     <  6:48155/149‑1 (MQ=255)
                                                         aCCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGc                                                                                     <  1:37879/149‑1 (MQ=255)
                                                           cATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCAACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGccc                                                                                   >  5:46695/1‑149 (MQ=255)
                                                             tgtgGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTAt                                                                                                        >  8:1176/1‑126 (MQ=255)
                                                             tgtgGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTAt                                                                                                        <  7:1176/126‑1 (MQ=255)
                                                             tgtgGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCcgc                                                                                 <  4:55864/149‑1 (MQ=255)
                                                                  cATGACAACTTCCTCTACAACCGCGACGCTGTCCTGCGAACGCATGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCTACGCGATAGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGACCGCGAGCa                                                                            >  6:284917/1‑149 (MQ=255)
                                                                      aCATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCAttt                                                                         <  8:229584/148‑1 (MQ=255)
                                                                      aCATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCAttt                                                                         <  8:132881/148‑1 (MQ=255)
                                                                      aCATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCAttt                                                                         <  6:90615/148‑1 (MQ=255)
                                                                       cATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCa                                                                       <  6:205821/149‑1 (MQ=255)
                                                                          cTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGt                                                                     >  8:302894/1‑148 (MQ=255)
                                                                           ttCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTa                                                                    >  2:377023/1‑148 (MQ=255)
                                                                           ttCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTa                                                                    >  6:61136/1‑148 (MQ=255)
                                                                             ccTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTa                                                                                                         >  2:400580/1‑109 (MQ=255)
                                                                             ccTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTa                                                                                                         <  1:400580/109‑1 (MQ=255)
                                                                             ccTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTaaa                                                                  <  6:291095/148‑1 (MQ=255)
                                                                              cTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTaaaa                                                                 <  3:112247/148‑1 (MQ=255)
                                                                               tCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATc                                                               >  1:23911/1‑149 (MQ=255)
                                                                                 gCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCt                                                             >  1:291912/1‑149 (MQ=255)
                                                                                  ccAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCACTCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCt                                                             <  1:314335/149‑1 (MQ=255)
                                                                                  ccAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTg                                                            <  4:363720/149‑1 (MQ=255)
                                                                                     aCCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCGGGTGAATACAcccag                                                                                                <  8:276735/109‑3 (MQ=255)
                                                                                     aCCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCGGGTGAATACAcccag                                                                                                >  7:276735/1‑107 (MQ=255)
                                                                                     aCCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGAGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGcc                                                         >  1:172829/1‑149 (MQ=255)
                                                                                               gCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGc                                                 <  6:318061/149‑1 (MQ=255)
                                                                                                   tCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTAt                                             <  2:23911/149‑1 (MQ=255)
                                                                                                     cTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTAtgt                                           >  3:95193/1‑149 (MQ=255)
                                                                                                     cTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTAtgt                                           >  3:318269/1‑149 (MQ=255)
                                                                                                           aCGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGc                                     >  3:185614/1‑149 (MQ=255)
                                                                                                              cTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTAt                                  <  2:404139/149‑1 (MQ=255)
                                                                                                               ttGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTAt                                  >  1:252541/1‑148 (MQ=255)
                                                                                                                 gATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATg                                 <  7:259627/147‑1 (MQ=255)
                                                                                                                 gATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGaa                               >  5:302676/1‑149 (MQ=255)
                                                                                                                  aTACATTTCG‑‑CAGTGATCTCGACGCCCTTCGCGACGCGATCGCTGAAGGGGATGTGCACCAGCTGCTGTGTGTATTCACCCGAGCACGGGTTGCCAGAGAGCATTTCAGTAAAATCCTGGCACGCCGTGACTATGTGGACGCTATGAAc                              >  4:247239/1‑149 (MQ=255)
                                                                                                                  aTACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAAc                              >  1:163170/1‑149 (MQ=255)
                                                                                                                  aTACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAAc                              <  6:302676/149‑1 (MQ=255)
                                                                                                                   tACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACg                             >  7:267616/1‑149 (MQ=255)
                                                                                                                      aTTTCG‑‑CAGTGATCTCGCCGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGcc                           >  6:57891/1‑148 (MQ=255)
                                                                                                                      aTTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGcc                           >  1:257371/1‑148 (MQ=255)
                                                                                                                       tttCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGACCGCCaa                         >  6:201779/1‑149 (MQ=255)
                                                                                                                       tccccttcAGCGATCGCGTCGCGCAGGGCGTCGAGATCACTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGc                            >  2:60052/8‑148 (MQ=37)
                                                                                                                         tCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCAAc                        <  8:68951/148‑1 (MQ=255)
                                                                                                                         tCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCAAc                        >  7:66664/1‑148 (MQ=255)
                                                                                                                           g‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCAACGAt                     >  7:344068/1‑149 (MQ=255)
                                                                                                                                 tGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCAACGATCTGa                 >  4:139762/1‑149 (MQ=255)
                                                                                                                                      tcGACGCCCTGCGGGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCAACGATCTGAGTTTc            >  4:58561/1‑149 (MQ=255)
                                                                                                                                      tcGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCAACGATCTGATTTTc            >  6:78880/1‑149 (MQ=255)
                                                                                                                                      tcGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCAACGATCTGATTTTc            >  7:268928/1‑149 (MQ=255)
                                                                                                                                      tcGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCAACGATCTGATTTTc            >  7:330997/1‑149 (MQ=255)
                                                                                                                                       cGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCAACGATCTGATTTTc            >  4:341341/1‑148 (MQ=255)
                                                                                                                                        gACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCAACGATCTGATTTTc            >  8:26282/1‑147 (MQ=255)
                                                                                                                                        gACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCAACGATCTGATTTTc            <  8:267616/147‑1 (MQ=255)
                                                                                                                                              cTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCAACGATCTGATTTTCCAGGccc     >  6:342430/1‑148 (MQ=255)
                                                                                                                                                 cgcgacgcgaTCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCAACGATCTGATTTTCCAGGCCCAAc  <  1:377023/148‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                             
TTCCGGTACGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCG‑‑CAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCAACGATCTGATTTTCCAGGCCCAAC  >  NC_002947/1977000‑1977288

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