Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 317,143 C→A 20.9% pseudogene (27/684 nt) ykgA ← pseudogene, AraC family;putative regulator; Not classified; putative ARAC‑type regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913317,1430CA20.9% 327.1 / 65.9 119pseudogene (27/684 nt)ykgApseudogene, AraC family;putative regulator; Not classified; putative ARAC‑type regulatory protein
Reads supporting (aligned to +/- strand):  ref base C (45/49);  new base A (11/14);  total (56/63)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 121 total aligned reads displayed.

CATGCATGAAATTTCGGAACAGAAGCTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATAGAGAGGTTTCAGTCCTACATTATTAATGA  >  NC_000913/316995‑317291
                                                                                                                                                    |                                                                                                                                                    
catgcatgAAATTTCGGAACAGAAGCTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATtc                                                                                                                                                      <  4:781553/149‑1 (MQ=255)
      tgAAATTTCGGAACAGAAGCTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGtt                                                                                                                                                >  4:708812/1‑149 (MQ=255)
        aaaTTTCGGAACAGAAGCTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTaa                                                                                                                                              <  1:431752/149‑1 (MQ=255)
         aaTTTCGGAACAGAAGCTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAAt                                                                                                                                             >  4:657697/1‑149 (MQ=255)
          aTTTCGGAACAGAAGCTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCGTTCTCCGTTAAtt                                                                                                                                            <  1:389892/149‑1 (MQ=255)
          aTTTCGGAACAGAAGCTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAAtt                                                                                                                                            >  4:750333/1‑149 (MQ=255)
            ttCGGAACAGAAGCTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCt                                                                                                                                          <  3:57914/149‑1 (MQ=255)
             tCGGAACAGAAGCTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCtt                                                                                                                                         <  2:377699/149‑1 (MQ=255)
              cGGAACAGAAGCTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCttt                                                                                                                                        <  1:701355/149‑1 (MQ=255)
              cGGAACAGAAGCTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTATCCGTTAATTCttt                                                                                                                                        <  4:584286/149‑1 (MQ=255)
              cGGAAAAGAAGCTGGATGTTGCGGCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCttt                                                                                                                                        <  4:34397/149‑1 (MQ=255)
                 aaCAGAAGCTGAATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTATCCGTTAATTCTTTGTc                                                                                                                                     >  2:150259/1‑149 (MQ=255)
                  aCAGAAGCTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCa                                                                                                                                    <  3:708812/149‑1 (MQ=255)
                  aCAGAAGCTGAATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTATCCGTTAATTCTTTGTCa                                                                                                                                    <  1:150259/149‑1 (MQ=255)
                    aGAAGCTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCata                                                                                                                                  <  1:432584/149‑1 (MQ=255)
                      aaGCTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATc                                                                                                                                <  4:677649/149‑1 (MQ=255)
                      aaGCTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGCGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATc                                                                                                                                <  2:36450/149‑1 (MQ=255)
                        gCTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGc                                                                                                                              <  2:663837/149‑1 (MQ=255)
                         cTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTATCCGTTAATTCTTTGTCATATCGca                                                                                                                             <  2:21637/149‑1 (MQ=255)
                            gATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACAt                                                                                                                          <  4:360980/149‑1 (MQ=255)
                               gTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACAtttt                                                                                                                       >  4:563339/1‑149 (MQ=255)
                                     tCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTtat                                                                                                                 <  3:657697/149‑1 (MQ=255)
                                     tCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTtat                                                                                                                 <  1:113280/149‑1 (MQ=255)
                                         ctgTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACAt                                                                                                                          <  3:268911/136‑1 (MQ=255)
                                              gcgaGATTTCTGTGCGATATCTTCGGTTGAAATAGGGTGCTCAAGATTGCAGTCACTCCACTCCAGGAGCTTCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGc                                                                                                        >  4:575114/4‑149 (MQ=255)
                                              gCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTc                                                                                                                      <  1:699394/135‑1 (MQ=255)
                                              gCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTc                                                                                                                      >  2:699394/1‑135 (MQ=255)
                                                cAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAAt                                                                                                                                             >  1:74361/1‑110 (MQ=255)
                                                cAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAAt                                                                                                                                             <  2:74361/110‑1 (MQ=255)
                                                   aTTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTaa                                                                                                   >  2:94013/1‑149 (MQ=255)
                                                    tttCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGt                                                                                                                                      <  2:525125/113‑1 (MQ=255)
                                                    tttCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGt                                                                                                                                      >  1:525125/1‑113 (MQ=255)
                                                      tCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAgtg                                                                                                <  4:44522/149‑1 (MQ=255)
                                                        tgtgCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTa                                                                                              >  4:426381/1‑149 (MQ=255)
                                                        tgtgCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTa                                                                                              >  3:106501/1‑149 (MQ=255)
                                                           gCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCa                                                                                           >  4:551336/1‑149 (MQ=255)
                                                            cGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTACCTGATCATCTCATTATCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCaa                                                                                          <  2:890931/149‑1 (MQ=255)
                                                             gATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCaaa                                                                                         <  2:741671/149‑1 (MQ=255)
                                                                atCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTATCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAAt                                                                                      <  4:814773/149‑1 (MQ=255)
                                                                 tCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATa                                                                                     <  1:680106/149‑1 (MQ=255)
                                                                   ttCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTGCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGt                                                                                                                                      <  2:511789/98‑1 (MQ=255)
                                                                   ttCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTGCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGt                                                                                                                                      >  1:511789/1‑98 (MQ=255)
                                                                   ttCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGc                                                                                   >  3:261370/1‑149 (MQ=255)
                                                                      gATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAg                                                                                <  2:795161/149‑1 (MQ=255)
                                                                       aTTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAgg                                                                               >  4:36103/1‑149 (MQ=255)
                                                                       aTTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTATCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAgg                                                                               >  2:17207/1‑149 (MQ=255)
                                                                         tGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGc                                                                             <  1:400358/149‑1 (MQ=255)
                                                                          gAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCg                                                                            <  1:391357/149‑1 (MQ=255)
                                                                            aaTAGGGTGCTCAAGATTGCACTCAATCCACTCCACGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATCTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGaa                                                                          <  3:632470/149‑1 (MQ=255)
                                                                             aTAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAg                                                                         <  3:66192/149‑1 (MQ=255)
                                                                                   tGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCaaaa                                                                   <  2:701592/149‑1 (MQ=255)
                                                                                   tGCTCAAGAGTGCACTCAATCCAATCCACGAGCTGCTGTAGAATCTTCTTCCTGATCATCTCATTATCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCGAAAATAGCAAGGGCGAAGGCaaaa                                                                   <  3:780733/149‑1 (MQ=255)
                                                                                    gCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAg                                                                  <  2:393257/149‑1 (MQ=255)
                                                                                      tCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCa                                                                >  2:314020/1‑149 (MQ=255)
                                                                                      tCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCa                                                                >  1:78117/1‑149 (MQ=255)
                                                                                        aaGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACt                                                              <  3:472960/149‑1 (MQ=255)
                                                                                          gATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTgg                                                            >  1:177102/1‑149 (MQ=255)
                                                                                           aTTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACt                                                              >  4:440606/1‑146 (MQ=255)
                                                                                           aTTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACt                                                              <  3:440606/146‑1 (MQ=255)
                                                                                                aCTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCAc                                                               >  4:213357/1‑140 (MQ=255)
                                                                                                 cTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTATCCGTTaa                                                                                                                                              >  2:568725/1‑60 (MQ=255)
                                                                                                 cTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTATCCGTTaa                                                                                                                                              <  1:568725/60‑1 (MQ=255)
                                                                                                  tCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCt                                                    >  4:208059/1‑149 (MQ=255)
                                                                                                    aaTCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGc                                                  <  1:577672/149‑1 (MQ=255)
                                                                                                         aCTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTATCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGt                                             >  3:617253/1‑149 (MQ=255)
                                                                                                             cAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTc                                         >  2:718463/1‑149 (MQ=255)
                                                                                                             cAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTATCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTc                                         <  1:143583/149‑1 (MQ=255)
                                                                                                              aGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTcc                                        <  3:467223/149‑1 (MQ=255)
                                                                                                                    tgctgTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTACTTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGaaa                                  >  4:535567/1‑149 (MQ=255)
                                                                                                                      ctgTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAAt                                >  4:622807/1‑149 (MQ=255)
                                                                                                                       tgTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATa                               <  4:106501/149‑1 (MQ=255)
                                                                                                                        gTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATag                              >  2:140766/1‑149 (MQ=255)
                                                                                                                         tAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATaga                             <  3:750333/149‑1 (MQ=255)
                                                                                                                         tAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATaga                             <  1:94013/149‑1 (MQ=255)
                                                                                                                         tAGAATCTTCTGCCTGATCATCTCATTATCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATaga                             <  2:623803/149‑1 (MQ=255)
                                                                                                                            aaTCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATagagag                          >  4:303980/1‑149 (MQ=255)
                                                                                                                            aaTCTTCTGCCTGATCATCTCATTATCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCGTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATagagag                          <  2:786261/149‑1 (MQ=255)
                                                                                                                            aaTCTTCTGCCTGATCATCTCATTATCCGTTAATTCGTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATagagag                          >  2:871623/1‑149 (MQ=255)
                                                                                                                              tcttctGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATAGAGAGGt                        >  1:767293/1‑149 (MQ=255)
                                                                                                                              tcttctGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATAGAGAGGt                        >  2:486713/1‑149 (MQ=255)
                                                                                                                               cttcttCCTGATCATCTCAATCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATAGAGAGGtt                       <  1:792098/149‑1 (MQ=255)
                                                                                                                                  ctGCCTGATCATCTCATTATCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATAGAGAGGTTTCa                    <  1:898574/149‑1 (MQ=255)
                                                                                                                                    gCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATAGAGAGGTTTCAGt                  >  1:145132/1‑149 (MQ=255)
                                                                                                                                      cTGATCATCTCATTATCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATAGAGAGGTTTCAGTcc                >  4:162548/1‑149 (MQ=255)
                                                                                                                                       tGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATAGAGAGGTTCCAGTCCt               >  3:373891/1‑149 (MQ=255)
                                                                                                                                       tGATCATCTCATTATCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATAGAGAGGTTTCAGTCCt               >  4:465622/1‑149 (MQ=255)
                                                                                                                                         atcatcTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATAGAGAGGTTTCAGTCCTAc             >  3:82223/1‑149 (MQ=255)
                                                                                                                                             tctcATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAg                                                                  <  3:435650/92‑1 (MQ=255)
                                                                                                                                             tctcATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAg                                                                  >  4:435650/1‑92 (MQ=255)
                                                                                                                                             tctcATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATAGAGAGGTTTCAGTCCTACatta         >  1:592403/1‑149 (MQ=255)
                                                                                                                                              ctcATTATCCGTTAATTCTTTGTCATATCGCACATTTTCGTTTTACCCTGGGCCCTAAGTGTATCAAAAATa                                                                                     <  2:480702/72‑1 (MQ=255)
                                                                                                                                              ctcATTATCCGTTAATTCTTTGTCATATCGCACATTTTCGTTTTACCCTGGGCCCTAAGTGTATCAAAAATa                                                                                     >  1:480702/1‑72 (MQ=255)
                                                                                                                                               tcATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTccc                                       >  3:821478/1‑117 (MQ=255)
                                                                                                                                               tcATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTccc                                       <  4:821478/117‑1 (MQ=255)
                                                                                                                                               tcATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATAGAGAGGTTTCAGTCCTACattatt       >  1:621928/1‑149 (MQ=255)
                                                                                                                                               tcATTATCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGc                                                                       >  4:493832/1‑85 (MQ=255)
                                                                                                                                               tcATTATCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGc                                                                       <  3:493832/85‑1 (MQ=255)
                                                                                                                                                  ttCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGaa                                                                          <  3:201128/79‑1 (MQ=255)
                                                                                                                                                  ttCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGaa                                                                          >  4:201128/1‑79 (MQ=255)
                                                                                                                                                    ctcCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATAGAGAGGTTTCAGTCCTACATTATTAATGa  <  3:551336/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATGCATGAAATTTCGGAACAGAAGCTGGATGTTGCGTCTGCTGTAGCCAGATTTCTGTGCGATATCTTCGATTGAAATAGGGTGCTCAAGATTGCACTCAATCCACTCCAGGAGCTGCTGTAGAATCTTCTGCCTGATCATCTCATTCTCCGTTAATTCTTTGTCATATCGCACATTTTCGTTATACCCTGGGCCCTAAGTGTATCAAAAATAGCAAGGGCGAAGGCAAAAGCACTGGTAGCGTCTGCTTCGTTGTCCCCGAAAATAGAGAGGTTTCAGTCCTACATTATTAATGA  >  NC_000913/316995‑317291

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