Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,706,338 G→C 7.3% A135P (GCC→CCC) ‡ PP_1500 → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,706,3380GC 7.3% 206.4 / 7.1 69A135P (GCC→CCC) ‡PP_1500lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base C (3/2);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.91e-01

GGCGGTAACCTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATGGTGTGCCGG  >  NC_002947/1706191‑1706485
                                                                                                                                                   |                                                                                                                                                    
ggcggTAACCTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGc                                                                                                                                                     >  6:330115/1‑149 (MQ=255)
         cTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTgca                                                                                                                                            <  5:230680/149‑1 (MQ=255)
             cAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGcagcag                                                                                                                                        >  3:269010/1‑149 (MQ=255)
               gggCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGgc                                                                                                                                      <  3:96801/149‑1 (MQ=255)
                   aCCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGccc                                                                                                                                  <  2:300738/149‑1 (MQ=255)
                              aGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAggg                                                                                                                       >  3:90140/1‑149 (MQ=255)
                                   cTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAgg                                                                                                                        <  4:200369/143‑1 (MQ=255)
                                   cTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAgg                                                                                                                        >  3:200369/1‑143 (MQ=255)
                                       cTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCt                                                                                                              >  8:124588/1‑149 (MQ=255)
                                          ccGGAGCAGGGCTACCAGGTGCTGATGAGCGATTACCGCGGCTATGGCTTGGCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTgg                                                                                                            <  8:299631/148‑1 (MQ=255)
                                            ggAGCATGGCTACCAGGGGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCCGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGTCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTAc                                                                                                         <  3:204135/149‑1 (MQ=255)
                                            ggAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTAc                                                                                                         >  2:8355/1‑149 (MQ=255)
                                               gCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTActgc                                                                                                      <  1:8355/149‑1 (MQ=255)
                                                 aGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTActgctg                                                                                                    <  1:180172/149‑1 (MQ=255)
                                                  gggCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTgg                                                                                                   <  5:330115/149‑1 (MQ=255)
                                                  gggCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTgg                                                                                                   <  4:269010/149‑1 (MQ=255)
                                                   ggCTACCAGGTACTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAgg                                                                                                                                       <  3:65408/112‑1 (MQ=255)
                                                   ggCTACCAGGTACTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAgg                                                                                                                                       >  4:65408/1‑112 (MQ=255)
                                                    gCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTCCACCAGGCGCCCCAGGTCAAGGGCAAGCCGCTGGTACTGCTgggg                                                                                                 >  5:199253/1‑149 (MQ=255)
                                                     caacgaGGTGCTGTAGATCGATTACCGCGGGTATTGGTTGGTTCCGGGCCATCCGATCCTTGCGGAGGTTTACCAGGACATTGCTTCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCTAGCCGCTGGTACTGCTGGGGc                                                                                                <  6:199253/144‑1 (MQ=255)
                                                     cTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCAGGCCATGGCTgcagcaat                                                                                                                                       >  3:328406/1‑108 (MQ=39)
                                                     cTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCAGGCCATGGCTgcagcaat                                                                                                                                       <  4:328406/110‑3 (MQ=39)
                                                        ccAGGTGCTGATTCTCGATTACCGCGGCTATGGCTGGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCAGGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCaga                                                                                             <  4:45534/149‑1 (MQ=255)
                                                        ccAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCa                                                                                               >  6:338880/1‑147 (MQ=255)
                                                               cTGAAGATCGATTCCCGCGGCTATGGCTTTTTTCACGGACAGCCGAGCCTGCCTGAGGTGTACCATGACATTGCTGCAGCCATTTCCTGTCTGCATCAGGCTCACGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGACAGAGCCTgg                                                                                       <  8:145470/148‑1 (MQ=255)
                                                                   tgatCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCt                                                                                                                                                   >  6:281180/1‑84 (MQ=255)
                                                                   tgatCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCt                                                                                                                                                   <  5:281180/84‑1 (MQ=255)
                                                                    gatCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTg                                                                                >  3:250109/1‑149 (MQ=255)
                                                                      tCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGc                                                                               >  7:146763/1‑148 (MQ=255)
                                                                        gATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCAt                                                                            <  5:338880/149‑1 (MQ=255)
                                                                             ccGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCca                                                                       >  8:176109/1‑149 (MQ=255)
                                                                             ccGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCca                                                                       >  7:146922/1‑149 (MQ=255)
                                                                              cgcgGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCcat                                                                      >  7:170171/1‑149 (MQ=255)
                                                                                cgGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCcattc                                                                    <  7:176109/149‑1 (MQ=255)
                                                                                     aTGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTAc                                                               >  7:41601/1‑149 (MQ=255)
                                                                                      tGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTAc                                                               <  8:146763/148‑1 (MQ=255)
                                                                                        gCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGCGCCTGGGCGGTGCCATGGCCATTCATTACCTg                                                            >  8:265956/1‑149 (MQ=255)
                                                                                        gCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTg                                                            >  6:43259/1‑149 (MQ=255)
                                                                                                 aGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCAc                                                   >  4:243222/1‑149 (MQ=255)
                                                                                                   ggCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCAccc                                                 >  6:270841/1‑149 (MQ=255)
                                                                                                   ggCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCAccc                                                 >  7:133161/1‑149 (MQ=255)
                                                                                                    gCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAGGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCGTGGCCATTCATTACCTGGCGGCCCGcccc                                                >  8:215691/1‑149 (MQ=255)
                                                                                                     ccAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCAcccc                                                >  2:239951/1‑148 (MQ=255)
                                                                                                     ccAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCg                                               <  5:43259/149‑1 (MQ=255)
                                                                                                        gCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGag                                             >  8:259865/1‑148 (MQ=255)
                                                                                                           gAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTATCTGGCGGCCCACCCCGagcagc                                         <  7:215691/149‑1 (MQ=255)
                                                                                                            aGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGagcagc                                         <  1:239951/149‑1 (MQ=255)
                                                                                                               cTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCa                                     <  7:124588/149‑1 (MQ=255)
                                                                                                                 gCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAgc                                   <  4:90140/149‑1 (MQ=255)
                                                                                                                  ccGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAg                                    <  5:270841/147‑1 (MQ=255)
                                                                                                                   cGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCcattcatt                                                                 <  8:321855/117‑1 (MQ=255)
                                                                                                                   cGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCcattcatt                                                                 >  7:321855/1‑117 (MQ=255)
                                                                                                                    ggAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCt                                >  6:92655/1‑149 (MQ=255)
                                                                                                                    ggAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCt                                >  7:67414/1‑149 (MQ=255)
                                                                                                                      aGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTc                              <  4:250109/149‑1 (MQ=255)
                                                                                                                         tgtACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAgc                           <  7:323889/149‑1 (MQ=255)
                                                                                                                             ccAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGccc                       >  4:101201/1‑149 (MQ=255)
                                                                                                                               aGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTg                     <  7:259865/149‑1 (MQ=255)
                                                                                                                                 gACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGgt                   <  5:92655/149‑1 (MQ=255)
                                                                                                                                    aTTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGtt                >  3:238001/1‑149 (MQ=255)
                                                                                                                                     ctgctgcAGCCAGGCCATGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTc               >  7:218472/2‑149 (MQ=255)
                                                                                                                                        ctgcAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGAt            <  7:237889/149‑1 (MQ=255)
                                                                                                                                         tgcAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATg           >  2:33215/1‑149 (MQ=255)
                                                                                                                                          gcAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATgg          >  2:180980/1‑149 (MQ=255)
                                                                                                                                          gcAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATgg          >  1:174956/1‑149 (MQ=255)
                                                                                                                                           cAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATg           <  8:41601/147‑1 (MQ=255)
                                                                                                                                           cAGCCAGGCCATGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTgg                                                                                                   >  7:201952/1‑60 (MQ=37)
                                                                                                                                           cAGCCAGGCCATGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTgg                                                                                                   <  8:201952/60‑1 (MQ=37)
                                                                                                                                              ccATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATGGtgtg      >  5:332981/1‑149 (MQ=255)
                                                                                                                                                   gCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATGGTGTGCCgg  <  1:33215/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGCGGTAACCTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATGGTGTGCCGG  >  NC_002947/1706191‑1706485

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: