Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,879,137 A→C 7.6% P212P (CCT→CCG) ‡ cas1 ← multifunctional nuclease Cas1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,879,1370AC 7.6% 310.2 / 16.0 106P212P (CCT→CCG) ‡cas1multifunctional nuclease Cas1
Reads supporting (aligned to +/- strand):  ref base A (46/52);  new base C (5/3);  total (51/55)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

Only 100 of 108 total aligned reads displayed.

GCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGAC  >  NC_000913/2878988‑2879285
                                                                                                                                                     |                                                                                                                                                    
gCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAaga                                                                                                                                                      <  2:741099/150‑1 (MQ=255)
   aaTGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCTTAAAAAAAGGAAAGAGGc                                                                                                                                                   <  2:1629983/150‑1 (MQ=255)
    aTGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCt                                                                                                                                                  <  2:517438/150‑1 (MQ=255)
     tGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCtt                                                                                                                                                 >  1:645467/1‑150 (MQ=255)
        tttACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTcc                                                                                                                                              <  1:1144926/150‑1 (MQ=255)
          tactacTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTg                                                                                                                                            <  1:244887/150‑1 (MQ=255)
            ctactGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGCCTCTTTCCtttg                                                                                                                                          <  1:1132672/150‑4 (MQ=255)
            ctactGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTa                                                                                                                                          <  2:814463/150‑1 (MQ=255)
            ctactGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTa                                                                                                                                          <  2:464951/150‑1 (MQ=255)
            ctactGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTa                                                                                                                                          <  1:292027/150‑1 (MQ=255)
            ctactGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTa                                                                                                                                          <  2:1463919/150‑1 (MQ=255)
            atactGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTa                                                                                                                                          <  1:401277/149‑1 (MQ=255)
              actGCGAAAAATACCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAACCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGCCTCTTTCCtttgtt                                                                                                                                        <  1:1633262/150‑6 (MQ=255)
               cttcGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGc                                                                                                                                       <  2:862780/147‑1 (MQ=255)
                tGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGCCTCTTTCCtttgttta                                                                                                                                      <  1:829799/150‑8 (MQ=255)
                 gCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGcac                                                                                                                                     >  1:1114902/1‑150 (MQ=255)
                 gCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGCTCCGGCTCACGAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGcac                                                                                                                                     <  2:886853/150‑1 (MQ=255)
                  cGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGcaca                                                                                                                                    >  1:786703/1‑150 (MQ=255)
                   gAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACaa                                                                                                                                   >  1:573444/1‑150 (MQ=255)
                      aaaTATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAAcc                                                                                                                                <  2:1001307/150‑1 (MQ=255)
                            cccTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAg                                                                                                                          <  2:474841/150‑1 (MQ=255)
                            cccTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTATCCTGTATGCACAAACCCAATAg                                                                                                                          >  2:878600/1‑150 (MQ=255)
                                gCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTgg                                                                                                                      <  2:747032/150‑1 (MQ=255)
                                 cACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGt                                                                                                                     <  2:482770/150‑1 (MQ=255)
                                   cGCCAAACTCACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGACGCTTTCCTGTAAGCACAAACCAAATAGCTGGtcc                                                                                                                   >  2:46952/1‑148 (MQ=255)
                                   cGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAATCTTTCTGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGc                                                                                                                   <  2:415491/150‑1 (MQ=255)
                                   cGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGc                                                                                                                   >  1:809064/1‑150 (MQ=255)
                                    gCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCa                                                                                                                  <  2:1172357/150‑1 (MQ=255)
                                      cAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATa                                                                                                                >  2:1431980/1‑150 (MQ=255)
                                        aaCGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAAc                                                                                                              >  2:1449437/1‑150 (MQ=255)
                                         aCGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAAcc                                                                                                             <  1:365100/150‑1 (MQ=255)
                                          cGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCa                                                                                                            <  1:796189/150‑1 (MQ=255)
                                            gACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGc                                                                                                          <  1:984736/150‑1 (MQ=255)
                                              cTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTg                                                                                                        >  2:1151663/1‑150 (MQ=255)
                                               ttCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGc                                                                                                       >  2:1239050/1‑150 (MQ=255)
                                                tCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACaa                                                                                                                                   >  2:13139/1‑121 (MQ=255)
                                                tCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACaa                                                                                                                                   <  1:13139/121‑1 (MQ=255)
                                                tCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCa                                                                                                      >  1:232250/1‑150 (MQ=255)
                                                  ccGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAg                                                                                                    >  1:37054/1‑150 (MQ=255)
                                                    ggTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTa                                                                                                  >  2:898821/1‑150 (MQ=255)
                                                     gTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTAt                                                                                                 <  2:718885/150‑1 (MQ=255)
                                                        cGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATcgc                                                                                              <  1:1431980/150‑1 (MQ=255)
                                                        cGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATcgc                                                                                              >  1:1387586/1‑150 (MQ=255)
                                                             cACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCtt                                                                                         >  2:608777/1‑150 (MQ=255)
                                                               ccAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCa                                                                                       >  2:97307/1‑150 (MQ=255)
                                                               ccAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCa                                                                                       >  1:1147011/1‑150 (MQ=255)
                                                                cAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAg                                                                                      >  1:1384714/1‑150 (MQ=255)
                                                                cAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAg                                                                                      <  2:1342931/150‑1 (MQ=255)
                                                                 aGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGt                                                                                     <  2:1081357/150‑1 (MQ=255)
                                                                   ggTTACGACGCGCTATCTCAAAAGCTTTCGTTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTa                                                                                   <  1:1039542/150‑1 (MQ=255)
                                                                      tACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGc                                                                                <  1:1239050/150‑1 (MQ=255)
                                                                        cgacgCGCTATCTCAAAAGCTTTCGGTACTACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCg                                                                              >  2:109921/1‑150 (MQ=255)
                                                                         gacgCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGt                                                                             >  1:71102/1‑150 (MQ=255)
                                                                          acgCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGta                                                                            <  1:951517/150‑1 (MQ=255)
                                                                           cgcgCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGtat                                                                           <  2:476703/150‑1 (MQ=255)
                                                                            gcgcTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCGGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGtata                                                                          <  2:1429401/150‑1 (MQ=255)
                                                                              gcTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATaaa                                                                        >  1:134838/1‑150 (MQ=255)
                                                                              gcTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATaaa                                                                        >  2:229708/1‑150 (MQ=255)
                                                                                tATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATg                                                                                                                                        >  1:391595/1‑84 (MQ=255)
                                                                                tATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATg                                                                                                                                        <  2:391595/84‑1 (MQ=255)
                                                                                    tcAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGaa                                                                  <  1:191363/150‑1 (MQ=255)
                                                                                    tcAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGaa                                                                  >  1:654012/1‑150 (MQ=255)
                                                                                     cAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAg                                                                 >  2:1209217/1‑150 (MQ=255)
                                                                                      aaaaGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGt                                                                >  2:856969/1‑150 (MQ=255)
                                                                                       aaaGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGtt                                                               >  2:171630/1‑150 (MQ=255)
                                                                                            tttCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAgc                                                          <  2:809064/150‑1 (MQ=255)
                                                                                            tttCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAgc                                                          <  2:1387586/150‑1 (MQ=255)
                                                                                               cGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCt                                                       >  2:1413271/1‑150 (MQ=255)
                                                                                                ggTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTa                                                      >  2:857195/1‑150 (MQ=255)
                                                                                                ggTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTa                                                      <  1:637550/150‑1 (MQ=255)
                                                                                                  tACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAAt                                                    >  2:548009/1‑150 (MQ=255)
                                                                                                   acaacaGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATg                                                   <  1:1099562/150‑1 (MQ=255)
                                                                                                      acaGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCAt                                                >  2:297074/1‑150 (MQ=255)
                                                                                                         gtgtCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTgg                                             <  2:343209/150‑1 (MQ=255)
                                                                                                          tgtCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGt                                            <  1:97307/150‑1 (MQ=255)
                                                                                                             cAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGa                                         >  2:322966/1‑150 (MQ=255)
                                                                                                              aaaTTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGAt                                        <  1:856969/150‑1 (MQ=255)
                                                                                                                        atgTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGccc                              <  1:836874/150‑1 (MQ=255)
                                                                                                                               cAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTccc                       <  2:37054/150‑1 (MQ=255)
                                                                                                                                   atCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTc                   >  2:1627982/1‑150 (MQ=255)
                                                                                                                                     cGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCtt                 <  2:1384714/150‑1 (MQ=255)
                                                                                                                                      gTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCttt                >  1:988626/1‑150 (MQ=255)
                                                                                                                                       tgcatacaGGAAAGCCTCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACAAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGAAGCGCTAATGAATTGGTTGATCGTATAGCCATTTTACCAGTCTTTc               >  2:834336/8‑150 (MQ=255)
                                                                                                                                       tAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTc               >  2:1503584/1‑150 (MQ=255)
                                                                                                                                        aaaCAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCGCTTTTCCCAGTCTTTCg              <  2:982787/150‑1 (MQ=255)
                                                                                                                                         catacaGGAAAGCCTCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCGGTCTTTCgg             >  1:501308/6‑150 (MQ=255)
                                                                                                                                         aaCAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCgg             <  2:134838/150‑1 (MQ=255)
                                                                                                                                         aaCAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCgg             <  2:1293385/150‑1 (MQ=255)
                                                                                                                                         aaCAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCgg             >  2:860520/1‑150 (MQ=255)
                                                                                                                                          atacaGGAAAGCCTCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGa            >  1:444413/5‑150 (MQ=255)
                                                                                                                                          aCAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCTCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCATTCTTTCGGa            <  2:573444/150‑1 (MQ=255)
                                                                                                                                           tacaGGAAAGCCTCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGTCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGAt           >  2:696899/4‑150 (MQ=255)
                                                                                                                                           cAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGAt           <  2:1147011/150‑1 (MQ=255)
                                                                                                                                            acaGGAAAGCCTCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATc          >  2:313846/3‑150 (MQ=255)
                                                                                                                                            aaaGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATc          <  1:857195/150‑1 (MQ=255)
                                                                                                                                             aaGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCg         >  2:1288785/1‑150 (MQ=255)
                                                                                                                                             aaGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCATTCGGATCg         >  2:1186130/1‑150 (MQ=255)
                                                                                                                                               ggAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTa       >  1:1164606/1‑150 (MQ=255)
                                                                                                                                                  aaGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCg    >  2:1215322/1‑150 (MQ=255)
                                                                                                                                                    gagGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGAc  <  2:654012/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGAC  >  NC_000913/2878988‑2879285

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