Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 758,336 C→A 5.7% intergenic (+4/+24) PP_0649 → / ← PP_0650 MutT/nudix family protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947758,3360CA 5.7% 221.4 / 3.2 70intergenic (+4/+24)PP_0649/PP_0650MutT/nudix family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/34);  new base A (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

GTGCATTTTGCAGAAGATGATGGAGGCCATGAATTCAGATTCTTTTGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGATCATGCCCATTCATGCGAGCCGATCA  >  NC_002947/758195‑758482
                                                                                                                                             |                                                                                                                                                  
gTGCATTTTGCAGAAGATGATGGAGGCCATGAATTCAGATTCTTTTGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACa                                                                                                                                             <  4:119685/149‑1 (MQ=255)
   cATTTTGCAGAAGATGATGGAGGCCATGAATTCAGATTCTTTTGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGc                                                                                                                                           <  4:167748/148‑1 (MQ=255)
          gtgaagaTGATGGAGGCCATGAATTCAGAGTCTTTTGGGATCCTCTCCTGAGCGAGCCCGCAGAAAATTGGCATCAGGTGTTCAAGGATGCTGTGAGCTTTATTCGGAAAAGGCTCACTGATGCCTAAGTGCGCTTACAGCTTTCCCTg                                                                                                                                   <  7:283366/147‑1 (MQ=255)
          cAGAAGATGATGGAGGCCATGAATTCAGATTCTTTTGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGATCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGTTTTCCCt                                                                                                                                    <  2:20552/148‑1 (MQ=255)
          cAGAAGATGATGGAGGCCATGAATTCAGATTCTTTTGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTg                                                                                                                                   >  2:360162/1‑149 (MQ=255)
           agaagaTGATGGAGGCCATGAATTCAGATTCTTTTGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTGAAGGATGGTTTGAGCTGTATTCGGAAACGGCTCACTGTAGCCGAAGTGCGCTTACAGCTTTCACTgc                                                                                                                                  >  4:225456/1‑148 (MQ=255)
           agaagaTGATGGAGGCCATGAATTCAGATTCTTTTGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGAGCTTACAGCTTTCCCTg                                                                                                                                   <  1:146368/148‑1 (MQ=255)
           agaagaTGATGGAGGCCATGAATTCAGATTCTTTTGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGAGCTTACAGCTTTCCCTg                                                                                                                                   <  1:151975/148‑1 (MQ=255)
           agaagaTGATGGAGGCCATGAATTCAGATTCTTTTGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGAGCTTACAGCTTTCCCTg                                                                                                                                   >  2:146368/1‑148 (MQ=255)
           agaagaTGATGGAGGCCATGAATTCAGATTCTTTTGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGAGCTTACAGCTTTCCCTg                                                                                                                                   >  2:151975/1‑148 (MQ=255)
                  gatgGAGGCCATGAATTCAGATTCTTTTGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGAATCAGGTTTTCTAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAATCGTAAGTGCGCTTACAGCTTTCCCTGGGGTAGt                                                                                                                            >  6:333086/1‑148 (MQ=255)
                   atgGAGGCCATGAATTCAGATTCTTTTGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTc                                                                                                                           <  4:11761/148‑1 (MQ=255)
                   atgGAGGCCATGAATTCAGATTCTTTTGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTc                                                                                                                           <  2:307205/148‑1 (MQ=255)
                             tGAATTCAGATTCTTTTGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATc                                                                                                                <  1:244372/149‑1 (MQ=255)
                                    aGATTCTTTTGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGAt                                                                                                         >  6:86014/1‑149 (MQ=255)
                                       ttCTTTTGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGTCATCATGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGAtgtg                                                                                                      <  5:159374/149‑1 (MQ=255)
                                           tttGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTgg                                                                                                                                  <  5:35768/117‑1 (MQ=255)
                                           tttGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTgg                                                                                                                                  >  6:35768/1‑117 (MQ=255)
                                               gCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGtagtag                                                                                               <  1:360162/148‑1 (MQ=255)
                                                 aTCCTCTCATGAGCGAGCCCTCTGAAAATTGGAATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTg                                                                                            >  1:336619/1‑149 (MQ=255)
                                                   ccTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGcc                                                                                          >  6:190139/1‑149 (MQ=255)
                                                    ctctCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGccc                                                                                         <  2:31246/149‑1 (MQ=255)
                                                     tctcATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGccc                                                                                         >  4:90037/1‑148 (MQ=255)
                                                         aTGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGtagtag                                                                                               >  1:100846/1‑138 (MQ=255)
                                                         aTGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGtagtag                                                                                               <  2:100846/138‑1 (MQ=255)
                                                           gagcgagcCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTggg                                                                                                                                 >  2:209296/1‑102 (MQ=255)
                                                           gagcgagcCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTggg                                                                                                                                 <  1:209296/102‑1 (MQ=255)
                                                             gcgagcCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGgta                                                                                                   >  5:73673/1‑130 (MQ=255)
                                                             gcgagcCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGgta                                                                                                   <  6:73673/130‑1 (MQ=255)
                                                                     tctGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGgta                                                                                                   <  7:260815/122‑1 (MQ=255)
                                                                     tctGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGgta                                                                                                   >  8:260815/1‑122 (MQ=255)
                                                                       tGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTgg                                                                                                                                  >  7:91765/1‑89 (MQ=255)
                                                                       tGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTgg                                                                                                                                  <  8:91765/89‑1 (MQ=255)
                                                                       tGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCa                                                                                                                          <  8:107396/97‑1 (MQ=255)
                                                                       tGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCACTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGc                                                                      >  3:102978/1‑149 (MQ=255)
                                                                         aaaaTTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTg                                                                       <  3:90037/146‑1 (MQ=255)
                                                                            aTTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCg                                                                  >  6:232879/1‑148 (MQ=255)
                                                                                 cATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGgt                                                                                                    <  2:219886/109‑1 (MQ=255)
                                                                                 cATCAGGTTTTCAAGGATGCTTTGAGCTTTAATCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGgt                                                                                                    >  1:219886/1‑109 (MQ=255)
                                                                                  aTCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGgtagta                                                                                                <  6:404/112‑1 (MQ=255)
                                                                                  aTCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGgtagta                                                                                                >  5:404/1‑112 (MQ=255)
                                                                                   tCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTg                                                                                                           <  8:243512/100‑1 (MQ=255)
                                                                                   tCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTg                                                                                                           >  7:243512/1‑100 (MQ=255)
                                                                                   tCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGATGCcgc                                                          <  2:336619/149‑1 (MQ=255)
                                                                                   tCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCcgc                                                          >  2:192633/1‑149 (MQ=255)
                                                                                                  tGCTTTCAGCTTTATTCGGAAAAGCCTCACTGAAGCCTAAGAGCGCTTACAGCTTTCCCTGGGGTACTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTg                                           <  3:225456/149‑1 (MQ=255)
                                                                                                     tttGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGc                                                                                           >  4:128038/1‑98 (MQ=255)
                                                                                                         aGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAgg                                    >  5:201942/1‑149 (MQ=255)
                                                                                                         aGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAgg                                    <  5:190139/149‑1 (MQ=255)
                                                                                                              tATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATAc                               >  4:207930/1‑149 (MQ=255)
                                                                                                                ttCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACg                              >  3:66368/1‑148 (MQ=255)
                                                                                                                  cGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGa                            <  5:232879/148‑1 (MQ=255)
                                                                                                                   ggAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGAt                           >  1:324307/1‑148 (MQ=255)
                                                                                                                   ggAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGAt                           <  1:192633/148‑1 (MQ=255)
                                                                                                                    gAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTgg                                                                                                     >  8:144699/1‑73 (MQ=255)
                                                                                                                    gAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTgg                                                                                                     <  7:144699/73‑1 (MQ=255)
                                                                                                                     aaaaGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGATCAt                        <  6:201942/149‑1 (MQ=255)
                                                                                                                       aaGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGATCATGc                      <  3:207930/149‑1 (MQ=255)
                                                                                                                          gCTCACTGAAGCCTAAGTGCGCTTACAGCCTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGATCATGCCCa                   >  8:335137/1‑149 (MQ=255)
                                                                                                                              aCTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGATCATGCCCATTCa               >  5:47568/1‑149 (MQ=255)
                                                                                                                                tGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGATCATGCCCATTc                <  6:185314/146‑1 (MQ=255)
                                                                                                                                tGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGATCATGCCCATTc                >  5:185314/1‑146 (MQ=255)
                                                                                                                                    gCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGATCATGCCCAt                  >  2:170269/1‑140 (MQ=255)
                                                                                                                                    gCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGATCATGCCCAt                  <  1:170269/140‑1 (MQ=255)
                                                                                                                                       tAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTTCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGATCATGCCCATTCATGCGAGcc       >  2:106052/1‑148 (MQ=255)
                                                                                                                                       tAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGATCATGCCCATTCATGCGAGcc       >  6:166193/1‑148 (MQ=255)
                                                                                                                                         aGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGATCATGc                      <  6:342003/131‑1 (MQ=255)
                                                                                                                                         aGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGATCATGc                      >  5:342003/1‑131 (MQ=255)
                                                                                                                                          gTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGATCATGCCCATTCATGCGAGCCGATc   <  5:86014/149‑1 (MQ=255)
                                                                                                                                           tGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTTCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGATCATGCCCATTCATGCGAGCCGATCa  <  1:106052/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GTGCATTTTGCAGAAGATGATGGAGGCCATGAATTCAGATTCTTTTGGCATCCTCTCATGAGCGAGCCCTCTGAAAATTGGCATCAGGTTTTCAAGGATGCTTTGAGCTTTATTCGGAAAAGGCTCACTGAAGCCTAAGTGCGCTTACAGCTTTCCCTGGGGTAGTCAGGCCGGCATCCACTGATGTGGTAGTAGTTGCCCAATCACACTGGCCCGCTGCGTCGGAAGCCGCGTCAGCACCTCTTTGAGGTAGGCATACGGATCATGCCCATTCATGCGAGCCGATCA  >  NC_002947/758195‑758482

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: