Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 746,774 C→T 5.8% G332G (GGG→GGA)  abrB ← regulator of aidB expression; inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913746,7740CT 5.8% 221.6 / 3.2 69G332G (GGG→GGA) abrBregulator of aidB expression; inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (32/33);  new base T (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

ATCGAGATGGCGAACCGTGCGAACGTTGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGTCATAAAATCAATGTGCAT  >  NC_000913/746629‑746920
                                                                                                                                                 |                                                                                                                                                  
aTCGAGATGGCGAACCGTGCGAACGTTGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCt                                                                                                                                                 <  2:143287/149‑1 (MQ=255)
    agaTGGCGAACCGTGCGAACGTTGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCa                                                                                                                                             >  4:348747/1‑149 (MQ=255)
          cGAACCGTGCGAACGTTGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCaaaa                                                                                                                                       <  2:364072/149‑1 (MQ=255)
            aaCCGTGCGAACGTTGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATa                                                                                                                                     <  2:339881/149‑1 (MQ=255)
                         ttGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTa                                                                                                                        <  3:509604/149‑1 (MQ=255)
                             ggCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGt                                                                                                                    >  4:478798/1‑149 (MQ=255)
                              gcagcaTCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGtt                                                                                                                   <  1:171829/149‑1 (MQ=255)
                               cagcaTCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGACCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGttt                                                                                                                  <  1:7853/149‑1 (MQ=255)
                                agcaTCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTg                                                                                                                 <  2:95875/149‑1 (MQ=255)
                                      atTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGc                                                                                                           >  3:718753/1‑149 (MQ=255)
                                       ttGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGcc                                                                                                          >  1:662239/1‑149 (MQ=255)
                                           gaAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATaa                                                                                                      >  1:30944/1‑149 (MQ=255)
                                            aaaaCCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAAt                                                                                                     <  3:299880/149‑1 (MQ=255)
                                               aCCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAg                                                                                                  >  2:192827/1‑149 (MQ=255)
                                                 cACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGTCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAgcg                                                                                                >  3:796883/1‑149 (MQ=255)
                                                    gCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCa                                                                                             <  2:360196/149‑1 (MQ=255)
                                                         cATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCg                                                                                        <  4:314276/149‑1 (MQ=255)
                                                          aTCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGttt                                                                                                                  >  2:435730/1‑122 (MQ=255)
                                                          aTCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGttt                                                                                                                  <  1:435730/122‑1 (MQ=255)
                                                             tcGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCAt                                                                                    <  4:718753/149‑1 (MQ=255)
                                                              cGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCAtt                                                                                   <  3:653231/149‑1 (MQ=255)
                                                               gCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTg                                                                                  >  4:388043/1‑149 (MQ=255)
                                                                cccGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGc                                                                                 >  4:196783/1‑149 (MQ=255)
                                                                 ccGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCt                                                                                >  2:544184/1‑149 (MQ=255)
                                                                   gTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCAt                                                                                    >  1:667818/1‑143 (MQ=255)
                                                                   gTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCAt                                                                                    <  2:667818/143‑1 (MQ=255)
                                                                    ttttACTGGTGCCCTGGGTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGTCATTGCTccc                                                                             <  2:34356/149‑1 (MQ=255)
                                                                     tttACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCt                                                                            <  3:669501/149‑1 (MQ=255)
                                                                         cTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTgcgg                                                                        <  1:282339/149‑1 (MQ=255)
                                                                                 cTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGTAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCGTCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGAcc                                                                <  4:692526/149‑1 (MQ=255)
                                                                                  tGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCg                                                               <  3:388043/149‑1 (MQ=255)
                                                                                      tGCCAGCACTATGCCGACCGCTTCGGCGCATAGGTTGAAGTAAACCGCGCAATGTCAAACCCTGTCAGCAAAATACTGAACAGGCGTAGGGGTAGCATCGCCATAATGAGCGCCATATAGGCATTGCTCCCTGCGGCGATGACCGCCAc                                                           >  4:726090/1‑149 (MQ=255)
                                                                                      tGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCAc                                                           <  3:196783/149‑1 (MQ=255)
                                                                                      tGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCAc                                                           <  1:748222/149‑1 (MQ=255)
                                                                                          aGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTa                                                       >  2:564617/1‑149 (MQ=255)
                                                                                             aCTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCa                                                    >  3:738876/1‑149 (MQ=255)
                                                                                              cTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGAt                                                                    >  4:555211/1‑132 (MQ=255)
                                                                                              cTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGAt                                                                    <  3:555211/132‑1 (MQ=255)
                                                                                              cTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCaa                                                   >  3:788218/1‑149 (MQ=255)
                                                                                              cTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCaa                                                   >  1:47590/1‑149 (MQ=255)
                                                                                              cTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCaa                                                   <  3:478798/149‑1 (MQ=255)
                                                                                                       ccGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGTCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCg                                                                                        >  4:645157/1‑103 (MQ=255)
                                                                                                       ccGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGTCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCg                                                                                        <  3:645157/103‑1 (MQ=255)
                                                                                                              ggCGCATAGGTTGAAATAAACCGCGCAATGGCAGGTCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCtggt                                   <  4:796883/149‑1 (MQ=255)
                                                                                                                   aTAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGa                              <  2:30944/149‑1 (MQ=255)
                                                                                                                     aGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAgg                            >  3:167505/1‑149 (MQ=255)
                                                                                                                        ttGCAAAAAACCGCGCAATGGCAGGCCCTCTCAGCAAAATACTGTACAGACATCGGGCTGGCATCGCCATAATCAGCGCCATATAGGCATTGCTGCCTGCGGCGAGGACCGCCACGGTACCAGGCCCGCCGGAGCTGGTGGCGCGGTAg                         >  2:66070/1‑149 (MQ=255)
                                                                                                                           aaaTAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCa                      <  1:592415/149‑1 (MQ=255)
                                                                                                                            aaTAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAg                     >  1:257609/1‑149 (MQ=255)
                                                                                                                            aaTAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAg                     >  1:4669/1‑149 (MQ=255)
                                                                                                                            aaTAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAg                     <  4:602820/149‑1 (MQ=255)
                                                                                                                            aaTAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAg                     >  3:602820/1‑149 (MQ=255)
                                                                                                                            aaTAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAg                     <  2:4669/149‑1 (MQ=255)
                                                                                                                             aTAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCTTAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGt                    >  4:546505/1‑149 (MQ=255)
                                                                                                                              tAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGTc                   <  2:257609/149‑1 (MQ=255)
                                                                                                                              tAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGTc                   <  1:544184/149‑1 (MQ=255)
                                                                                                                                aaCCGCGAAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGTCAt                 <  4:738876/149‑1 (MQ=255)
                                                                                                                                 aCCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGTCATa                >  4:155938/1‑149 (MQ=255)
                                                                                                                                 aCCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGTCATa                >  1:504556/1‑149 (MQ=255)
                                                                                                                                   cgcgCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGTCATaaa              >  2:225757/1‑149 (MQ=255)
                                                                                                                                   cgcgCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGTCATaaa              >  4:693659/1‑149 (MQ=255)
                                                                                                                                        aaTGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGTCATAAAATCaa         >  2:370980/1‑149 (MQ=255)
                                                                                                                                        aaTGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGTCATAAAATCaa         >  4:457/1‑149 (MQ=255)
                                                                                                                                         aTGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGTCATAAAATCaa         >  4:170280/1‑148 (MQ=255)
                                                                                                                                         aTGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGTCATAAAATCaa         <  3:170280/148‑1 (MQ=255)
                                                                                                                                             cAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGTCATAAAATCAATGTGc    <  3:348747/149‑1 (MQ=255)
                                                                                                                                               ggCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGTCATAAAATCAATGTGCAt  >  2:657999/1‑149 (MQ=255)
                                                                                                                                                gCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGTCATAAAATCAATGTGc    >  4:424575/1‑146 (MQ=255)
                                                                                                                                                gCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGTCATAAAATCAATGTGc    <  3:424575/146‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ATCGAGATGGCGAACCGTGCGAACGTTGTGGCAGCATCATTGAGAAAACCACGCTGTCATCTCGCCCGTTTTACTGGTGCCCTGGCTGCCAGCACTAGGCCGACCGCTTCGGCGCATAGGTTGAAATAAACCGCGCAATGGCAGGCCCTGTCAGCAAAATACTGAACAGGCGTAGGGTTTGCATCGCCATAATGAGCGCCATATCGGCATTGCTCCCTGCGGCGATGACCGCCACGGTATCAAGCCCGCCGGGGCTGGTGGCGAGGTAGGCAGTCATAAAATCAATGTGCAT  >  NC_000913/746629‑746920

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