Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,483,798 G→C 5.5% P709R (CCG→CGG)  PP_3091 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,483,7980GC 5.5% 272.2 / 3.1 84P709R (CCG→CGG) PP_3091membrane protein
Reads supporting (aligned to +/- strand):  ref base G (38/41);  new base C (2/3);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGGTCGTTGAGGCTCAGCGCCGGGAACTGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATG  >  NC_002947/3483658‑3483946
                                                                                                                                            |                                                                                                                                                    
ccAGGTCGTTGAGGCTCAGCGCCGGGAACTGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAacgac                                                                                                                                              <  8:321122/149‑1 (MQ=255)
  aGGTCGTTGAGGCTCAGCGCCGGGAACTGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACgg                                                                                                                                            >  3:68013/1‑149 (MQ=255)
  aGGGCGGTGAGGCTCAGCGCCGGGAACGGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCATCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGAGCAGCTGTGCGTTGCCGCGCAacgaaca                                                                                                                                            <  8:303577/149‑4 (MQ=255)
   ggTCGTTGAGGCTCAGCGCCGGGAACTGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGATATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGGTCGTTGCCGCGCAacgaacag                                                                                                                                           <  1:28712/149‑5 (MQ=255)
        ttGAGGCTCAGCGCCGGGAACTGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTg                                                                                                                                      >  8:12509/1‑149 (MQ=255)
             gCTCAGCGCCGGGAACTGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCt                                                                                                                                 >  1:157458/1‑149 (MQ=255)
             gCTCAGCGCCGGGAACTGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCt                                                                                                                                 >  5:255771/1‑149 (MQ=255)
               tCAGCGCCGGGAACTGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGACCGCCTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGc                                                                                                                               >  7:130248/1‑149 (MQ=255)
               tCAGCGCCGGGAACTGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGc                                                                                                                               >  7:363805/1‑149 (MQ=255)
               gcAGCCCCGGCAACTGCCGGCTCGACTCGAGCTTCAGGCTCATGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAAAGACTTTAGGTTCTGCCGGGGTTCGTTGTTCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGc                                                                                                                               <  3:39687/148‑1 (MQ=255)
                 aGCGCCGGGAACTGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGc                                                                                                                               <  4:68013/147‑1 (MQ=255)
                 aGCGCCGGGAACTGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCt                                                                                                                              <  4:155282/148‑1 (MQ=255)
                  gcgcCGGGAACTGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTc                                                                                                                             >  1:338506/1‑148 (MQ=255)
                      cGGGAACTGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGACAACGACGGCGCCTGGCCCTGCTCCAGc                                                                                                                         >  3:390277/1‑148 (MQ=255)
                       gggAACTGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCaa                                                                                                                       <  2:49075/149‑1 (MQ=255)
                        ggACCTGCCGGCTGGACTCGCGCTTCAGTCTGCTTTATGCTCGTTCGACCTGATATGTGGTGATAATGTACTAGTTGATGTTGTGCGGTGTTTCGGTGTTCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAAc                                                                                                                      <  8:212812/149‑1 (MQ=255)
                           aCTGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGc                                                                                                                   <  8:107660/149‑1 (MQ=255)
                           aCTGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGACAACGACGGCGCCTGGCCCTGCTCCAGCAACTGc                                                                                                                   <  4:390277/149‑1 (MQ=255)
                                      gactcgaGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCTGCAACTGCACGTACAGCg                                                                                                         <  8:363805/148‑1 (MQ=255)
                                      gactcgaGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTAGAGCgg                                                                                                        <  6:211332/149‑1 (MQ=255)
                                       actcgaGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAacgacg                                                                                                                                             <  2:224168/111‑1 (MQ=255)
                                       actcgaGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAacgacg                                                                                                                                             >  1:224168/1‑111 (MQ=255)
                                       actcgaGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCggg                                                                                                       <  4:353251/149‑1 (MQ=255)
                                        ctcgaGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCgggg                                                                                                      >  4:184355/1‑149 (MQ=255)
                                              cTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCAtttt                                                                                                <  4:79449/149‑1 (MQ=255)
                                               tGCAGGCACAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTc                                                                                               >  8:175323/1‑149 (MQ=255)
                                                gCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTCGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCCGCAACGACGGCGCCTGGCCCTGCACCAGCAACTCCACGTACAGCGCGGCATATTCt                                                                                              >  8:204233/1‑149 (MQ=255)
                                                    gcgcAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGa                                                                                          >  4:34284/1‑149 (MQ=255)
                                                      gcAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATc                                                                                        <  7:12509/149‑1 (MQ=255)
                                                               cgcgATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGc                                                                                <  3:34284/148‑1 (MQ=255)
                                                               cgcgATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCg                                                                               >  3:372978/1‑149 (MQ=255)
                                                               agcgATCGACGAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATGAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGc                                                                                <  8:130248/147‑1 (MQ=255)
                                                                      gACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCg                                                                        >  4:77312/1‑149 (MQ=255)
                                                                       aCCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCgg                                                                       >  1:96431/1‑149 (MQ=255)
                                                                        ccAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTATAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCggg                                                                      <  2:157458/149‑1 (MQ=255)
                                                                        ccAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCggg                                                                      <  1:182846/149‑1 (MQ=255)
                                                                           gCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCggg                                                                                                       <  4:214774/113‑1 (MQ=255)
                                                                           gCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCggg                                                                                                       >  3:214774/1‑113 (MQ=255)
                                                                                cTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTc                                                              <  2:338506/149‑1 (MQ=255)
                                                                                 tGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCg                                                             >  7:22570/1‑149 (MQ=255)
                                                                                   gagaTCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTTTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGc                                                            <  1:265773/148‑1 (MQ=255)
                                                                                    agttgaTGCACGACTTGCGGTTCTGCCGGGGTTCGTTGATCAGCTGTTCGGTGGGCGGCAACGTCGGCGCCTGGCCCTGGTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCt                                                          <  7:204233/144‑1 (MQ=255)
                                                                                     gttcttGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTTGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTg                                                         <  3:138864/144‑1 (MQ=255)
                                                                                      aTCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCa                                                                                  <  7:126110/123‑1 (MQ=255)
                                                                                      aTCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCa                                                                                  >  8:126110/1‑123 (MQ=255)
                                                                                       tCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTgcg                                                       <  2:211917/149‑1 (MQ=255)
                                                                                       tCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTgcg                                                       <  4:356481/149‑1 (MQ=255)
                                                                                           gAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTTGACCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTg                                                         >  4:138864/1‑143 (MQ=255)
                                                                                           gAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTc                                                              >  1:303628/1‑138 (MQ=255)
                                                                                           gAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTc                                                              <  2:303628/138‑1 (MQ=255)
                                                                                           gAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGgcgc                                                   >  6:46746/1‑149 (MQ=255)
                                                                                           gAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGgcgc                                                   >  5:136830/1‑149 (MQ=255)
                                                                                            aaCGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCt                                                  >  3:101833/1‑149 (MQ=255)
                                                                                                 cTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGccc                                             <  4:372978/149‑1 (MQ=255)
                                                                                                 cTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGccc                                             <  6:136830/149‑1 (MQ=255)
                                                                                                    gAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAgg                                          >  4:375641/1‑149 (MQ=255)
                                                                                                    gAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAg                                           <  2:96431/148‑1 (MQ=255)
                                                                                                      ggTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGc                                               <  5:328395/142‑1 (MQ=255)
                                                                                                      ggTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGc                                               >  6:328395/1‑142 (MQ=255)
                                                                                                       gTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCt                                       >  3:194451/1‑149 (MQ=255)
                                                                                                        ttCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTg                                      <  1:357148/149‑1 (MQ=255)
                                                                                                        ttCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTg                                      <  7:175323/149‑1 (MQ=255)
                                                                                                        ttCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTg                                      <  6:255771/149‑1 (MQ=255)
                                                                                                                  tttCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTg                            <  3:184355/149‑1 (MQ=255)
                                                                                                                     cGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTgg                           >  1:39887/1‑147 (MQ=255)
                                                                                                                      gTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAg                        >  7:324821/1‑149 (MQ=255)
                                                                                                                      gTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAg                        >  3:52683/1‑149 (MQ=255)
                                                                                                                       ttGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGc                                                                                                                               <  8:259068/45‑1 (MQ=255)
                                                                                                                       ttGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGc                                                                                                                               >  7:259068/1‑45 (MQ=255)
                                                                                                                       ttGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGc                       >  3:296279/1‑149 (MQ=255)
                                                                                                                        tGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGc                       >  3:220398/1‑148 (MQ=255)
                                                                                                                           tCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGc                                                                          <  2:292144/94‑1 (MQ=255)
                                                                                                                           tCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGc                                                                          >  1:292144/1‑94 (MQ=255)
                                                                                                                           tCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGcgac                    >  5:29805/1‑148 (MQ=255)
                                                                                                                              gCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGcga                     <  2:203952/144‑1 (MQ=255)
                                                                                                                              gCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGcga                     >  1:203952/1‑144 (MQ=255)
                                                                                                                               gccgtCGTTGCCGCGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAgg               >  1:134161/5‑149 (MQ=255)
                                                                                                                               cTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGcgac                    <  2:149101/144‑1 (MQ=255)
                                                                                                                               cTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGcgac                    >  1:149101/1‑144 (MQ=255)
                                                                                                                                tGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGc              >  3:298880/1‑149 (MQ=255)
                                                                                                                                   tCGTTGCCGCGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCgg           <  2:134161/149‑1 (MQ=255)
                                                                                                                                    cGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCgg           <  4:101833/148‑1 (MQ=255)
                                                                                                                                     gTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGt         >  2:301052/1‑149 (MQ=255)
                                                                                                                                      ttGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTc        <  4:220398/149‑1 (MQ=255)
                                                                                                                                          gcgGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCAGCGGCGCAGGCCCAGGCCTGGGGCACGCTGGCAGCGCCGAAGGCCGGGTCGGGa    >  8:287624/1‑149 (MQ=255)
                                                                                                                                           cgGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGAt   <  8:22570/149‑1 (MQ=255)
                                                                                                                                            ggCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATg  <  8:324821/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
CCAGGTCGTTGAGGCTCAGCGCCGGGAACTGCCGGCTCGACTCGAGCTGCAGGCGCAGGTATTCGCGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATG  >  NC_002947/3483658‑3483946

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: