Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,237,611 C→G 7.1% A14G (GCC→GGC)  dadA → D‑amino acid dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,237,6110CG 7.1% 228.1 / 7.9 70A14G (GCC→GGC) dadAD‑amino acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (41/24);  new base G (2/3);  total (43/27)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

CACCATAAAAAAGCCGCATGTTGAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAGGTGTGCCTTTAAAAG  >  NC_000913/1237470‑1237760
                                                                                                                                             |                                                                                                                                                     
cACCATAAAAAAGCCGCATGTTGAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCt                                                                                                                                               >  2:1198191/1‑150 (MQ=255)
cACCATAAAAAAGCCGCATGTTGAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGCTGGCAACGCCt                                                                                                                                               <  1:88425/150‑1 (MQ=255)
 aCCATAAAAAAGCCGCATGTTGAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTg                                                                                                                                              <  1:373254/150‑1 (MQ=255)
 aCCATAAAAAAGCCGCATGTTGAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTg                                                                                                                                              <  1:477475/150‑1 (MQ=255)
  ccATAAAAAAGCCGCATGTTGAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGCTGGCAACGCCtac                                                                                                                                             <  1:799853/150‑3 (MQ=255)
   cATAAAAAAGCCGCATGTTGAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGt                                                                                                                                            <  2:1431262/150‑1 (MQ=255)
    aTAAAAAAGCCGCATGTTGAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTa                                                                                                                                           <  2:683245/150‑1 (MQ=255)
      aaaaaaGCCGCATGTTGAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGCTGGCAACGCCtaccaca                                                                                                                                         <  2:119110/150‑7 (MQ=255)
         aaaGCCGCATGTTGAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTaa                                                                                                                                      >  2:963906/1‑150 (MQ=255)
              cGCATGTTGAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATcag                                                                                                                                 <  1:680226/150‑1 (MQ=255)
                  tGTTGAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATcaggcag                                                                                                                             >  1:1911228/1‑150 (MQ=255)
                  tGTTGAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATCCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATcaggcag                                                                                                                             <  2:1030313/150‑1 (MQ=255)
                   gTTGAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATcaggcagg                                                                                                                            <  2:453567/150‑1 (MQ=255)
                   gTTGAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATcaggcagg                                                                                                                            <  2:794305/150‑1 (MQ=255)
                    ttGAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGa                                                                                                                           >  1:1506952/1‑150 (MQ=255)
                      gAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACa                                                                                                                         <  1:227861/150‑1 (MQ=255)
                        ataataTTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATg                                                                                                                       >  1:837665/1‑150 (MQ=255)
                         taataTTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGa                                                                                                                      >  2:200139/1‑150 (MQ=255)
                         taataTTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGa                                                                                                                      <  1:1198191/150‑1 (MQ=255)
                           atatTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAgg                                                                                                                    <  1:490913/150‑1 (MQ=255)
                             atTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTc                                                                                                                  >  2:1939486/1‑150 (MQ=255)
                             atTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTc                                                                                                                  >  1:1991818/1‑150 (MQ=255)
                                ttCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCAcc                                                                                                               <  2:454185/150‑1 (MQ=255)
                                         gTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGAt                                                                                                      >  2:35636/1‑150 (MQ=255)
                                          ttATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATc                                                                                                     >  1:431227/1‑150 (MQ=255)
                                            aTCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCgg                                                                                                   >  1:1966231/1‑150 (MQ=255)
                                              cAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGa                                                                                                 >  1:1315355/1‑150 (MQ=255)
                                                     tgATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGTGAATTGGTTTGGTAGGCGTTGCCAGAGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTTATCGTGAGCCgtgg                                                                                          <  2:843795/150‑1 (MQ=255)
                                                     tgATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCgggg                                                                                          >  1:1408259/1‑150 (MQ=255)
                                                       attagattaTTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGgc                                                                                        >  1:926670/1‑150 (MQ=255)
                                                        ttagattaTTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGgcc                                                                                       <  1:200139/150‑1 (MQ=255)
                                                        ttagattaTTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGgcc                                                                                       <  1:963906/150‑1 (MQ=255)
                                                           gattaTTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGgccgcc                                                                                    >  2:1235103/1‑150 (MQ=255)
                                                           gattaTTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGgccgcc                                                                                    >  2:1232588/1‑150 (MQ=255)
                                                             ttattattCTTTTACTTTATCTACCGTAATCGGACTGGCTATGCGATTTGTCATACTCCCAACAGCTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAAACAGGCAGTACATGAGGTCACCCTCATTGATCTGGAGCCGGGGCCCGCCCt                                                                                  >  2:373167/1‑150 (MQ=255)
                                                                         tACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGc                                                                      >  2:1632962/1‑150 (MQ=255)
                                                                         tACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGc                                                                      >  2:1414604/1‑150 (MQ=255)
                                                                          aCTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGcc                                                                     >  2:247413/1‑150 (MQ=255)
                                                                           cTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGCGGCCGCCCTGGAAACCAGTGCcg                                                                    >  2:351342/1‑150 (MQ=255)
                                                                                   aCCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTACCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGcc                                                            >  1:176162/1‑150 (MQ=255)
                                                                                    cgatcctCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCg                                                           <  1:1998497/144‑1 (MQ=255)
                                                                                         aTCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCaa                                                      >  2:1456137/1‑150 (MQ=255)
                                                                                          tCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCaaa                                                     >  1:1436030/1‑150 (MQ=255)
                                                                                                  gCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTccccc                                             >  1:946390/1‑150 (MQ=255)
                                                                                                   cTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCg                                            >  1:1832074/1‑150 (MQ=255)
                                                                                                      tGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGat                                         <  2:431227/150‑1 (MQ=255)
                                                                                                      tGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGat                                         >  1:964842/1‑150 (MQ=255)
                                                                                                      tGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGat                                         >  2:1605015/1‑150 (MQ=255)
                                                                                                        cGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGatat                                       >  1:1938170/1‑150 (MQ=255)
                                                                                                        cGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGatat                                       <  2:1966231/150‑1 (MQ=255)
                                                                                                                    cTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTgg                           >  1:256450/1‑150 (MQ=255)
                                                                                                                      gggAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGggc                         <  1:1414604/150‑1 (MQ=255)
                                                                                                                          aGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCCCCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCa                     >  2:1807856/1‑150 (MQ=255)
                                                                                                                           gtggtgtggtAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCAc                    >  2:1527581/1‑150 (MQ=255)
                                                                                                                           gtggtgtggtAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCAc                    <  2:1199773/150‑1 (MQ=255)
                                                                                                                              gtgtggtAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAg                 <  2:926670/150‑1 (MQ=255)
                                                                                                                              gtgtggtAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAg                 >  1:66524/1‑150 (MQ=255)
                                                                                                                               tgtggtAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAgg                >  2:955028/1‑150 (MQ=255)
                                                                                                                               tgtggtAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGACACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAgg                >  1:985904/1‑150 (MQ=255)
                                                                                                                               agtaccaGGCGCTGGCAACGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAgg                >  1:1364424/7‑150 (MQ=255)
                                                                                                                               agtaccaGGCGCTGGCAACGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTCCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAgg                >  1:1040302/7‑150 (MQ=255)
                                                                                                                                 tggtAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAGgtg              >  1:1640365/1‑150 (MQ=255)
                                                                                                                                  ggtAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAGgtgt             >  2:1507259/1‑150 (MQ=255)
                                                                                                                                   gtAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAGGtgtg            <  1:1235103/150‑1 (MQ=255)
                                                                                                                                    tAGGCGTTGCCAGCGCCTGGTACTTTTATCTGGCTGGACATGTGGTGACCGTCTTTGATCGGGTGCCGGGGGCCGCCCTGGTAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAGGTGTGc           <  2:828120/150‑1 (MQ=255)
                                                                                                                                    tAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAGGTGTGc           >  2:1001996/1‑150 (MQ=255)
                                                                                                                                     aGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAGGTGTGcc          >  2:1787596/1‑150 (MQ=255)
                                                                                                                                         gTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAGGTGTGCCTTTa      >  2:1520046/1‑150 (MQ=255)
                                                                                                                                            gCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAGGTGTGCCTTTaaaa   >  1:1595755/1‑150 (MQ=255)
                                                                                                                                            gCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAGGTGTGCCTTTaaaa   <  1:35636/150‑1 (MQ=255)
                                                                                                                                             ccAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAGGTGTGCCTTTAAAAg  >  2:627837/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
CACCATAAAAAAGCCGCATGTTGAATAATATTTTCAACTGAGTTATCAAGATGTGATTAGATTATTATTCTTTTACTGTATCTACCGTTATCGGAGTGGCTATGCGAGTTGTCATACTGGGAAGTGGTGTGGTAGGCGTTGCCAGCGCCTGGTACTTAAATCAGGCAGGACATGAGGTCACCGTCATTGATCGGGAGCCGGGGGCCGCCCTGGAAACCAGTGCCGCGAATGCCGGGCAAATCTCCCCCGGATATGCTGCACCGTGGGCGGCACCAGGTGTGCCTTTAAAAG  >  NC_000913/1237470‑1237760

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