Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 873,380 A→T 6.2% S134S (TCA→TCT)  yliE → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913873,3800AT 6.2% 1051.4 / 5.3 343S134S (TCA→TCT) yliEputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (171/135);  new base T (2/32);  total (173/170)
Fisher's exact test for biased strand distribution p-value = 5.34e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 343 total aligned reads displayed.

TGATTCTTAAAATGGAAATGATGCCATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCA  >  NC_000913/873232‑873528
                                                                                                                                                    |                                                                                                                                                    
  aTTCTTAAAATGGAAATGATGCCATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAggg                                                                                                                                                   <  1:7028463/150‑1 (MQ=255)
   ttCTTAAAATGGAAATGATGCCATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGc                                                                                                                                                  <  1:1236362/150‑1 (MQ=255)
    tCTTAAAATGGAAATGATGCCATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGACGACGTGTCAGGGCg                                                                                                                                                 >  2:6811063/1‑150 (MQ=255)
     cTTAAAATGGAAATGATGCCATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCgg                                                                                                                                                <  1:3488484/150‑1 (MQ=255)
       tAAAATGGAAATGATGCCATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGt                                                                                                                                              <  2:3911859/150‑1 (MQ=255)
          aaTGGAAATGATGCCATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTAt                                                                                                                                           >  2:775001/1‑150 (MQ=255)
            tGGAAATGATGCCATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTAttt                                                                                                                                         >  1:6100351/1‑150 (MQ=255)
            tGGAAATGATGCCATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTAttt                                                                                                                                         >  1:931338/1‑150 (MQ=255)
              gAAATGATGCCATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTAttttt                                                                                                                                       <  2:7472882/150‑1 (MQ=255)
              gAAATGATGCCATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTAttttt                                                                                                                                       <  1:3907028/150‑1 (MQ=255)
              gAAATGATGCCATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTAttttt                                                                                                                                       <  1:1352331/150‑1 (MQ=255)
              gAAATGATGCCATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTAttttt                                                                                                                                       >  2:5851015/1‑150 (MQ=255)
                aaTGATGCCATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGc                                                                                                                                     <  1:3791739/150‑1 (MQ=255)
                   gatgCCATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGa                                                                                                                                  <  2:7538990/150‑1 (MQ=255)
                      gCCATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGAttt                                                                                                                               >  1:3131472/1‑150 (MQ=255)
                        cATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTa                                                                                                                             <  2:6030058/150‑1 (MQ=255)
                         aTTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTaa                                                                                                                            <  2:7710705/150‑1 (MQ=255)
                         aTTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTaa                                                                                                                            >  2:4550691/1‑150 (MQ=255)
                         aTTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTaa                                                                                                                            >  2:4684765/1‑150 (MQ=255)
                         aTTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTaa                                                                                                                            <  2:4865742/150‑1 (MQ=255)
                          tttttAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAAc                                                                                                                           <  2:1365619/150‑1 (MQ=255)
                           ttttAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAAcc                                                                                                                          <  2:5564781/150‑1 (MQ=255)
                               aGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATc                                                                                                                      <  2:3390941/150‑1 (MQ=255)
                               aGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATc                                                                                                                      <  1:4582835/150‑1 (MQ=255)
                                gACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCg                                                                                                                     <  2:4437129/150‑1 (MQ=255)
                                   tCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGaaa                                                                                                                  >  2:2570333/1‑150 (MQ=255)
                                     ggTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACg                                                                                                                >  1:5755649/1‑150 (MQ=255)
                                     ggTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACg                                                                                                                >  2:5705342/1‑150 (MQ=255)
                                      gTGGGACTGTTAATTGATGATAATAAATATTATATTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACATGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGa                                                                                                               <  1:8095429/150‑1 (MQ=255)
                                      gTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGa                                                                                                               >  2:1508842/1‑150 (MQ=255)
                                             tGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGt                                                                                                        >  2:5064444/1‑150 (MQ=255)
                                             tGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGt                                                                                                        >  2:802171/1‑150 (MQ=255)
                                             tGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGt                                                                                                        >  2:2715542/1‑150 (MQ=255)
                                             tGGTAATTGATTATAATAAATATTATCTTTTTTCGCGGTGGGCTAATGATAAAATCGTTGTTTATAATCAGGAACAATTAAATGGACCGCTTGTCGACGAGTCAGGGCGGTTTATTTTTGCCGATTTTAACCAATCGAAACGACCGTGGt                                                                                                        <  2:6475783/150‑1 (MQ=255)
                                               gTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAGTCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCg                                                                                                      >  1:109068/1‑150 (MQ=255)
                                               gTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCg                                                                                                      >  1:518768/1‑150 (MQ=255)
                                               gTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCg                                                                                                      >  2:1401196/1‑150 (MQ=255)
                                                tACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCCGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCgg                                                                                                     >  1:4661528/1‑150 (MQ=255)
                                                   ttGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCAATCGAATCGAGCGTGGTCGGTgg                                                                                                  >  2:4941923/1‑150 (MQ=255)
                                                   ttGATGATAATAAATATTATCTATTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACAGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATATTAACCCATCGAAACGACCGTGGTCGGTgg                                                                                                  >  2:3133453/1‑150 (MQ=255)
                                                       tgatAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTc                                                                                              >  1:4195640/1‑150 (MQ=255)
                                                        tataataAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCa                                                                                             <  2:3268071/149‑1 (MQ=255)
                                                        gatAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCa                                                                                             >  1:7390839/1‑150 (MQ=255)
                                                          taataaATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGa                                                                                           <  2:4055423/150‑1 (MQ=255)
                                                          taataaATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGa                                                                                           <  2:3192920/150‑1 (MQ=255)
                                                            ataaataTTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATg                                                                                         >  2:1079207/1‑150 (MQ=255)
                                                             taaataTTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGa                                                                                        >  1:6568844/1‑150 (MQ=255)
                                                              aaataTTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGa                                                                                                                                  <  1:5603341/107‑1 (MQ=255)
                                                                atatTCTCTTTGTTGGCGGCGGGCGAAGGATAAAATGGGGGTTTAGCATCAGGAACACGTAATTGGCCCGCTGGTCGACGAGTCAGGGGGGGTTATTTTGGCCGATTTTACCCCATCGACACGACCGTGGTCGGTGGCTTCAGATGActc                                                                                     <  2:3645255/150‑1 (MQ=255)
                                                                       cTTTTTTGGGGGCGGGCGATTGATAAAATGGTTGTTTTTCATCTGGAACATGTAATGGGACGGGTGGTGGAGGTGTCTGGGGGGGTTTTTTTTGCCGATTTTAACCCTTCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTaacaac                                                                              <  2:6350647/150‑1 (MQ=255)
                                                                         tttttCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAg                                                                            >  2:8153492/1‑150 (MQ=255)
                                                                          ttttCGCGGTGGGCTATTGATAAATTGTTTTTTAATAATCTTTAACATTTAAAGGGACCGTGTGTCGTGGAGTCTGGGGGGGTTTTTTTTGCCGATTTTACCCATTCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGc                                                                           <  1:7337545/150‑1 (MQ=255)
                                                                           tttCGCGGAGGGCGAATTATAAAATCGTTTTTTATCATAAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCt                                                                          <  2:6370943/150‑1 (MQ=255)
                                                                              cgcgGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGa                                                                       >  1:6456346/1‑150 (MQ=255)
                                                                                cgGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAAt                                                                     >  2:22716/1‑150 (MQ=255)
                                                                                  gAGGGCTAATGATAACCTCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCTGGGCGGGTTTTTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATcc                                                                   <  2:1903881/150‑1 (MQ=255)
                                                                                    gggCTATTGATAAAATCGTTGTTTATCATCTGGAACATGTAAATGGACCGCTTGTCGACGAGTCAGGGGGGTTTTTTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCgg                                                                 <  1:4731715/150‑1 (MQ=255)
                                                                                    gggCTAATTTTAAAATCTTTTTTTTTCATCAGGAAAATTTAATTGGCCCGCTTTTCGACGTGTCTGGGCGGTTTATTTTTGCCTATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCgg                                                                 <  2:3556728/150‑1 (MQ=255)
                                                                                      gctaaTTATAAAATGTTTTTTTATAATATTGAAAAATTAATTGTACCGCTTGTCGAGGAGTCTGGGGGGGTTATTTTTGCCTATTTTAACCAATCTAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCa                                                               <  1:1484446/147‑1 (MQ=255)
                                                                                       cGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCAt                                                              >  2:1290013/1‑150 (MQ=255)
                                                                                       cGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACACATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCAt                                                              >  2:2239126/1‑150 (MQ=255)
                                                                                        gAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATa                                                             >  1:3354462/1‑150 (MQ=255)
                                                                                        gAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATa                                                             >  2:2041004/1‑150 (MQ=255)
                                                                                        gAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATa                                                             >  2:7427123/1‑150 (MQ=255)
                                                                                             aTAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAAtt                                                        <  2:3428647/150‑1 (MQ=255)
                                                                                              tAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCCTACAATTg                                                       >  1:7706632/1‑150 (MQ=255)
                                                                                                  aTCGTTTTTTATCATAAGTAACAAGTAAATGGCCCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCTAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCttt                                                   <  2:2851211/150‑1 (MQ=255)
                                                                                                  aTCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCttt                                                   <  2:5354110/150‑1 (MQ=255)
                                                                                                  aTCGTTGTTTAACATCAGCAACAAGTAAATGCACCGCTTGACCACGAGACAGGGCGGCTTATATTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTACAGATGACTCAAACAACAGCTGGAATCCGGCATACAATTGCttt                                                   <  2:1192946/150‑1 (MQ=255)
                                                                                                    cGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGa                                                 <  2:7211868/150‑1 (MQ=255)
                                                                                                      ttgttTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGTGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATc                                               <  2:2424154/150‑1 (MQ=255)
                                                                                                      ttgttTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCTGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATc                                               <  2:904736/150‑1 (MQ=255)
                                                                                                      ttgttTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTTCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATc                                               >  1:2438097/1‑150 (MQ=255)
                                                                                                      ttgttTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATc                                               >  2:5972600/1‑150 (MQ=255)
                                                                                                        gttTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGt                                             >  2:2803623/1‑150 (MQ=255)
                                                                                                         tttATCACAAGGAAAAATTAAATGGACCGCTTGTCGACGAGTCTGGGCGGGTTTTTTTTGCCTATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTc                                            <  1:6393876/150‑1 (MQ=255)
                                                                                                            ataatcaGGAACAAGTAAAGGGACCGCTTGTCGACGAGCCTGGGCGGGTTATTTTTGCCTATTTTAACCAATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCgg                                         <  1:3213995/147‑1 (MQ=255)
                                                                                                             tcatcaGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCggg                                        >  1:4809165/1‑150 (MQ=255)
                                                                                                              catcaGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGt                                       >  2:6650995/1‑150 (MQ=255)
                                                                                                              catcaGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGt                                       >  2:7007706/1‑150 (MQ=255)
                                                                                                                tcaGGAACAAGTAAATGGACCGCTTGTCGTCGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGTTGACTCTAACAACAGCTGGTATCCGGCATACTATTGCTTTGATCGTCCGGGTat                                     >  2:6544352/1‑149 (MQ=255)
                                                                                                                  aGGAACAAGTAAATGGACCGCTTGTCGACGTGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTaaaa                                   <  1:4487850/150‑1 (MQ=255)
                                                                                                                    taaCAAGTAATTGGACCGCTTGTCGACGTGTCTGGGCGGTTTATTTTTGCCTATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGATaaaaaa                                 <  1:5832330/149‑1 (MQ=255)
                                                                                                                     aaCAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTCAAAAAt                                >  2:4256946/1‑150 (MQ=255)
                                                                                                                     aaCAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAAt                                >  1:6013599/1‑150 (MQ=255)
                                                                                                                        aaGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTa                             >  1:1541023/1‑150 (MQ=255)
                                                                                                                          gTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTAtc                           >  1:5041936/1‑150 (MQ=255)
                                                                                                                            aaaTGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTAtctc                         <  2:7390839/150‑1 (MQ=255)
                                                                                                                             aaTGGCCCGCTTGTCGCCGCGTCTGGGCGGTTTTTTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATctct                        <  1:2155760/150‑1 (MQ=255)
                                                                                                                                ggACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCtttt                     >  2:1708470/1‑150 (MQ=255)
                                                                                                                                   ccGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACg                  <  1:6777979/150‑1 (MQ=255)
                                                                                                                                    cGCTTGTCGACGTGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGc                 <  1:4414521/150‑1 (MQ=255)
                                                                                                                                    cGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGc                 >  1:7832654/1‑150 (MQ=255)
                                                                                                                                    cGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGc                 >  2:7941751/1‑150 (MQ=255)
                                                                                                                                      cTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTa               >  2:5291173/1‑150 (MQ=255)
                                                                                                                                        tGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTaca             <  2:4775346/150‑1 (MQ=255)
                                                                                                                                                aGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACg     >  2:5653395/1‑150 (MQ=255)
                                                                                                                                                 gTCTGGGGGGGTTATTTTTGCCGATTTTAACCCTTCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACgg    <  1:7873856/150‑1 (MQ=255)
                                                                                                                                                 gTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACgg    >  2:5527620/1‑150 (MQ=255)
                                                                                                                                                   cAGGGCGGGTTATTTTTGCCGATTTTAAACCATAGAAACGACCGTGGTCGGTGGCTTCAGAAGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTATTTTACGCTACACAACAACGGCa  >  1:2547195/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGATTCTTAAAATGGAAATGATGCCATTTTTAGACTCGGTGGGACTGGTACTTGATGATAATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCA  >  NC_000913/873232‑873528

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑