Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 701,130 T→C 5.2% P155P (CCA→CCG)  nagC ← N‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913701,1300TC 5.2% 246.6 / 3.2 77P155P (CCA→CCG) nagCN‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (39/34);  new base C (2/2);  total (41/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 = 3.34e-01

GCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGCGGGTAATGTTCTTCTGCCAGCACTTTGCTGC  >  NC_000913/700985‑701278
                                                                                                                                                 |                                                                                                                                                    
gCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGc                                                                                                                                                   <  3:165482/149‑1 (MQ=255)
     gccagcgccagACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGAt                                                                                                                                              <  2:42452/149‑1 (MQ=255)
         gcgccagACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCAcc                                                                                                                                          <  2:379430/149‑1 (MQ=255)
            ccagACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGaa                                                                                                                                       >  3:291277/1‑149 (MQ=255)
                  tACGGATCTCGTGACCAACGAAACAGGTCCCTTTAAACCGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCGGAATATGCGGCATGTAAGGAATTTTGCCGCTGTCCGGGTAAACAAGCCCTGGCAGGATCACCGAAAGcgag                                                                                                                                 <  2:465422/149‑3 (MQ=255)
                    cGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCTCCAGTTTTCTACCTGAATATGCGGCATGTAATGTATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGAt                                                                                                                               <  2:301574/149‑1 (MQ=255)
                    cGGATATCATGACCCACGAAACAGGTCACTTAAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTACGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGAt                                                                                                                               <  3:350758/149‑1 (MQ=255)
                      gATATCGTGACCGACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCTGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCa                                                                                                                             >  4:355707/1‑149 (MQ=255)
                           cGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTcg                                                                                                                        >  1:48517/1‑149 (MQ=255)
                                 cAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGt                                                                                                                  >  4:393775/1‑149 (MQ=255)
                                 cAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGt                                                                                                                  >  1:287597/1‑149 (MQ=255)
                                 cAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGt                                                                                                                  >  1:74508/1‑149 (MQ=255)
                                 cAACGAAACAGGGCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGt                                                                                                                  >  1:86032/1‑149 (MQ=255)
                                   aCGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGttt                                                                                                                <  1:262435/149‑1 (MQ=255)
                                       aaCAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTgcgc                                                                                                            >  4:74016/1‑149 (MQ=255)
                                       aaCAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTgcgc                                                                                                            <  2:389875/149‑1 (MQ=255)
                                       aaCAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCCGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTgcgc                                                                                                            >  2:203828/1‑149 (MQ=255)
                                       aaCAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCCGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTgcgc                                                                                                            >  4:321791/1‑149 (MQ=255)
                                          aGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTgg                                                                                                         <  2:287597/149‑1 (MQ=255)
                                          aGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCAGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTgg                                                                                                         <  2:86032/149‑1 (MQ=255)
                                            gTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTa                                                                                                       <  3:149260/149‑1 (MQ=255)
                                                 tttAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTTGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTAt                                                                                                  <  2:443886/149‑1 (MQ=255)
                                                  ttAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATc                                                                                                 <  2:203970/149‑1 (MQ=255)
                                                      aaCGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGAt                                                                                                                               >  4:71096/1‑115 (MQ=255)
                                                      aaCGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGAt                                                                                                                               <  3:71096/115‑1 (MQ=255)
                                                      aaCGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTg                                                                                                          >  1:371851/1‑136 (MQ=255)
                                                      aaCGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTg                                                                                                          <  2:371851/136‑1 (MQ=255)
                                                      aaCGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATa                                                                                             >  3:121422/1‑149 (MQ=255)
                                                             ttcCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCCGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAg                                                                                      <  1:203828/149‑1 (MQ=255)
                                                             ttcCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCCGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAg                                                                                      <  3:321791/149‑1 (MQ=255)
                                                              tcCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGc                                                                                     <  4:291277/149‑1 (MQ=255)
                                                               ccAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAGGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCa                                                                                    >  1:448842/1‑149 (MQ=255)
                                                                 agagCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAAt                                                                                  >  4:69259/1‑149 (MQ=255)
                                                                        cTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTc                                                                           <  3:69259/149‑1 (MQ=255)
                                                                          aCCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCCATGGCATTc                                                                           <  4:248500/147‑1 (MQ=255)
                                                                          aCCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTc                                                                           >  3:248500/1‑147 (MQ=255)
                                                                            cAGCCCCCAGTTTTCTACATGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACa                                                                       >  4:186142/1‑149 (MQ=255)
                                                                                  ccAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATaaa                                                                                           <  3:275552/123‑1 (MQ=255)
                                                                                  ccAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATaaa                                                                                           >  4:275552/1‑123 (MQ=255)
                                                                                  ccAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCAt                                                                 >  1:311711/1‑149 (MQ=255)
                                                                                  ccAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCAt                                                                 <  3:459738/149‑1 (MQ=255)
                                                                                      ttttCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTc                                                             >  3:269862/1‑149 (MQ=255)
                                                                                          cTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGc                                                         <  4:148876/149‑1 (MQ=255)
                                                                                            aCCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGt                                                       >  1:182646/1‑149 (MQ=255)
                                                                                             ccTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGtt                                                      >  2:574828/1‑149 (MQ=255)
                                                                                              cTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGttt                                                     >  4:247978/1‑149 (MQ=255)
                                                                                              cTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGttt                                                     >  3:51317/1‑149 (MQ=255)
                                                                                              cTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGttt                                                     >  4:92280/1‑149 (MQ=255)
                                                                                              cTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGttt                                                     <  3:213345/149‑1 (MQ=255)
                                                                                              cTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGttt                                                     >  2:81372/1‑149 (MQ=255)
                                                                                              cTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGttt                                                     >  1:202600/1‑149 (MQ=255)
                                                                                              cTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGttt                                                     >  2:68519/1‑149 (MQ=255)
                                                                                              cTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAAAAGGGCATGTTCCAGCGttt                                                     >  4:117874/1‑149 (MQ=255)
                                                                                               tGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTg                                                    >  4:194073/1‑149 (MQ=255)
                                                                                                gAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGc                                                   >  1:492201/1‑149 (MQ=255)
                                                                                                gAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGc                                                   >  3:269360/1‑149 (MQ=255)
                                                                                                 aaTATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCt                                                  <  1:234183/149‑1 (MQ=255)
                                                                                                   tatGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTgg                                                <  3:468134/149‑1 (MQ=255)
                                                                                                     tGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGt                                              <  3:407312/149‑1 (MQ=255)
                                                                                                      gCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGAt                                                                                                                                              <  3:561241/52‑1 (MQ=255)
                                                                                                      gCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGAt                                                                                                                                              >  4:561241/1‑52 (MQ=255)
                                                                                                                aTGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCa                                   >  3:38242/1‑149 (MQ=255)
                                                                                                                  gAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGc                                 <  3:92280/149‑1 (MQ=255)
                                                                                                                  gAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGc                                 <  3:117874/149‑1 (MQ=255)
                                                                                                                     ttttGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGCggg                              <  1:345872/149‑1 (MQ=255)
                                                                                                                           cGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGCGGGTAATGt                        <  4:38242/149‑1 (MQ=255)
                                                                                                                             cTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGCGGGTAATGttc                      <  3:247978/149‑1 (MQ=255)
                                                                                                                             cTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAAGCGCGTTCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGCGGGTAATGttc                      <  1:446934/149‑1 (MQ=255)
                                                                                                                                   gggTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGCGGGTAATGTTCTTCTGc                <  3:146098/149‑1 (MQ=255)
                                                                                                                                   gggTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGCGGGTAATGTTCTTCTGc                >  2:528340/1‑149 (MQ=255)
                                                                                                                                   gggTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGCGGGTAATGTTCTTCTGc                <  1:299592/149‑1 (MQ=255)
                                                                                                                                     gTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGCGGGTAATGTTCTTCTGCCa              >  2:376028/1‑149 (MQ=255)
                                                                                                                                         acaaGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGCGGGTAATGTTCTTCTGCCAGCAc          >  1:206098/1‑149 (MQ=255)
                                                                                                                                         acaaGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGCGGGTAATGTTCTTCTGCCAGCAc          >  3:223826/1‑149 (MQ=255)
                                                                                                                                         acaaGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGCGGGTAATGTTCTTCTGCCAGCAc          >  1:210550/1‑149 (MQ=255)
                                                                                                                                           aaGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGCGGGTAATGTTCTTCTGCCAGCACtt        >  3:475751/1‑149 (MQ=255)
                                                                                                                                              cccTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGCGGGTAATGTTCTTCTGCCAGCACTTtgc     >  3:548730/1‑149 (MQ=255)
                                                                                                                                                 tGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGCGGGTAATGTTCTTCTGCCAGCACTTtgctgc  <  3:527690/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAACAAGCCCTGGCAGGATCACCGAAATCGCGATCAGCTCGCGCAGTTTGCGCTGGTAGCTATCAATAAACTGAGCAATGGCATTCAACAGGGCATGTTCCAGCGTTTGCTGGGTACGTTCCGGCAGCGGGTAATGTTCTTCTGCCAGCACTTTGCTGC  >  NC_000913/700985‑701278

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑