Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,296,545 G→A 7.3% Q125* (CAG→TAG) ‡ PP_2898 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,296,5450GA 7.3% 274.2 / 8.0 95Q125* (CAG→TAG) ‡PP_2898hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (46/42);  new base A (2/5);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

GATGATGTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAATGTACCGGTG  >  NC_002947/3296399‑3296694
                                                                                                                                                  |                                                                                                                                                     
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:2200612/1‑150 (MQ=255)
          aGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTgg                                                                                                                                          >  1:2183948/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                                                                                                                                               >  1:177928/1‑121 (MQ=255)
                                  cACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCCAGCCGTTgc                                                                                                                                               <  2:177928/121‑1 (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                                                                           >  1:2312339/1‑150 (MQ=255)
                                                                         gCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGcc                                                                           >  2:351854/1‑150 (MQ=255)
                                                                         gCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGcc                                                                           <  2:1008309/150‑1 (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                                                                                       <  1:709860/118‑1 (MQ=255)
                                                                                             gCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTc                                                                                       >  2:709860/1‑118 (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                                           >  2:1369295/1‑147 (MQ=255)
                                                                                                            gcgcCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAg                                           <  1:1369295/147‑1 (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                                                                                                                                         <  1:613731/42‑1 (MQ=255)
                                                                                                                       ccAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGc                                                                                                                                         >  2:613731/1‑42 (MQ=255)
                                                                                                                        cagcagGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGc                            >  2:2338068/1‑150 (MQ=255)
                                                                                                                        cagcagGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGc                            <  2:1053972/150‑1 (MQ=255)
                                                                                                                         agcagGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAg                                           >  1:640114/1‑134 (MQ=255)
                                                                                                                         agcagGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAg                                           <  2:640114/134‑1 (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)
                                                                                                                                                  ggCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAATGTACCGGTg  <  2:182410/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GATGATGTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAATGTACCGGTG  >  NC_002947/3296399‑3296694

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: