Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,062,856 G→T 7.5% G591W (GGG→TGG)  PP_1838 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,062,8560GT 7.5% 208.1 / 6.3 67G591W (GGG→TGG) PP_1838membrane protein
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base T (2/3);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

CTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCCGCGGATGTACCGCTACCAGTT  >  NC_002947/2062709‑2063004
                                                                                                                                                   |                                                                                                                                                    
cTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTgggg                                                                                                                                                     >  4:34204/1‑149 (MQ=255)
    gCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGgcgc                                                                                                                                                 <  3:34204/149‑1 (MQ=255)
     ccGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCCCCAGCACTGGTGgcgag                                                                                                                                                <  3:96402/149‑3 (MQ=255)
       gCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGAc                                                                                                                                              <  6:309925/149‑1 (MQ=255)
        cTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGAcc                                                                                                                                             <  7:62882/149‑1 (MQ=255)
         ttCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCt                                                                                                                                            >  7:179467/1‑149 (MQ=255)
              cGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCGCGCCCCCAGCACTGGTGgcgagcctt                                                                                                                                            >  3:195003/1‑138 (MQ=255)
              cGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCGCGCCCCCAGCACTGGTGgcgagcctt                                                                                                                                            <  4:195003/144‑7 (MQ=255)
                   cGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTg                                                                                                                                  >  6:298032/1‑149 (MQ=255)
                        ctgctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGTATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGc                                                                                                                                          >  3:257054/1‑136 (MQ=255)
                        ctgctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGTATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGc                                                                                                                                          <  4:257054/136‑1 (MQ=255)
                        ctgctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAAcc                                                                                                                                    >  3:252505/1‑142 (MQ=255)
                        ctgctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAAcc                                                                                                                                    <  4:252505/142‑1 (MQ=255)
                           ctgATTGCCCCGGGCATTCAGGAGAATTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGCGTGGCAAGGCTCGCCACCAGTGCTCGGGGCGCGACCTGCTCAACCTGCCCGAggg                                                                                                                          >  1:168347/1‑149 (MQ=255)
                             gATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGt                                                                                                                         >  6:242194/1‑148 (MQ=255)
                             gATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTg                                                                                                                        >  7:235078/1‑149 (MQ=255)
                              aTTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCTGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCTAGGCTGGCCACCAGTGCGGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTg                                                                                                                        <  5:267448/148‑1 (MQ=255)
                                tGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGAc                                                                                                                      >  1:75072/1‑148 (MQ=255)
                                    ccGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAgg                                                                                                                  <  2:286455/148‑1 (MQ=255)
                                    ccGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAgg                                                                                                                  >  1:286455/1‑148 (MQ=255)
                                     cGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGc                                                                                                                >  5:235176/1‑149 (MQ=255)
                                         cATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAg                                                                                                            <  6:235176/149‑1 (MQ=255)
                                         cATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAg                                                                                                            <  8:235078/149‑1 (MQ=255)
                                            tCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAggg                                                                                                                 >  5:145834/1‑141 (MQ=255)
                                            tCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAggg                                                                                                                 <  6:145834/141‑1 (MQ=255)
                                             cAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCa                                                                                                         <  5:298032/148‑1 (MQ=255)
                                              atgagaATTTGGGTGGGGTGAATCAAACCGTGGGTAACCATGTCGTCATAGTGCCGACGATCATTGGCCGCGTGGGTGGCAATTCTCGCTACCATTGCTGGGGTCGAGACCTGCTCAACCTGCCCGAGGGTTACCAGGGCTTAGGCAt                                                                                                        <  2:294756/146‑1 (MQ=255)
                                              aGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCAt                                                                                                        >  4:165072/1‑148 (MQ=255)
                                              aGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCCACATCGTGCCGACAATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCAt                                                                                                        >  6:116885/1‑148 (MQ=255)
                                               ggAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGa                                                                                                      >  4:278486/1‑149 (MQ=255)
                                                gagaAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGAt                                                                                                     >  2:93550/1‑149 (MQ=255)
                                                   aaGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCaa                                                                                                  <  1:93550/149‑1 (MQ=255)
                                                     gTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAg                                                                                                 >  8:273092/1‑148 (MQ=255)
                                                     gTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAg                                                                                                 <  7:273092/148‑1 (MQ=255)
                                                      ttCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGcc                                                                                               <  1:39531/149‑1 (MQ=255)
                                                            gcggtgAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCggg                                                                                         >  3:101685/1‑149 (MQ=255)
                                                             cggtgAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGc                                                                                        <  8:179467/149‑1 (MQ=255)
                                                                      acacCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAg                                                                                  >  7:30577/1‑146 (MQ=255)
                                                                       cacCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAg                                                                                  >  3:222513/1‑145 (MQ=255)
                                                                         ccGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGAtt                                                                            >  2:116158/1‑149 (MQ=255)
                                                                              ggTACCCAGGTCGACATCGTGCCAACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCACGCCGTCGGGCAGCGAGCAGATTGTCgg                                                                       >  7:166114/1‑149 (MQ=255)
                                                                                 aCCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCgggag                                                                    >  1:275431/1‑147 (MQ=255)
                                                                                 aCCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCgg                                                                       <  2:275431/146‑1 (MQ=255)
                                                                                  cccAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTg                                                                   <  3:165072/149‑1 (MQ=255)
                                                                                   ccAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTgg                                                                  >  3:40928/1‑149 (MQ=255)
                                                                                     aGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTg                                                                <  2:75072/149‑1 (MQ=255)
                                                                                          gACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGccc                                                                                                                               <  2:195740/81‑1 (MQ=255)
                                                                                          gACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGccc                                                                                                                               >  1:195740/1‑81 (MQ=255)
                                                                                          gACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAgg                                                            <  8:30577/148‑1 (MQ=255)
                                                                                                 tGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGt                                                    <  4:40928/149‑1 (MQ=255)
                                                                                                           catGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATc                                           >  2:120900/1‑148 (MQ=255)
                                                                                                                gccgccTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTcc                                     <  4:101685/149‑1 (MQ=255)
                                                                                                                 ccgccTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTcc                                     >  5:42760/1‑148 (MQ=255)
                                                                                                                  cgccTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCaa                                   <  3:61350/149‑1 (MQ=255)
                                                                                                                  cgccTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCa                                    >  4:61350/1‑148 (MQ=255)
                                                                                                                    ccTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAg                                  >  7:32393/1‑148 (MQ=255)
                                                                                                                        ggTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGGTGGTGCAGGGTGACCGTATCCTGATCGATTCCAAGGACAt                             <  5:116885/149‑1 (MQ=255)
                                                                                                                        ggTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACAt                             >  2:215035/1‑149 (MQ=255)
                                                                                                                          tGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGa                           <  4:222513/149‑1 (MQ=255)
                                                                                                                                ggCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCcgc                     <  3:278486/149‑1 (MQ=255)
                                                                                                                                 gCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCTACCTGCCCGTGGGTGACCAGGGCGTTGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTGCAAGGACATGAGCCcgc                     <  2:168347/148‑1 (MQ=255)
                                                                                                                                 gCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATc                                                 >  3:174312/1‑120 (MQ=255)
                                                                                                                                 gCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATc                                                 <  4:174312/120‑1 (MQ=255)
                                                                                                                                  cTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCcgcg                    >  8:23418/1‑148 (MQ=255)
                                                                                                                                     gCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCCGCGGATg                >  2:250011/1‑149 (MQ=255)
                                                                                                                                      ccaccaGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGGGACCGTATCCTGATCGAGTCCAAGGACATGAGCCCGCGGATg                >  6:187228/1‑148 (MQ=255)
                                                                                                                                           agcacTGGTGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGa                                             <  3:158421/110‑1 (MQ=255)
                                                                                                                                           agcacTGGTGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGa                                             >  4:158421/5‑114 (MQ=255)
                                                                                                                                                   gggCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCCGCGGATGTACCGCTACCAGtt  <  1:120900/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCCGCGGATGTACCGCTACCAGTT  >  NC_002947/2062709‑2063004

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