Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,062,318 G→T 5.4% P674H (CCC→CAC)  PP_3581 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,062,3180GT 5.4% 239.5 / 3.5 74P674H (CCC→CAC) PP_3581membrane protein
Reads supporting (aligned to +/- strand):  ref base G (34/36);  new base T (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

CGATGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAACAATCGCTGCATCGT  >  NC_002947/4062173‑4062465
                                                                                                                                                 |                                                                                                                                                   
cgATGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCGCATGAACTGCTCGGCCGGCACGGGCt                                                                                                                                                  >  7:56935/1‑149 (MQ=255)
 gATGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTGGGCCAGCCCGTGccg                                                                                                                                                 >  3:174955/1‑147 (MQ=255)
  aTGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTgg                                                                                                                                                <  3:265848/149‑1 (MQ=255)
       gagCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCcag                                                                                                                                           >  4:194403/1‑149 (MQ=255)
       gagCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCAGCCCGTGCCGGccga                                                                                                                                            <  8:4157/148‑3 (MQ=255)
         gcgcATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCcagca                                                                                                                                         <  3:194403/149‑1 (MQ=255)
            catcatCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTGGGCCAGCCCGTGCCGGCCgagcagtt                                                                                                                                      <  4:174955/149‑3 (MQ=255)
                atcGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGcc                                                                                                                                  <  8:12140/149‑1 (MQ=255)
                 tcGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGccc                                                                                                                                 >  5:126699/1‑149 (MQ=255)
                      gCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCaa                                                                                                                            >  4:223733/1‑149 (MQ=255)
                        cGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCAGGGCCTGCAACGTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAAc                                                                                                                           >  1:120363/1‑148 (MQ=255)
                        cGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCACGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCACCg                                                                                                                          >  2:10252/1‑149 (MQ=255)
                         gTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGa                                                                                                                         <  2:272342/149‑1 (MQ=255)
                         gTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGa                                                                                                                         <  1:359065/149‑1 (MQ=255)
                            gAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTc                                                                                                                      <  1:318346/149‑1 (MQ=255)
                                  gAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGTCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATcc                                                                                                                >  4:99624/1‑149 (MQ=255)
                                  gAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATcc                                                                                                                <  2:197097/149‑1 (MQ=255)
                                     aTGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGa                                                                                                                         >  6:293145/1‑137 (MQ=255)
                                            ccTTCTAACCCGCTTTTACGTTGCGGCATCATGTCTCGCCGGCTCCATGCAGTACCGGGCCCATGGCCTGAAACTTGCTATTTACTTGCTCGGCCGGAAGGGGCTGGCCAAGC‑TGTATGCCTGCAACGAGTCACATCCCAGACGGGTCa                                                                                                     <  6:141863/149‑1 (MQ=255)
                                              ttCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGG‑CCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAg                                                                                                    >  3:250382/1‑148 (MQ=255)
                                                cATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGTCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAg                                                                                                    <  8:68782/147‑1 (MQ=255)
                                                cATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAg                                                                                                    <  8:73757/147‑1 (MQ=255)
                                                   aCCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGACGGCGCGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAg                                                                                                >  7:88035/1‑148 (MQ=255)
                                                    cccGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTc                                                                                              <  5:322599/149‑1 (MQ=255)
                                                     ccGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTcc                                                                                             >  5:101318/1‑149 (MQ=255)
                                                         cTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCctgct                                                                                         >  3:361083/1‑149 (MQ=255)
                                                           ttACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCtgctgc                                                                                       >  1:185663/1‑149 (MQ=255)
                                                            tACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCtgctgc                                                                                       <  7:291909/148‑1 (MQ=255)
                                                             aCCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCtt                                                                                     >  8:95548/1‑149 (MQ=255)
                                                               cTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGt                                                                                   <  6:126699/149‑1 (MQ=255)
                                                                    gcgATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGtt                                                                              <  3:223733/149‑1 (MQ=255)
                                                                    gcgATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGtt                                                                              <  7:95548/149‑1 (MQ=255)
                                                                        tGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCa                                                                                                     >  2:301445/1‑122 (MQ=255)
                                                                        tGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCa                                                                                                     <  1:301445/122‑1 (MQ=255)
                                                                          cAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAAc                                                                         <  2:185663/148‑1 (MQ=255)
                                                                           aGTCTCGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGCTTCAAccc                                                                       <  2:120363/149‑1 (MQ=255)
                                                                           aGTCTAGCCCGCAGCAGTGAGTTACGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTAGGTTTCAAccc                                                                       >  3:80350/1‑149 (MQ=255)
                                                                            gTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAAcccc                                                                      >  1:334871/1‑149 (MQ=255)
                                                                             tCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCCCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCtgctgc                                                                                       >  5:267977/1‑131 (MQ=255)
                                                                             tCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCCCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCtgctgc                                                                                       <  6:267977/131‑1 (MQ=255)
                                                                              cTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGa                                                                                                  >  7:277069/1‑119 (MQ=255)
                                                                              cTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGa                                                                                                  <  8:277069/119‑1 (MQ=255)
                                                                                 gcgcGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGG‑CCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTc                                                                   <  4:250382/143‑1 (MQ=255)
                                                                                         cagTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATa                                                         >  3:199041/1‑149 (MQ=255)
                                                                                            tGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATCAACTGCTCGGCCGGCACGGCCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACg                                                       >  8:207188/1‑148 (MQ=255)
                                                                                              aGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCa                                                     <  1:10252/148‑1 (MQ=255)
                                                                                                 tCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGcc                                                 <  6:101318/149‑1 (MQ=255)
                                                                                                    gggCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCaa                                              >  3:191867/1‑149 (MQ=255)
                                                                                                       cccATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGc                                           >  6:43910/1‑149 (MQ=255)
                                                                                                           tGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTgg                                       <  4:361083/149‑1 (MQ=255)
                                                                                                               cTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCgg                                                                 >  4:34045/1‑119 (MQ=255)
                                                                                                               cTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCgg                                                                 <  3:34045/119‑1 (MQ=255)
                                                                                                               cTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGcc                                     >  7:325776/1‑147 (MQ=255)
                                                                                                               cTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCa                                   <  2:334729/149‑1 (MQ=255)
                                                                                                                    aCTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGc                               >  5:310678/1‑148 (MQ=255)
                                                                                                                      ttGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGcgac                            >  2:300572/1‑149 (MQ=255)
                                                                                                                            attaaCTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGc                      <  7:207188/146‑1 (MQ=255)
                                                                                                                            aTGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGAtt                        <  2:334871/147‑1 (MQ=255)
                                                                                                                            aTGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGc                      >  7:87119/1‑149 (MQ=255)
                                                                                                                            aTGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGACACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGAtt                        <  1:300572/147‑1 (MQ=255)
                                                                                                                             tGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGt                                                                                                       >  6:344666/1‑67 (MQ=255)
                                                                                                                             tGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGt                                                                                                       <  5:344666/67‑1 (MQ=255)
                                                                                                                               aaCTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCg                    >  3:162405/1‑148 (MQ=255)
                                                                                                                                 cTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGc                                                  <  7:144856/116‑1 (MQ=255)
                                                                                                                                 cTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGc                                                  >  8:144856/1‑116 (MQ=255)
                                                                                                                                  tGCTGGGCCAGCCCGTGCCGGCCGAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGaaa                 >  6:285971/1‑148 (MQ=255)
                                                                                                                                  tGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGaaa                 >  8:276950/1‑148 (MQ=255)
                                                                                                                                    cTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAAc                <  8:325776/147‑1 (MQ=255)
                                                                                                                                     tCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAACaa              >  4:107885/1‑148 (MQ=255)
                                                                                                                                      cggccggcACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAACAATc            <  8:87119/149‑1 (MQ=255)
                                                                                                                                       ggccggcACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAACAATc            <  6:310678/148‑1 (MQ=255)
                                                                                                                                        gccggcACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGaa                  >  2:13324/1‑141 (MQ=255)
                                                                                                                                        gccggcACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGaa                  <  1:13324/141‑1 (MQ=255)
                                                                                                                                         ccggcACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAACAATCGCt         <  4:162405/149‑1 (MQ=255)
                                                                                                                                                 ggCTGGCCCAGCATGTTGGCCTGCAACGTTTAACATCTGAGACTTGTCATGAAGTTCTGCTGCTTGTCGGTTTCAATCCCCTGGGAAACTATACGCTGGGCCTACGTGTGGCCCAGTGCGACGATTGGTGAAACAATCTCTGCATGGt  <  8:88035/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGATGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAACAATCGCTGCATCGT  >  NC_002947/4062173‑4062465

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