Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 58,747 T→A 10.6% intergenic (‑658/+30) PP_5424 ← / ← PP_0050 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294758,7470TA10.6% 179.6 / 9.9 66intergenic (‑658/+30)PP_5424/PP_0050hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/26);  new base A (5/2);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GGTGAGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGCGCTTTC  >  NC_002947/58600‑58894
                                                                                                                                                     |                                                                                                                                                   
ggTGAGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCgct                                                                                                                                                     <  2:32256/148‑1 (MQ=255)
 gTGAGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCgctgc                                                                                                                                                   >  3:102084/1‑149 (MQ=255)
    aGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTGACGGGGTCTGGCCTATGACCTATTTTGTTAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCtgctgc                                                                                                                                                >  1:163143/1‑149 (MQ=255)
    aGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCtgctgc                                                                                                                                                <  2:13138/149‑1 (MQ=255)
     gTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGGCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCtgctgc                                                                                                                                                <  3:224949/148‑1 (MQ=255)
      tAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCgg                                                                                                                                              >  1:36553/1‑149 (MQ=255)
           aTTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCttt                                                                                                                                         >  6:37216/1‑149 (MQ=255)
           aTTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACCGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGGCACTAAAAAGCCCGCAGAAGCGGGCttt                                                                                                                                         >  3:2575/1‑149 (MQ=255)
             ttCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGGCACTAAAAAGCCCGCAGCAGCGGGCtttt                                                                                                                                        <  4:2575/148‑1 (MQ=255)
               ctcGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                      <  5:156253/148‑1 (MQ=255)
               ctcGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                      <  3:265657/148‑1 (MQ=255)
               ctcGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                      <  7:115627/148‑1 (MQ=255)
                 cGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTg                                                                                                                                   >  4:269724/1‑149 (MQ=255)
                 cGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTg                                                                                                                                   <  5:37216/149‑1 (MQ=255)
                  ggCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGc                                                                                                                                  <  8:21918/149‑1 (MQ=255)
                   gCTTTTCTGTGAAAAGCACCCCTTCATTTCGACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGGTGCGGGCTTTTTAGTg                                                                                                                                   >  1:28683/1‑148 (MQ=255)
                   gCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGc                                                                                                                                  >  8:224784/1‑148 (MQ=255)
                   gCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGc                                                                                                                                  <  7:224784/148‑1 (MQ=255)
                     ttttCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTg                                                                                                                               >  1:13278/1‑149 (MQ=255)
                            tgAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCaaa                                                                                                                           <  4:102084/146‑1 (MQ=255)
                             gAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACaa                                                                                                                        >  2:93009/1‑148 (MQ=255)
                               aaaGCACACCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCa                                                                                                                     <  8:117400/149‑1 (MQ=255)
                                      cccTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCg                                                                                                              >  3:3611/1‑149 (MQ=255)
                                        cTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGt                                                                                                            >  2:85577/1‑149 (MQ=255)
                                                cacaAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATg                                                                                                     <  4:3611/148‑1 (MQ=255)
                                                  acaAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCGTTTTAGTGCCTGCAAACAATCa                                                                                                                     >  8:271746/1‑131 (MQ=255)
                                                                gAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGc                                                                                     <  6:177988/149‑1 (MQ=255)
                                                                     ttCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACGAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAAACAGTAGCCGGTGAGAATGTTCACGATCACACCGATgg                                                                                  >  1:229195/1‑148 (MQ=255)
                                                                        cAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCAAGATCACACCGCTGGCAATg                                                                             >  6:124559/1‑149 (MQ=255)
                                                                         aaTTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATgg                                                                            <  1:85577/149‑1 (MQ=255)
                                                                         aaTTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATg                                                                             <  6:102047/148‑1 (MQ=255)
                                                                         aaTTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATg                                                                             >  5:102047/1‑148 (MQ=255)
                                                                                 ggggTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCAGCAGCGGGCTTTTTa                                                                                                                                      >  2:254991/1‑83 (MQ=255)
                                                                                 ggggTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCAGCAGCGGGCTTTTTa                                                                                                                                      <  1:254991/83‑1 (MQ=255)
                                                                                    gTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACa                                                                 >  1:230420/1‑149 (MQ=255)
                                                                                      cTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATa                                                               <  2:229195/149‑1 (MQ=255)
                                                                                           cTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGc                                                          <  1:93009/149‑1 (MQ=255)
                                                                                           cTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGc                                                          >  6:240716/1‑149 (MQ=255)
                                                                                              tCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCa                                                        <  7:164396/148‑1 (MQ=255)
                                                                                                aCCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAGAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATgg                                                                            <  4:52139/126‑1 (MQ=255)
                                                                                                aCCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAGAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATgg                                                                            >  3:52139/1‑126 (MQ=255)
                                                                                                aCCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCAtt                                                      >  6:287059/1‑148 (MQ=255)
                                                                                                  cTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTca                                                    >  6:167636/1‑148 (MQ=255)
                                                                                                  cTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTca                                                    >  2:285860/1‑148 (MQ=255)
                                                                                                   tATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTcaca                                                  >  6:225160/1‑149 (MQ=255)
                                                                                                    aTTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGc                                                                                                                              <  3:85082/72‑1 (MQ=255)
                                                                                                    aTTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGc                                                                                                                              >  4:85082/1‑72 (MQ=255)
                                                                                                        tgtgAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACa                                              <  2:230420/148‑1 (MQ=255)
                                                                                                        tgtgAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACa                                              <  5:167636/148‑1 (MQ=255)
                                                                                                        tgtgAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACa                                              <  2:163143/148‑1 (MQ=255)
                                                                                                         gtgAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGc                                            >  5:155281/1‑149 (MQ=255)
                                                                                                          tgAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGc                                            >  5:2537/1‑148 (MQ=255)
                                                                                                             aCCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCAtt                                                      >  2:28763/1‑135 (MQ=255)
                                                                                                             aCCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCAtt                                                      <  1:28763/135‑1 (MQ=255)
                                                                                                               cTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTgg                                                                                  <  7:126543/105‑1 (MQ=255)
                                                                                                               cTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTgg                                                                                  >  8:126543/1‑105 (MQ=255)
                                                                                                                tACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGc                                            >  1:125914/1‑142 (MQ=255)
                                                                                                                  cAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACgatg                                    <  6:2537/148‑1 (MQ=255)
                                                                                                                    aaGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACgatgat                                  >  7:159802/1‑149 (MQ=255)
                                                                                                                        tgtAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCAc                               <  4:200820/148‑1 (MQ=255)
                                                                                                                          tAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGa                             <  4:66645/148‑1 (MQ=255)
                                                                                                                                gCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTcac                                                   <  7:53665/120‑1 (MQ=255)
                                                                                                                                gCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTcac                                                   >  8:53665/1‑120 (MQ=255)
                                                                                                                                gCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGgcg                      >  5:272840/1‑149 (MQ=255)
                                                                                                                                         tAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGt             >  7:149130/1‑149 (MQ=255)
                                                                                                                                         tAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGACAACGATGATCACGAGGTGGCGATTACAGGt             >  2:94279/1‑149 (MQ=255)
                                                                                                                                         tAAAAAGCCCGCAGCAGCGGGCTTTTTAGTGTCTGCAACCAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGt             >  7:278269/1‑149 (MQ=255)
                                                                                                                                           aaaaGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTgg           >  7:232503/1‑149 (MQ=255)
                                                                                                                                           aaaaGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATCTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTgg           >  3:172002/1‑149 (MQ=255)
                                                                                                                                           aaaaGCCCGCAGCAGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTgg           >  8:78432/1‑149 (MQ=255)
                                                                                                                                             aaGCCCGCAGCAGCGGGCTTTTTAGTGTCTGCAAACAATCAGTAGCAGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGgcg         >  7:35721/1‑149 (MQ=255)
                                                                                                                                                    ctgctgCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGCGCTTTc  <  8:159802/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGTGAGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGCGCTTTC  >  NC_002947/58600‑58894

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: