Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,296,544 T→C 7.4% Q125R (CAG→CGG) ‡ PP_2898 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,296,5440TC 7.4% 272.4 / 8.2 95Q125R (CAG→CGG) ‡PP_2898hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (46/41);  new base C (2/5);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

TCGCGATGATGTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAATGTAC  >  NC_002947/3296395‑3296689
                                                                                                                                                     |                                                                                                                                                 
tCGCGATGATGTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCt                                                                                                                                                   >  2:1268495/1‑150 (MQ=255)
    gatgatGTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGcc                                                                                                                                               >  1:182410/1‑150 (MQ=255)
       gatgTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGtt                                                                                                                                            <  2:1964199/150‑1 (MQ=255)
          gtgCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTgcg                                                                                                                                         >  2:944781/1‑150 (MQ=255)
           tgCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTgcgc                                                                                                                                        >  2:1791976/1‑150 (MQ=255)
           tgCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCCAGCCGTTgcac                                                                                                                                        <  1:562643/150‑3 (MQ=255)
           tgCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGCGCAACGGCCAGCCGTTgcac                                                                                                                                        <  1:1080519/150‑3 (MQ=255)
           tgCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGCGCAACGGCCAGCCGTTgcac                                                                                                                                        <  2:1897493/150‑3 (MQ=255)
             cAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGCGGGGTGGCGCGCATGCGCGGCGGCAACGTGCGCCAGTCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTg                                                                                                                                      >  2:266799/1‑150 (MQ=255)
              aGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTgg                                                                                                                                     >  1:2183948/1‑150 (MQ=255)
              aGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTgg                                                                                                                                     >  1:2200612/1‑150 (MQ=255)
              aGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTgg                                                                                                                                     >  1:1008309/1‑150 (MQ=255)
               gTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCGAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGc                                                                                                                                    >  2:199278/1‑150 (MQ=255)
                   gTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGAt                                                                                                                                >  1:1424976/1‑150 (MQ=255)
                    tACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATg                                                                                                                               >  1:1064349/1‑150 (MQ=255)
                    tACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATg                                                                                                                               <  1:2342138/150‑1 (MQ=255)
                          ACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTgg                                                                                                                         <  2:2057451‑M1/150‑1 (MQ=255)
                            ccctcaagcctatctgggtagcgccagtttgctgagcagctacccgtgcgattgaggggtgGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTc                                                                                                                       >  2:971365‑M2/62‑150 (MQ=255)
                               TCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGc                                                                                                                    <  1:1053972‑M1/150‑1 (MQ=255)
                                      cACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCCAGCCGTTgc                                                                                                                                          <  2:177928/121‑1 (MQ=255)
                                      cACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCCAGCCGTTgc                                                                                                                                          >  1:177928/1‑121 (MQ=255)
                                         gggTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTggcgg                                                                                                          <  2:1156203/150‑1 (MQ=255)
                                              gctgctCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGAGGATGCACGGCTGCAACGCGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTa                                                                                                     <  1:2379829/150‑1 (MQ=255)
                                               ctgctCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTAGGGCTGGCGGATTGCCTGGGCTGCAAAGTGGCGGCGGTAc                                                                                                    >  2:2092452/1‑150 (MQ=255)
                                                tgctCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGTCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACa                                                                                                   >  1:1135379/1‑150 (MQ=255)
                                                tgctCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACa                                                                                                   >  2:835911/1‑150 (MQ=255)
                                                 gctCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAg                                                                                                  >  1:86901/1‑150 (MQ=255)
                                                     aGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGTCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTAcc                                                                                              <  2:1746774/150‑1 (MQ=255)
                                                     aGCAAACTGGCGCTACCCAGATAGGCATGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGGGGTGCAGTAcc                                                                                              >  2:1126766/1‑150 (MQ=255)
                                                       cAAACTGGCGCTACCCAGATAGGCGTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGGTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTAcccg                                                                                            >  2:434780/1‑148 (MQ=255)
                                                        aaaCTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTgg                                                                                           <  2:1617606/150‑1 (MQ=255)
                                                         aaCTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGtt                                                                                                                                            >  2:1986806/1‑100 (MQ=255)
                                                         aaCTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGtt                                                                                                                                            <  1:1986806/100‑1 (MQ=255)
                                                           cTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTc                                                                                        <  1:1268495/150‑1 (MQ=255)
                                                              gcgcTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTgg                                                                                           >  1:1571653/1‑144 (MQ=255)
                                                              gcgcTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTgg                                                                                           <  2:1571653/144‑1 (MQ=255)
                                                              gcgcTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGa                                                                                     <  2:2200612/150‑1 (MQ=255)
                                                              gcgcTACACAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGa                                                                                     <  2:669229/150‑1 (MQ=255)
                                                                   aCCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGcc                                                                                                                                               <  1:2133996/87‑1 (MQ=255)
                                                                   aCCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGcc                                                                                                                                               >  2:2133996/1‑87 (MQ=255)
                                                                   aCCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAg                                                                                >  1:1551320/1‑150 (MQ=255)
                                                                    tacccgtgcgattgaggggtgGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGt                                                                               <  1:971365‑M2/129‑1 (MQ=255)
                                                                    cccAGATAGTCTTGAGGGGTGGCGCTCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGt                                                                               <  1:2333483/150‑1 (MQ=255)
                                                                    cccAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGt                                                                               <  2:1817856/150‑1 (MQ=255)
                                                                     ccAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGt                                                                                                               <  1:270182/117‑1 (MQ=255)
                                                                     ccAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACAAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTTCGCTGGCCGATGACCTGGTCGGCAAAGt                                                                                                               >  2:270182/1‑117 (MQ=255)
                                                                      cAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTgg                                                                             >  1:1030602/1‑150 (MQ=255)
                                                                          tAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCtgt                                                                         <  1:835911/150‑1 (MQ=255)
                                                                             gCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGcc                                                                      >  2:351854/1‑150 (MQ=255)
                                                                             gCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGcc                                                                      <  2:1008309/150‑1 (MQ=255)
                                                                             gCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGcc                                                                      >  1:2312339/1‑150 (MQ=255)
                                                                               ttgaggggtgGCGCGCATGCTCGGCGGCAACGTGCGCCAAGCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCGGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTATGATTCAGTGGCTGTGCCTg                                                                    >  2:948362‑M2/11‑150 (MQ=255)
                                                                                tgaggggtgGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTgg                                                                   <  1:606037‑M2/141‑1 (MQ=255)
                                                                                  aggggtgGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTc                                                                 <  2:2312339‑M2/143‑1 (MQ=255)
                                                                                   ggggtgGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCg                                                                <  1:271290‑M2/144‑1 (MQ=255)
                                                                                      gtgGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGcc                                                             <  1:351854‑M2/147‑1 (MQ=255)
                                                                                           gcgcATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCCAGCCGTTgca                                                                                                                                         >  2:581200/1‑68 (MQ=39)
                                                                                           gcgcATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCCAGCCGTTgca                                                                                                                                         <  1:581200/69‑2 (MQ=39)
                                                                                             gcATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGt                                                                               >  1:1586616/1‑125 (MQ=255)
                                                                                             gcATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGt                                                                               <  2:1586616/125‑1 (MQ=255)
                                                                                                tGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCg                                                   <  1:2092452/150‑1 (MQ=255)
                                                                                                 gCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTc                                                                                  >  2:709860/1‑118 (MQ=255)
                                                                                                 gCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTc                                                                                  <  1:709860/118‑1 (MQ=255)
                                                                                                 ccTCGGCGGCAACGTGCGCCAATCGCCCCGCCGGCACCGCGCGTGCCTCGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGt                                                  <  2:1592838/149‑1 (MQ=255)
                                                                                                   tCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTGCCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGttt                                                >  1:1238433/1‑150 (MQ=255)
                                                                                                   tCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGttt                                                >  1:200293/1‑150 (MQ=255)
                                                                                                         gcAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCggg                                          <  1:199278/150‑1 (MQ=255)
                                                                                                          cAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGc                                         >  2:837036/1‑150 (MQ=255)
                                                                                                           aaCGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGcc                                        >  2:681089/1‑150 (MQ=255)
                                                                                                            aCGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCa                                       <  2:1424976/150‑1 (MQ=255)
                                                                                                                gcgcCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAg                                      <  1:1369295/147‑1 (MQ=255)
                                                                                                                gcgcCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAg                                      >  2:1369295/1‑147 (MQ=255)
                                                                                                                  gcCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAg                                                                                                                >  1:825259/1‑71 (MQ=255)
                                                                                                                  gcCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAg                                                                                                                <  2:825259/71‑1 (MQ=255)
                                                                                                                   ccAATCGACCAGCAGGCCCAGCGCGTGCAACCGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCggg                                <  2:1946624/150‑1 (MQ=255)
                                                                                                                     aaTCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGATCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCa                              >  2:1238492/1‑150 (MQ=255)
                                                                                                                        cGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGc                           <  1:898424/150‑1 (MQ=255)
                                                                                                                          aCCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCgg                         >  2:57628/1‑150 (MQ=255)
                                                                                                                           ccAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGc                                                                                                                                    >  2:613731/1‑42 (MQ=255)
                                                                                                                           ccAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGc                                                                                                                                    <  1:613731/42‑1 (MQ=255)
                                                                                                                            cagcagGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGc                       <  2:1053972/150‑1 (MQ=255)
                                                                                                                            cagcagGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGc                       >  2:2338068/1‑150 (MQ=255)
                                                                                                                             agcagGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAg                                      <  2:640114/134‑1 (MQ=255)
                                                                                                                             agcagGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAg                                      >  1:640114/1‑134 (MQ=255)
                                                                                                                                 ggCACAGCGCGTGCAGCGGCGGGCCGTTGCTCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGTTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCGGTgg                  <  2:493212/150‑1 (MQ=255)
                                                                                                                                 ggCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTgg                  <  2:1030602/150‑1 (MQ=255)
                                                                                                                                 ggCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTgg                  <  1:57628/150‑1 (MQ=255)
                                                                                                                                 ggCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCCGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTgg                  >  2:652121/1‑150 (MQ=255)
                                                                                                                                  gCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTggg                 >  1:248926/1‑150 (MQ=255)
                                                                                                                                    acaGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGcc               >  2:106388/1‑150 (MQ=255)
                                                                                                                                       gcgcGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGTGCCCGc            >  1:1530054/1‑150 (MQ=255)
                                                                                                                                             gCAACGGCTGACCGTTCCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAAt      >  2:1799202/1‑150 (MQ=255)
                                                                                                                                              cAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAATg     >  2:1659765/1‑150 (MQ=255)
                                                                                                                                               aaCGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAATGt    >  1:1621140/1‑150 (MQ=255)
                                                                                                                                                aCGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGTGCCCGCTGAAATGTa   <  2:1530054/150‑1 (MQ=255)
                                                                                                                                                 cGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAATGTAc  <  1:1028898/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TCGCGATGATGTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAATGTAC  >  NC_002947/3296395‑3296689

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