Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,483,443 T→A 5.5% L89Q (CTG→CAG)  puuA‑I → glutamate‑putrescine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,483,4430TA 5.5% 224.7 / 4.6 73L89Q (CTG→CAG) puuA‑Iglutamate‑putrescine ligase
Reads supporting (aligned to +/- strand):  ref base T (38/31);  new base A (2/2);  total (40/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 = 4.64e-01

CGCTACCCCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTGGCATT  >  NC_002947/2483295‑2483591
                                                                                                                                                     |                                                                                                                                                     
cgcTACCCCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTg                                                                                                                                                      <  2:210576/150‑1 (MQ=255)
   tACCCCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTTCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTg                                                                                                                                                   >  1:18569/1‑150 (MQ=255)
       ccATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCgg                                                                                                                                               >  1:605478/1‑150 (MQ=255)
        cATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCggg                                                                                                                                              >  2:1358395/1‑150 (MQ=255)
          tGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTACCCGGCACCAGGCGCCgggca                                                                                                                                            <  1:888251/150‑3 (MQ=255)
                  cTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAg                                                                                                                                    >  1:1226383/1‑150 (MQ=255)
                  cTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAg                                                                                                                                    >  2:715658/1‑150 (MQ=255)
                   tGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGc                                                                                                                                   >  1:647942/1‑150 (MQ=255)
                     gAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCg                                                                                                                                 <  1:1403038/150‑1 (MQ=255)
                       aaaaGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGt                                                                                                                               >  1:1450563/1‑150 (MQ=255)
                        aaaGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTg                                                                                                                              >  2:1120190/1‑150 (MQ=255)
                           gTCGCAGCCGGTGCGCCTCTGAAGCTGCCGGAGAACTGCTTGCTGCGGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTgcgc                                                                                                                           >  1:35710/1‑150 (MQ=255)
                              gCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTgg                                                                                                                       <  1:1120190/150‑1 (MQ=255)
                                   cGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGaaaa                                                                                                                  <  2:704598/150‑1 (MQ=255)
                                     gTGCGCCTCTGACGCGGCCGCAGAACTGCGTGCTGCTGGGTACCCAGAGCGGCGTGTACGCGGCGGTGTTGTGTTTTTTCCATGACGGCGGCCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAgg                                                                                                                <  1:115982/150‑1 (MQ=255)
                                       gcgcCTCTGAAGCTGCCGCAGAACTGCGTGTTGTTCGGTCCCAAGGGGCGGGTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTTGGAAAAGGaa                                                                                                              >  1:660290/1‑150 (MQ=255)
                                       gcgcCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGaa                                                                                                              <  1:715658/150‑1 (MQ=255)
                                           ctctGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGa                                                                                                          >  2:949596/1‑150 (MQ=255)
                                             ctGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATc                                                                                                        <  1:949596/150‑1 (MQ=255)
                                               gAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCg                                                                                                                                 >  1:964556/1‑124 (MQ=255)
                                               gAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCg                                                                                                                                 <  2:964556/124‑1 (MQ=255)
                                                  gCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCa                                                                                                   <  2:647942/150‑1 (MQ=255)
                                                       cGCAGAACTGCGTGATGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGCGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGc                                                                                              <  2:18569/150‑1 (MQ=255)
                                                            aaCTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATc                                                                                         <  2:1226383/150‑1 (MQ=255)
                                                                         tgGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCg                                                                                                                                 >  2:759733/1‑98 (MQ=255)
                                                                         tgGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCg                                                                                                                                 <  1:759733/98‑1 (MQ=255)
                                                                            gTACCCAGGGCGGCCTTTAACCGATTGGTGCCTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCAGGTACGCTGAAGCCGGTGCAGTCAGGAGTTGAAGTTATCGCCCAGATGCTGATCACGTCCGACGGCACTg                                                                         <  1:1518744/150‑1 (MQ=255)
                                                                              aCCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGCCGCCGCTATGCGTGGGCAATCGTTGAAACGTCCAGAGGCTGATCACGTCCGACGGCACTGCa                                                                       >  2:357401/1‑150 (MQ=255)
                                                                                ccAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCaaa                                                                     >  2:1091065/1‑150 (MQ=255)
                                                                                ccAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCaaa                                                                     <  1:1153032/150‑1 (MQ=255)
                                                                                 cAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAg                                                                    <  2:605478/150‑1 (MQ=255)
                                                                                       ggcCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGCCTCGGAAGCTGCCCGGCGCCTCATGCTGGGTACGTTGAATCCGGTGCGCTGGGAAAAGGAAGTGATGGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATc                                                              <  2:1004676/150‑1 (MQ=255)
                                                                                       ggcCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGaa                                                                                                              <  1:1584226/102‑1 (MQ=255)
                                                                                       ggcCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGaa                                                                                                              >  2:1584226/1‑102 (MQ=255)
                                                                                       ggcCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCAGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTggg                                                                                                                      >  2:265076/1‑94 (MQ=255)
                                                                                       ggcCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCAGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTggg                                                                                                                      <  1:265076/94‑1 (MQ=255)
                                                                                         ccTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACgg                                                                              >  1:814185/1‑132 (MQ=255)
                                                                                         ccTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACgg                                                                              <  2:814185/132‑1 (MQ=255)
                                                                                           tttACTCGAGTGGTGACTATTGCTTCGCCGACGGCGACCCGAATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCCGGGAAACGGAAGTGATCGCCCCGATGCTGATCACGTCCGCCGGCACTGCGAAGCCAATCGAAt                                                          >  2:1408577/1‑150 (MQ=255)
                                                                                             tACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCgg                                                                                                                                               >  2:254349/1‑64 (MQ=255)
                                                                                             tACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCgg                                                                                                                                               <  1:254349/64‑1 (MQ=255)
                                                                                               cccGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGa                                                                <  2:1543082/140‑1 (MQ=255)
                                                                                               cccGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGa                                                                >  1:1543082/1‑140 (MQ=255)
                                                                                                ccGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGaa                                                                                                              >  2:985178/1‑93 (MQ=255)
                                                                                                ccGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGaa                                                                                                              <  1:985178/93‑1 (MQ=255)
                                                                                                 cGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGc                                                    >  1:786528/1‑150 (MQ=255)
                                                                                                       gTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCCAATTCGACCCTAGTg                                              >  2:1271462/1‑150 (MQ=255)
                                                                                                       gTGACTATTGCTTCGCCGACGGCGACCCGGATGCTAACCGGCTCCTGGGGCCGGGTACGCTTAAGCCGGAACGCTGGGAAAAGGAAGTGATCGCCCATATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGCATGCGAGGCTCGTg                                              >  1:1349277/1‑150 (MQ=255)
                                                                                                        tGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATc                                                                                         <  2:620695/106‑1 (MQ=255)
                                                                                                        tGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATc                                                                                         >  1:620695/1‑106 (MQ=255)
                                                                                                                 gCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTgg                                     >  2:1370376/1‑150 (MQ=255)
                                                                                                                   ttCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTAGTGCTGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGgcg                                   >  1:1549175/1‑150 (MQ=255)
                                                                                                                        cGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCgtgt                              >  2:938681/1‑150 (MQ=255)
                                                                                                                              cGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGGACGCTCAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGCGCCTCGTGAAGTGCTGGCGCGTGgtcttcc                        >  2:234171/1‑149 (MQ=255)
                                                                                                                                  ccGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGt                    >  1:978714/1‑150 (MQ=255)
                                                                                                                                  ccGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGt                    >  2:1426876/1‑150 (MQ=255)
                                                                                                                                     gATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTg                 <  1:1358395/150‑1 (MQ=255)
                                                                                                                                          tgcCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCaa            >  1:571883/1‑150 (MQ=255)
                                                                                                                                          tacccGGCACCAGGCGCCGGGCACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGtt                             >  2:1544636/3‑133 (MQ=255)
                                                                                                                                          tacccGGCACCAGGCGCCGGGCACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGtt                             <  1:1544636/131‑1 (MQ=255)
                                                                                                                                            cccGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTg                                                                                            >  1:978734/1‑68 (MQ=255)
                                                                                                                                            cccGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTg                                                                                            <  2:978734/68‑1 (MQ=255)
                                                                                                                                            cccGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGcc              >  1:340048/1‑146 (MQ=255)
                                                                                                                                            cccGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGcc              <  2:340048/146‑1 (MQ=255)
                                                                                                                                            cccGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGc          <  1:357401/150‑1 (MQ=255)
                                                                                                                                            cccGGCACCAGGCGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGc          >  1:488794/1‑150 (MQ=255)
                                                                                                                                               ggCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTg       <  1:1370376/150‑1 (MQ=255)
                                                                                                                                                gcgcCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTgg      <  2:1450563/150‑1 (MQ=255)
                                                                                                                                                  gcCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAGGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGc                                                    >  1:1491441/1‑102 (MQ=255)
                                                                                                                                                  gcCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAGGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGc                                                    <  2:1491441/102‑1 (MQ=255)
                                                                                                                                                    cTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCgtg                               >  2:881115/1‑121 (MQ=255)
                                                                                                                                                    cTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCgtg                               <  1:881115/121‑1 (MQ=255)
                                                                                                                                                    cTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTGGCAtt  >  1:1428395/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCTACCCCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTGGCATT  >  NC_002947/2483295‑2483591

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