Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,062,315 A→C 5.4% V675G (GTG→GGG)  PP_3581 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,062,3150AC 5.4% 231.0 / 3.6 74V675G (GTG→GGG) PP_3581membrane protein
Reads supporting (aligned to +/- strand):  ref base A (33/36);  new base C (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 = 9.61e-01

CGATGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAACAATCGCT  >  NC_002947/4062173‑4062458
                                                                                                                                              |                                                                                                                                               
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                                                                       <  7:95548/149‑1 (MQ=255)
                                                                    gcgATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGtt                                                                       <  3:223733/149‑1 (MQ=255)
                                                                        tGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCa                                                                                              <  1:301445/122‑1 (MQ=255)
                                                                        tGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCa                                                                                              >  2:301445/1‑122 (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                                                          <  3:34045/119‑1 (MQ=255)
                                                                                                               cTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCgg                                                          >  4:34045/1‑119 (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)
                                                                                                                                              |                                                                                                                                               
CGATGTAGAGCGCATCATCGAGGCCGTCGAGGCGGAAATGGATGCCTTCATACCCGCCTTTACCTTGCGCGATGCAGTCTAGCCCGCAGCAGTGAGTTCCGGGCCCATGGCCTGCAACTTGCTCATGAACTGCTCGGCCGGCACGGGCTGGCCCAGCAGGTAGCCCTGCAACGAGTCACATCCCAGACGGGTCAGGAAGTCCTGCTGCTTGTCGGTTTCAACCCCCTCGGCAACGATACGCAGGCCCAACGCCTGGCCCAGTGCGACGATTGCCGAAACAATCGCT  >  NC_002947/4062173‑4062458

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