Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,444,632 C→A 6.3% Q86H (CAG→CAT sulA ← cell division inhibitor SulA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,444,6320CA 6.3% 201.1 / 3.1 63Q86H (CAG→CATsulAcell division inhibitor SulA
Reads supporting (aligned to +/- strand):  ref base C (25/34);  new base A (2/2);  total (27/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 = 9.62e-01

CATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCTGAACAGCTCGGGCTGGCT  >  NC_002947/2444487‑2444778
                                                                                                                                                 |                                                                                                                                                    
cATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTgg                                                                                                                                                   <  7:143983/149‑1 (MQ=255)
    gACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTgacc                                                                                                                                                 >  6:148482/1‑144 (MQ=255)
     aCGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAACGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCa                                                                                                                                               <  3:146461/148‑1 (MQ=255)
     aCGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAg                                                                                                                                              >  1:376183/1‑149 (MQ=255)
     aCGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAg                                                                                                                                              >  7:352777/1‑149 (MQ=255)
          gCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCtt                                                                                                                                          <  8:194741/148‑1 (MQ=255)
           cccAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCTGGCATCGCGCAGCCAAGCCTGACCCAGGCCTgg                                                                                                                                        <  1:52326/149‑1 (MQ=255)
           cccAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTg                                                                                                                                         <  1:201167/148‑1 (MQ=255)
            ccagccagCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGc                                                                                                                                       <  7:53566/149‑1 (MQ=255)
            ccagccagCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGc                                                                                                                                       >  7:36137/1‑149 (MQ=255)
               gccagcTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGACCCAGGCCTGGctgc                                                                                                                                    <  5:148482/149‑4 (MQ=255)
                 cagcTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTgg                                                                                                                                  <  8:201635/149‑1 (MQ=255)
                   gcTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTgggg                                                                                                                                <  1:348008/149‑1 (MQ=255)
                       accacGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGa                                                                                                                             <  6:99765/148‑1 (MQ=255)
                       accacGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCCCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGAt                                                                                                                            <  1:309246/149‑1 (MQ=255)
                         cacGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGAt                                                                                                                            >  6:39809/1‑147 (MQ=255)
                         cacGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGAt                                                                                                                            >  7:275186/1‑147 (MQ=255)
                          acGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAg                                                                                                                         <  1:202253/149‑1 (MQ=255)
                           cGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAgg                                                                                                                        <  5:39809/149‑1 (MQ=255)
                                tgGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTcag                                                                                                                   <  5:211129/149‑1 (MQ=255)
                                tgGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTcag                                                                                                                   <  5:206032/149‑1 (MQ=255)
                                      aGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTcagccag                                                                                                               <  4:231392/147‑1 (MQ=255)
                                        gcgcAGGTGCATAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATGGCTGCGGGGTTGGAGCAGCAGAATACGTTCACGGGTGAGCCCGGCATCTCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCTGGGTCAGCCAGCggg                                                                                                           <  8:352777/146‑1 (MQ=255)
                                        gCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCggg                                                                                                           <  6:186177/149‑1 (MQ=255)
                                          ccAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGcc                                                                                                         >  8:30646/1‑149 (MQ=255)
                                            aGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTc                                                                                                       >  5:137204/1‑149 (MQ=255)
                                             ggTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCg                                                                                                      <  2:376183/149‑1 (MQ=255)
                                               tGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGCTCAGGGTCAGCCAGCGGGCCTCGTc                                                                                                    >  4:190525/1‑149 (MQ=255)
                                                 cAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCct                                                                                                  >  4:336334/1‑149 (MQ=255)
                                                    aGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCtt                                                                                               <  7:12198/149‑1 (MQ=255)
                                                    aGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGGTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCtt                                                                                               <  1:163081/149‑1 (MQ=255)
                                                       cTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGGTGCAGCAGCAGAATACGTTCACGGTGGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGc                                                                                            <  2:119562/149‑1 (MQ=255)
                                                        tCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCt                                                                                           >  3:27955/1‑149 (MQ=255)
                                                          acaGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCt                                                                                           <  7:350594/147‑1 (MQ=255)
                                                             ggCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGtt                                                                                      <  8:29781/149‑1 (MQ=255)
                                                                   ttGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCgca                                                                                <  8:57541/149‑1 (MQ=255)
                                                                     gCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCgcag                                                                               >  8:146784/1‑148 (MQ=255)
                                                                      cAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCgcagca                                                                             >  3:101376/1‑149 (MQ=255)
                                                                         cgcgcTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACgg                                                                          >  8:180916/1‑149 (MQ=255)
                                                                         cgcgcTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACgg                                                                          >  5:132706/1‑149 (MQ=255)
                                                                           cgcTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACgggg                                                                        <  6:137204/149‑1 (MQ=255)
                                                                            gcTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACggggg                                                                       <  6:132706/149‑1 (MQ=255)
                                                                                  aTTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCcag                                                                   >  1:104514/1‑147 (MQ=255)
                                                                                   ttGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCcagcag                                                                >  3:162707/1‑149 (MQ=255)
                                                                                          ggggTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTgg                                                         <  6:167824/149‑1 (MQ=255)
                                                                                                gcagcagcaGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCgg                                                                                                            <  4:243150/92‑1 (MQ=255)
                                                                                                gcagcagcaGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCgg                                                                                                            >  3:243150/1‑92 (MQ=255)
                                                                                                gcagcagcaGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCttt                                                           >  2:250634/1‑141 (MQ=255)
                                                                                                gcagcagcaGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCttt                                                           <  1:250634/141‑1 (MQ=255)
                                                                                                    cagcagAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAg                                             >  6:182502/1‑149 (MQ=255)
                                                                                                     agcagAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACgggg                                                                        <  6:114439/123‑1 (MQ=255)
                                                                                                     agcagAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCt                                                             >  1:226009/1‑134 (MQ=255)
                                                                                                     agcagAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCt                                                             <  2:226009/134‑1 (MQ=255)
                                                                                                       cagAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAgggg                                          <  4:27955/149‑1 (MQ=255)
                                                                                                          aaTACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGgcgc                                       >  3:255286/1‑149 (MQ=255)
                                                                                                                   aCGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGc                               >  7:226112/1‑148 (MQ=255)
                                                                                                                     ggTTGAGCCCGGCATCGCGCAGCCAGGCATGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGa                             >  2:242335/1‑148 (MQ=255)
                                                                                                                     ggTTGAGCCCGGCATCGCGCAGCCAGGCATGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGa                             >  2:247660/1‑148 (MQ=255)
                                                                                                                           gCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCg                      <  4:255286/149‑1 (MQ=255)
                                                                                                                                      cgcAGCCAGGCATGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAAGGGGCGCCGCGCAGCGACAGTTCGCTGAACAGCt            <  1:247660/149‑1 (MQ=255)
                                                                                                                                      cgcAGCCAGGCATGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCACGGGGCGCCGCGCAGCGACAGTTCGCTGAACAGCt            <  1:242335/149‑1 (MQ=255)
                                                                                                                                         cgCCAAGCCTGACCCAGGCCTGGCGGTGGGGCGATCAGGGTCAGCCAGCGCGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCTGAACAGCTCggg        >  4:618/2‑149 (MQ=255)
                                                                                                                                               gCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCTGAACAGCTCGggctggct  >  1:165096/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCTGAACAGCTCGGGCTGGCT  >  NC_002947/2444487‑2444778

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: