Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,320,872 C→G 5.4% A277G (GCC→GGC)  atoS → sensory histidine kinase AtoS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,320,8720CG 5.4% 252.6 / 4.2 74A277G (GCC→GGC) atoSsensory histidine kinase AtoS
Reads supporting (aligned to +/- strand):  ref base C (33/37);  new base G (2/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

ACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATG  >  NC_000913/2320729‑2321021
                                                                                                                                                |                                                                                                                                                     
acaAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGa                                                                                                                                                 <  2:1152558/150‑1 (MQ=255)
 caAAATATTCCCCCTCGATTACCACAATTGCCCGGGGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGAc                                                                                                                                                <  1:1022690/150‑1 (MQ=255)
   aaaTATTCCCACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCg                                                                                                                                              <  1:926268/150‑1 (MQ=255)
      tatTCCCACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGGTCAATGGCAATGAcgcca                                                                                                                                            <  1:114176/150‑5 (MQ=255)
       atTCCCACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAg                                                                                                                                          <  1:1374079/150‑1 (MQ=255)
        ttCCCACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAAAGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAgg                                                                                                                                         >  1:1176946/1‑150 (MQ=255)
         tCCCACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGGTCAATGGCAATGAcgccatcg                                                                                                                                         <  2:23120/150‑8 (MQ=255)
         tCCCACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAggg                                                                                                                                        >  2:1466779/1‑150 (MQ=255)
         tCCCACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAggg                                                                                                                                        >  2:401779/1‑150 (MQ=255)
           ccACTCGATTACCACAATTTCCCTGTTAAATTTGTAAAATCATTCATAGTTTTAATATCCTCGCCCAGGGACTTCGTGAAACGCGTACACTTAACGATCTTATTATTTAAAACGCTTAAGATGGCGTCATTGCCATTGACCGCCAGGGTg                                                                                                                                      <  2:173580/150‑1 (MQ=255)
           ccACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGTCGTCATTGCCATTGACCGCCAGGGTg                                                                                                                                      <  2:1123479/150‑1 (MQ=255)
            cACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGa                                                                                                                                     >  2:1267843/1‑150 (MQ=255)
               tCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGt                                                                                                                                  <  1:913427/150‑1 (MQ=255)
               tCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGt                                                                                                                                  <  1:719458/150‑1 (MQ=255)
                      ccACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAaccacc                                                                                                                           >  2:1142592/1‑150 (MQ=255)
                      ccACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAaccacc                                                                                                                           >  1:258869/1‑150 (MQ=255)
                       cacaATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAAccacca                                                                                                                          >  2:1090984/1‑150 (MQ=255)
                        acaATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCAt                                                                                                                         <  1:851866/150‑1 (MQ=255)
                        acaATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCAt                                                                                                                         >  2:578818/1‑150 (MQ=255)
                        acaATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCAt                                                                                                                         <  1:781615/150‑1 (MQ=255)
                         caATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATg                                                                                                                        <  2:1482400/150‑1 (MQ=255)
                            ttGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAAc                                                                                                                     <  2:889355/150‑1 (MQ=255)
                              gCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAaccacc                                                                                                                           <  2:545502/142‑1 (MQ=255)
                              gCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAaccacc                                                                                                                           >  1:545502/1‑142 (MQ=255)
                              gCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAAccc                                                                                                                   <  1:326784/150‑1 (MQ=255)
                                          ggggCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGt                                                                                                       >  1:505128/1‑150 (MQ=255)
                                           gggCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGtt                                                                                                      <  1:1359802/150‑1 (MQ=255)
                                            ggCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTa                                                                                                     <  1:369904/150‑1 (MQ=255)
                                            ggCAAATCAGTCAGAGTGTTAATAACCTCGCACAGGCACTGCGTGAAACGAGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGAGTCAATGCCATTGCCCGCCAGGGTGATGTAACAACCATGAACCCAGCATCAGAACTTa                                                                                                     >  1:509811/1‑150 (MQ=255)
                                                aaTCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCAc                                                                                                 <  2:517658/150‑1 (MQ=255)
                                                            tgtTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACg                                                                                     <  1:275529/150‑1 (MQ=255)
                                                             gataataaCCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGc                                                                                    >  1:1121070/3‑150 (MQ=255)
                                                               taataaCCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCa                                                                                  >  2:689373/1‑150 (MQ=255)
                                                                aataaCCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACg                                                                                     <  2:1058529/146‑1 (MQ=255)
                                                                aataaCCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACg                                                                                     >  1:1058529/1‑146 (MQ=255)
                                                                aataaCCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCAt                                                                                 <  1:794030/150‑1 (MQ=255)
                                                                aataaCCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCAt                                                                                 >  2:803651/1‑150 (MQ=255)
                                                                 ataaCCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATg                                                                                >  2:1204542/1‑150 (MQ=255)
                                                                  taaCCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGa                                                                               >  2:1019482/1‑150 (MQ=255)
                                                                    aCCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAAc                                                                             <  2:548592/150‑1 (MQ=255)
                                                                    aCCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAAc                                                                             >  1:1074902/1‑150 (MQ=255)
                                                                         gCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATTGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTa                                                                        <  1:1028206/150‑1 (MQ=255)
                                                                                 aCTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGcc                                                                >  1:1046139/1‑150 (MQ=255)
                                                                                 aCTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGcc                                                                >  1:464699/1‑150 (MQ=255)
                                                                                 aCTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGcc                                                                <  1:306898/150‑1 (MQ=255)
                                                                                  cTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCt                                                               >  1:507280/1‑150 (MQ=255)
                                                                                         aaaCGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCa                                                        >  1:1502077/1‑150 (MQ=255)
                                                                                           aCGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCAtg                                                      <  1:1466779/150‑1 (MQ=255)
                                                                                           aCGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCAtg                                                      <  2:1095528/150‑1 (MQ=255)
                                                                                           aCGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCAtg                                                      <  2:1100909/150‑1 (MQ=255)
                                                                                                  cacTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCtt                                                              >  1:664782/1‑135 (MQ=255)
                                                                                                  cacTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCtt                                                              <  2:664782/135‑1 (MQ=255)
                                                                                                  cacTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCg                                               <  1:1173701/150‑1 (MQ=255)
                                                                                                    cTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGAc                                             >  1:618335/1‑150 (MQ=255)
                                                                                                    cTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGAc                                             >  1:1446323/1‑150 (MQ=255)
                                                                                                             cTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCTTGAAGTGTTAGGGCAGCCTTACTCCTTGTTGTTCGACAATACTCAg                                    <  2:1121070/150‑1 (MQ=255)
                                                                                                              tGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGt                                   >  1:991901/1‑150 (MQ=255)
                                                                                                                attattGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTc                                 <  1:490480/150‑1 (MQ=255)
                                                                                                                    ttGAAAACGCTGCCGATGGCGTCATTGCCATAGCCCGCCAGGGAGATGTAACCACCATGAACCCAGCAGCAGAAGTAATCACTGGCTATCAACGCCATGAACAGGTAGGGCAGCCTTACTCCATCAAGTTCGACAAAACTCACTTCTACa                             >  2:123736/1‑150 (MQ=255)
                                                                                                                     tGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAg                            >  2:1293482/1‑150 (MQ=255)
                                                                                                                     tGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAg                            >  2:109280/1‑150 (MQ=255)
                                                                                                                     tGAAAACGCTGCAGATGTCGTCATTGCCATTGACCGCAAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAg                            >  2:24115/1‑150 (MQ=255)
                                                                                                                        aaaCGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAGTcc                         <  1:1090984/150‑1 (MQ=255)
                                                                                                                             cTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAGTCCAGTAc                    <  1:578818/150‑1 (MQ=255)
                                                                                                                             cTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCACCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAGTCCAGTAc                    <  1:1142592/150‑1 (MQ=255)
                                                                                                                              tGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAGTCCAGTACt                   >  1:921268/1‑150 (MQ=255)
                                                                                                                               gCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAGTCCAGTACTg                  >  1:840747/1‑150 (MQ=255)
                                                                                                                                  cctggcggTCAATGGCAATGACGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAGTCCAGTACTGGAt               >  1:206980/8‑150 (MQ=255)
                                                                                                                                   ctggcggTCAATGGCAATGACGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAGTCCAGTACTGGATa              >  1:980952/7‑150 (MQ=255)
                                                                                                                                   aTGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAGTCCAGTACTGGATa              <  2:991901/150‑1 (MQ=255)
                                                                                                                                   aTGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAGTCCAGTACTGGATa              <  1:17692/150‑1 (MQ=255)
                                                                                                                                    tGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAGTCCAGTACTGGATAc             <  2:507280/150‑1 (MQ=255)
                                                                                                                                         gTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTg         <  2:1502077/150‑1 (MQ=255)
                                                                                                                                               gCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACAt   >  1:414834/1‑150 (MQ=255)
                                                                                                                                                ccATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATg  >  2:50647/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
ACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAATAACCTCGCCCAGGCACTGCGTGAAACGCGGACACTTAACGATCTGATTATTGAAAACGCTGCCGATGGCGTCATTGCCATTGACCGCCAGGGTGATGTAACCACCATGAACCCAGCAGCAGAAGTTATCACTGGCTATCAACGCCATGAACTGGTAGGGCAGCCTTACTCCATGTTGTTCGACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATG  >  NC_000913/2320729‑2321021

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: