Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 569,161 A→C 6.3% D87A (GAC→GCC)  ybcK → DLP12 prophage; putative recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913569,1610AC 6.3% 224.3 / 5.6 81D87A (GAC→GCC) ybcKDLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base A (36/39);  new base C (2/3);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

TTATTATCTTGATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATAC  >  NC_000913/569012‑569307
                                                                                                                                                     |                                                                                                                                                  
ttattaTCTTGATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGa                                                                                                                                                    >  1:234488/1‑150 (MQ=255)
  attatCTTGATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACa                                                                                                                                                  >  2:1096921/1‑150 (MQ=255)
    tatCTTGATACCATAACATATGAAGATTTACGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTAGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGATGAAAGTCTGGACAGa                                                                                                                                                <  2:581709/150‑1 (MQ=255)
    tatCTTGAAACCATAACATATGAAGATTTAGGATTAAGTGCAATCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATCCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGa                                                                                                                                                <  2:587046/150‑1 (MQ=255)
       cTTGATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGTTTATATATTGCCAGGAAATACATTGTTATTTGAAAGTCTGTACAGACtt                                                                                                                                             >  1:907949/1‑150 (MQ=255)
       cTTGATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACtt                                                                                                                                             >  2:108183/1‑150 (MQ=255)
        ttGATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACttt                                                                                                                                            <  2:620160/150‑1 (MQ=255)
          gATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCAGACTTTCaa                                                                                                                                         <  2:950431/150‑1 (MQ=255)
           aTACCGTAACAAAAGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTAGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaa                                                                                                                                         <  2:1133534/150‑1 (MQ=255)
           aTACCATAACATAGGCAGCTTTAGGATTAAGTGCATTCAAAGGACAGCATGCACTATCAGGAGCTTTTTAGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCaa                                                                                                                                         <  1:211489/150‑1 (MQ=255)
            tACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                                                                        <  1:333719/150‑2 (MQ=255)
            tACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                                                                        <  1:143509/150‑1 (MQ=255)
            tACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                                                                        <  1:1281248/150‑1 (MQ=255)
            aaCCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                                                                        <  1:250203/149‑1 (MQ=255)
            aaCCATAACATATGAAGATTTAGGATTAAGTGCATTAAAAGGAAAGCATGCAGAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCaac                                                                                                                                        <  1:227769/149‑2 (MQ=255)
              ccATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAagag                                                                                                                                      <  1:1151115/150‑1 (MQ=255)
              ccATAACATATGAACCTTTAGGAATAAGGGCATTCAAAGGAAAGCATGCACAATCAGGAGCTGTATCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGTAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAggag                                                                                                                                      <  2:58911/150‑1 (MQ=255)
                      tatGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTc                                                                                                                              >  2:1263924/1‑150 (MQ=255)
                        tGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTAGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCgg                                                                                                                            >  1:993642/1‑150 (MQ=255)
                        tGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGAAGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCgg                                                                                                                            >  2:369959/1‑150 (MQ=255)
                                tAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGa                                                                                                                    >  2:360723/1‑150 (MQ=255)
                                 aGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTGCGGAATTTTTAGATGCTATAGAGCATTGTTATATATTGCCAGGAACTACATTGTTATTTGAAAGTCTGGACATACTTTCAAGAGAAAAAGTCGGTGAAGCGAt                                                                                                                   <  1:809740/150‑1 (MQ=255)
                                  ggATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTGTTTCGGAAGTTTTAGATTCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGATAAAAAGTCGGTGAAGCGAtt                                                                                                                  >  1:237462/1‑150 (MQ=255)
                                   gATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGTTGAAGCGATTg                                                                                                                 >  1:395437/1‑150 (MQ=255)
                                   gATTAAGTGCATACAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAttt                                                                                                                 >  2:1245747/1‑149 (MQ=255)
                                    aTTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGa                                                                                                                >  2:27351/1‑150 (MQ=255)
                                        aGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGt                                                                                                            <  2:677880/150‑1 (MQ=255)
                                         gTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCAAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTc                                                                                                           >  2:725224/1‑150 (MQ=255)
                                           gCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTg                                                                                                         <  2:281755/150‑1 (MQ=255)
                                               tCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTAAAGAGAAAAAGTCGGTGAAGAGATTGAACGTCTGAAAt                                                                                                     >  1:175195/1‑150 (MQ=255)
                                                cAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAAtt                                                                                                    <  1:302849/150‑1 (MQ=255)
                                                    ggAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAtt                                                                                                <  1:268482/150‑1 (MQ=255)
                                                    ggAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGATATAGAGCATGGTTATATATTTACAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGAGATTGAACGTCTGAAATTGAtt                                                                                                >  2:138388/1‑150 (MQ=255)
                                                     gAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAttt                                                                                               <  2:587566/150‑1 (MQ=255)
                                                      aaaGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAtttt                                                                                              >  2:702612/1‑150 (MQ=255)
                                                       aaGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTg                                                                                             <  2:284947/150‑1 (MQ=255)
                                                        aGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGa                                                                                            >  1:610613/1‑150 (MQ=255)
                                                                 aaTCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTa                                                                                   <  1:1136000/150‑1 (MQ=255)
                                                                  aTCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTAt                                                                                  >  2:286965/1‑150 (MQ=255)
                                                                  aTCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGGCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTAt                                                                                  <  2:1173261/150‑1 (MQ=255)
                                                                    cAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTg                                                                                <  1:447185/150‑1 (MQ=255)
                                                                       gAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATg                                                                             >  2:251333/1‑150 (MQ=255)
                                                                         gCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCTGTGAAGCGATTGAACGTATGAAATTGATTTTGAATCACGGTATTGATGtt                                                                           >  1:774520/1‑150 (MQ=255)
                                                                         gCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGtt                                                                           >  2:874192/1‑150 (MQ=255)
                                                                         gCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGtt                                                                           <  1:295060/150‑1 (MQ=255)
                                                                            ttttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAATGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtata                                                                        <  2:197233/150‑1 (MQ=255)
                                                                                 ggAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTTATTTTGAATCACGGTATTGATGTTATAActct                                                                   >  2:223098/1‑150 (MQ=255)
                                                                                 ggAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAActct                                                                   >  2:509984/1‑150 (MQ=255)
                                                                                  gAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                                                                  <  1:1150853/150‑1 (MQ=255)
                                                                                   aaTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCttt                                                                 >  2:292429/1‑150 (MQ=255)
                                                                                     tttttAGATTCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGc                                                               <  2:529225/150‑1 (MQ=255)
                                                                                     tttttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGc                                                               <  1:985843/150‑1 (MQ=255)
                                                                                       tttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGa                                                             >  1:1134710/1‑150 (MQ=255)
                                                                                         tAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACa                                                           >  2:739209/1‑150 (MQ=255)
                                                                                            aTGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATa                                                        <  1:483397/150‑1 (MQ=255)
                                                                                                 ataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTAAAAAGAAAAAGTCGTTGAAGCGATTGAACGTCTGAAATTGATTTTGAAACACTGTATTGATCTTATAACACTTTGCTACAATACAGta                                                   >  1:252539/1‑149 (MQ=255)
                                                                                                   agagCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCta                                                 >  1:553743/1‑150 (MQ=255)
                                                                                                      gCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTataa                                              >  1:49868/1‑150 (MQ=255)
                                                                                                           gTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTg                                         <  2:1173171/150‑1 (MQ=255)
                                                                                                             tatataTTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGAc                                       <  2:574837/150‑1 (MQ=255)
                                                                                                               tatatTGCCAGGAACTACATTGTTCGTTGACAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGActc                                     <  1:74905/150‑1 (MQ=255)
                                                                                                                 aatTGCCAGGAACTACATTGTTAGTTGAAAGGCGGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGGCTGAAATTCATTTTGAAGCACGGTATGGATGTAATAACTCTTTGCGACACTCCCGTCTACAATATTGACTCtt                                   <  2:8408/149‑1 (MQ=255)
                                                                                                                   ttGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTg                                 <  1:847740/150‑1 (MQ=255)
                                                                                                                     gCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGaa                               >  2:105343/1‑150 (MQ=255)
                                                                                                                     gCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGaa                               >  1:816665/1‑150 (MQ=255)
                                                                                                                          taaCTACATTGTTAGTTGAAAGTCTTTACAGAATTTCAAGATAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTATGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGc                          <  1:66887/149‑1 (MQ=255)
                                                                                                                           aaCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTTAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGcc                         <  1:108183/150‑1 (MQ=255)
                                                                                                                              tagaTTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCata                      <  1:1289248/147‑1 (MQ=255)
                                                                                                                               aCATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCatat                     <  2:230553/150‑1 (MQ=255)
                                                                                                                                    gTTAGTTGAAAGTCTTGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCAtt                <  1:143883/150‑1 (MQ=255)
                                                                                                                                    gTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCAtt                <  2:1279744/150‑1 (MQ=255)
                                                                                                                                      tAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATtaa              <  1:322267/150‑1 (MQ=255)
                                                                                                                                        gTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGCCTCTTTGAATGAGCCATATTCATtaata            >  1:184224/1‑150 (MQ=255)
                                                                                                                                         ttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGTATCACGGTATTGATGTTATAACTGTTTGCGACAATACAGTCTATAATTTTGACTCTTTGAATGAGCCATTTTCATtaataa           >  1:1144358/1‑150 (MQ=255)
                                                                                                                                          tGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCTACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATaaa          >  2:149006/1‑150 (MQ=255)
                                                                                                                                          tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATaaa          <  1:286021/150‑1 (MQ=255)
                                                                                                                                             aaGTCTGGACAGAATTTCAAGAGAAAAAGTGGGTGAAGCGTTTGAACGTCTGAAATTGATTTTTAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGc       <  2:252539/150‑1 (MQ=255)
                                                                                                                                               gTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCa     >  2:1288132/1‑150 (MQ=255)
                                                                                                                                               gTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCa     >  1:1052071/1‑150 (MQ=255)
                                                                                                                                               gTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCa     >  2:159618/1‑150 (MQ=255)
                                                                                                                                               gTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTAGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCa     >  2:1249260/1‑150 (MQ=255)
                                                                                                                                                  tGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATAc  >  1:94221/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTATTATCTTGATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATAC  >  NC_000913/569012‑569307

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