Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,630,921 A→T 18.5% intergenic (+56/+37) yfgI → / ← guaA nalidixic acid resistance protein, putative periplasmic protein/GMP synthetase (glutamine aminotransferase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,630,9210AT18.5% 220.7 / 38.7 104intergenic (+56/+37)yfgI/guaAnalidixic acid resistance protein, putative periplasmic protein/GMP synthetase (glutamine aminotransferase)
Reads supporting (aligned to +/- strand):  ref base A (49/34);  new base T (8/11);  total (57/47)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

GACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTAC  >  NC_000913/2630773‑2631069
                                                                                                                                                     |                                                                                                                                                    
gACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTaa                                                                                                                                                     >  2:1149281/1‑150 (MQ=255)
  cGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACa                                                                                                                                                   <  2:1513811/150‑1 (MQ=255)
  cGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGAAGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACa                                                                                                                                                   <  2:907431/150‑1 (MQ=255)
    ccAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCATATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACaga                                                                                                                                                 >  1:1142464/1‑150 (MQ=255)
        cAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAATCAGAAGAGAACAAATCGGACAGCTTGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TTACTACAGAGGGt                                                                                                                                             <  1:369892/150‑1 (MQ=255)
        cAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TTACTACAGAGGGt                                                                                                                                             <  2:65775/150‑1 (MQ=255)
        cAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TTACTACAGAGGGt                                                                                                                                             <  1:222482/150‑1 (MQ=255)
        cAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TTACTACAGAGGGt                                                                                                                                             <  2:588583/150‑1 (MQ=255)
        cAGAAATATGACTGAAGGCGGGCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TTACTACAGAGGGt                                                                                                                                             <  1:203040/150‑1 (MQ=255)
            aaTATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TTACTACAGAGGGttttt                                                                                                                                         <  1:515226/150‑2 (MQ=255)
             atatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TTACTACAGAGGGtttttt                                                                                                                                        <  2:693626/150‑3 (MQ=255)
             atatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGt                                                                                                                                        <  2:948060/150‑1 (MQ=255)
             atatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGt                                                                                                                                        <  2:941929/150‑1 (MQ=255)
              tatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCAATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGtt                                                                                                                                       <  2:91066/150‑1 (MQ=255)
              tatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAACAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TTACTACAGAGGGtttttta                                                                                                                                       <  1:197818/150‑4 (MQ=255)
               atGTCTGAATGCGTTCTCTATCTGATGAAAAGAGTGGGCAAATCAGAAGAGATCAAATGGGAGAGGTGTTTTAATTAATCGATGTTAGTAACTTGAATCCTTTAATTCTTTAAGATAAAATACCCTCTGTTAGT‑ACAGAGGGttttttat                                                                                                                                      <  2:1551642/150‑5 (MQ=255)
                tGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGttca                                                                                                                                     >  2:678090/1‑150 (MQ=255)
                 gACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGttcat                                                                                                                                    <  2:602840/150‑1 (MQ=255)
                    tGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGttcattca                                                                                                                                 >  1:839316/1‑150 (MQ=255)
                     gAAGGCGGTCGCTATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTtcattcat                                                                                                                                <  1:168668/150‑1 (MQ=255)
                     gAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTtcattcat                                                                                                                                >  2:1270497/1‑150 (MQ=255)
                     gAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTtcattcat                                                                                                                                >  1:13579/1‑150 (MQ=255)
                           ggTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCa                                                                                                                          >  2:115367/1‑150 (MQ=255)
                            gTCGCGATCTGATGAAAAGCGTGGGCAAAGTAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAg                                                                                                                         >  2:1692660/1‑150 (MQ=255)
                            gTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAg                                                                                                                         <  1:854173/150‑1 (MQ=255)
                              cgcgATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAggg                                                                                                                       <  2:458825/150‑1 (MQ=255)
                                    cTGATGTAAAGCGTGGGCAAAGCAGAAGAGACCAAATCGTACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAAt                                                                                                                 <  2:1142464/150‑1 (MQ=255)
                                    cTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAAt                                                                                                                 >  1:1525423/1‑150 (MQ=255)
                                     tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCGATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATc                                                                                                                <  1:8590/150‑1 (MQ=255)
                                     tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATc                                                                                                                <  1:1270497/150‑1 (MQ=255)
                                     tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATc                                                                                                                >  2:190886/1‑150 (MQ=255)
                                     tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACACTCTG‑TTGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATc                                                                                                                >  1:846972/1‑150 (MQ=255)
                                      gatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCa                                                                                                               <  1:450416/150‑1 (MQ=255)
                                          aaaaGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTcc                                                                                                           >  2:717611/1‑150 (MQ=255)
                                           aaaGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAATATAAAAAACCCTCTG‑TAGTAAAAGAGGGTTTTGTTCATTCATAGTTCAGGGTCAAATCATTccc                                                                                                          <  1:738238/150‑1 (MQ=255)
                                             aGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCAc                                                                                                        <  1:1180854/150‑1 (MQ=255)
                                                     cAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGt                                                                                                >  2:113301/1‑150 (MQ=255)
                                                      aaaGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTa                                                                                               <  1:471005/150‑1 (MQ=255)
                                                        aGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAg                                                                                              <  1:573670/149‑1 (MQ=255)
                                                        aGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAg                                                                                              >  2:573670/1‑149 (MQ=255)
                                                          cAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACACTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTg                                                                                           >  1:1214951/1‑150 (MQ=255)
                                                           agaagaGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTgg                                                                                          >  2:187708/1‑150 (MQ=255)
                                                              agagAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAAAAGAGGGTTTTTTTCATTCATATTGAAGGGAAAAATAATTCACAATAAATGGTAGATtgcgg                                                                                       >  1:18529/1‑150 (MQ=255)
                                                                agaACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCt                                                                                     >  2:43229/1‑150 (MQ=255)
                                                                     acaTCGGACAGCTGGTTTACTTACTCGATGGTAGTCACTTCACTCCTCGACTGCTTGCAGATACACAACCCTCTG‑TCGTCCCAGAGGGTTTTGTGCATTCCTAGTGCCGGGGCACATCATTCCCACTCATTGGTAGCTGGCGGCTTGCCg                                                                                <  1:215481/148‑1 (MQ=255)
                                                                       aTCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAATTCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCt                                                                              <  1:115367/150‑1 (MQ=255)
                                                                       aTCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCt                                                                              <  2:72567/150‑1 (MQ=255)
                                                                          ggACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGAt                                                                           >  1:838221/1‑150 (MQ=255)
                                                                               gCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCAt                                                                      <  1:164486/150‑1 (MQ=255)
                                                                               gCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCAt                                                                      >  2:684093/1‑150 (MQ=255)
                                                                                cTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATa                                                                     <  2:65164/150‑1 (MQ=255)
                                                                                cTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGATGGCGGCTTGCCGCTGATGTCATa                                                                     >  2:38397/1‑150 (MQ=255)
                                                                                       aattaatCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCAcg                                                              >  1:1273478/1‑150 (MQ=255)
                                                                                           aatCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCggg                                                          <  2:839316/150‑1 (MQ=255)
                                                                                           aatCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGATGATGTCATACAACACGCggg                                                          >  1:881399/1‑150 (MQ=255)
                                                                                              cGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGaaa                                                       >  2:1022731/1‑150 (MQ=255)
                                                                                                aTGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATa                                                     >  1:1345381/1‑150 (MQ=255)
                                                                                                   ttAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCg                                                  >  1:1248083/1‑150 (MQ=255)
                                                                                                   ttAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCg                                                  >  1:1260056/1‑150 (MQ=255)
                                                                                                         aCTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACt                                            >  1:811138/1‑150 (MQ=255)
                                                                                                         aCTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCATGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACt                                            >  1:640988/1‑150 (MQ=255)
                                                                                                          cTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACAACGCGGGAAATACAGTTCACtt                                           >  1:1032572/1‑150 (MQ=255)
                                                                                                            tCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGGCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATTCACCACGCGGGAAATACCGTTCACTTCa                                         <  1:49736/150‑1 (MQ=255)
                                                                                                              aaTCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCAtt                                       <  1:981950/150‑1 (MQ=255)
                                                                                                              aaTCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCAtt                                       <  2:1463492/150‑1 (MQ=255)
                                                                                                               aTCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTg                                      <  1:574542/150‑1 (MQ=255)
                                                                                                                tCCTATAATTCTTTAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCAACTCAATGGTAGCTGGAGGATTGCAACTTATGTCATACAACACGCGGGAAATACCGTTCACTTCATTGa                                     >  1:1694518/1‑150 (MQ=255)
                                                                                                                tCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGa                                     <  1:1174041/150‑1 (MQ=255)
                                                                                                                 ccTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGAt                                    <  1:389502/150‑1 (MQ=255)
                                                                                                                 ccTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGAt                                    >  1:1613206/1‑150 (MQ=255)
                                                                                                                     taATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGTAAATACCGTTCACTTCATTGATAATg                                <  2:1273478/150‑1 (MQ=255)
                                                                                                                       aTTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCg                              >  1:1294205/1‑150 (MQ=255)
                                                                                                                       aTTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCg                              >  1:979477/1‑150 (MQ=255)
                                                                                                                        ttCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCgg                             <  2:1525423/150‑1 (MQ=255)
                                                                                                                          cTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGtt                           >  1:617449/1‑150 (MQ=255)
                                                                                                                          cTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGtt                           >  2:180867/1‑150 (MQ=255)
                                                                                                                               atgaacaaaaCCCTCTG‑TTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATTCCGTTCACTTCATTGATAATGCAGTTGGaaa                      >  2:850545/7‑150 (MQ=255)
                                                                                                                               aGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGaaa                      >  1:1163103/1‑150 (MQ=255)
                                                                                                                               aGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGaaa                      <  1:602226/150‑1 (MQ=255)
                                                                                                                                gATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAc                     >  2:204333/1‑150 (MQ=255)
                                                                                                                                  tAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAcgc                   >  1:1335942/1‑150 (MQ=255)
                                                                                                                                   aaaaaaCCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAcgcg                  >  2:1316130/1‑150 (MQ=255)
                                                                                                                                   aaaaaaCCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTAATAGTGCAGGGTCAAAACATTCACACTCAATCGTAGCTGGAGGATTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAcgcg                  >  1:64529/1‑150 (MQ=255)
                                                                                                                                    caaaaCCCTCTG‑TTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGa                 >  1:22684/2‑150 (MQ=255)
                                                                                                                                    aaaaaCCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTCGGAAACGCGa                 >  1:1001790/1‑150 (MQ=255)
                                                                                                                                       aaCCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCg              >  1:1479898/1‑150 (MQ=255)
                                                                                                                                        aCCCTCTG‑TTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTg                          >  1:691457/1‑137 (MQ=255)
                                                                                                                                        aCCCTCTG‑TTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGa             >  1:957798/1‑150 (MQ=255)
                                                                                                                                        aCCCTCTG‑TTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGa             >  1:411973/1‑150 (MQ=255)
                                                                                                                                        aCCCTCTG‑TTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGa             >  1:515587/1‑150 (MQ=255)
                                                                                                                                        aCCCTCTG‑TTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGa             >  2:331268/1‑150 (MQ=255)
                                                                                                                                        aCCCTCTG‑TTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACACAATGCTAGCTGGCTTCTTGCCGAATATGTCATACACAACGCGGGAAATAACTTTCACATCATTGATAATGCGGTTGGAAACGCGACCGa             >  1:283760/1‑150 (MQ=255)
                                                                                                                                          ccTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAgg           <  1:987446/150‑1 (MQ=255)
                                                                                                                                           ctctG‑TAGTAACAGAGGGTTTTGTTGATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGAAGCTGATGTCATACACCACCCGGGAAATTCCGTTCACTTCATTGATAATGCGGTTGGAACCGCGACCGAGGa          >  2:300335/1‑150 (MQ=255)
                                                                                                                                               g‑tagtaACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATc      >  1:1162528/1‑150 (MQ=255)
                                                                                                                                                 tagtaACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCg     <  1:1149281/150‑1 (MQ=255)
                                                                                                                                                   gtaACAGAGGGTTTTGTTCATTCATAGTGAAGGCTCAAATCATTCCCACTCAATGGTAGCTGGCGGCAAGCCGGTTATGTCATACACCACGCGGGAAATATCGTTCACTTCATTGATAATGCGCTAGGAAACGCGACCGAGGAAATCGTa   >  2:536106/1‑150 (MQ=255)
                                                                                                                                                    taACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAc                     >  2:1390196/1‑131 (MQ=255)
                                                                                                                                                    taACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAc                     <  1:1390196/131‑1 (MQ=255)
                                                                                                                                                    taACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTAc  <  2:13579/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACGCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTG‑TAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTAC  >  NC_000913/2630773‑2631069

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