Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,015,912 T→A 7.1% intergenic (+30/+40) rmf → / ← fabA ribosome modulation factor/beta‑hydroxydecanoyl thioester dehydrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,015,9120TA 7.1% 420.4 / 16.0 123intergenic (+30/+40)rmf/fabAribosome modulation factor/beta‑hydroxydecanoyl thioester dehydrase
Reads supporting (aligned to +/- strand):  ref base T (58/56);  new base A (4/5);  total (62/61)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

Only 100 of 123 total aligned reads displayed.

AGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACG  >  NC_000913/1015767‑1016059
                                                                                                                                                 |                                                                                                                                                   
aGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGc                                                                                                                                                  >  2:740668/1‑149 (MQ=255)
aGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGc                                                                                                                                                  >  1:950204/1‑149 (MQ=255)
aGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGc                                                                                                                                                  >  1:711575/1‑149 (MQ=255)
   ccGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCgg                                                                                                                                                >  4:737513/1‑148 (MQ=255)
   ccGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCgg                                                                                                                                                <  4:747057/148‑1 (MQ=255)
   ccGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGa                                                                                                                                               <  1:189790/149‑1 (MQ=255)
   ccGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGa                                                                                                                                               <  1:314246/149‑1 (MQ=255)
      gCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAgg                                                                                                                                             <  4:605265/148‑1 (MQ=255)
      gCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAgg                                                                                                                                             <  3:624767/148‑1 (MQ=255)
      gCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAgg                                                                                                                                             <  4:140042/148‑1 (MQ=255)
        aTCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCTCAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCCTTTAAAAGAAACCTCCGCATTGCGGAGGtt                                                                                                                                           <  4:648492/148‑1 (MQ=255)
        aTCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGtt                                                                                                                                           <  4:851003/148‑1 (MQ=255)
         tCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGttt                                                                                                                                          >  4:197560/1‑148 (MQ=255)
         tCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTc                                                                                                                                         >  2:463418/1‑149 (MQ=255)
          cGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAGGTTTct                                                                                                                                        <  1:783847/149‑2 (MQ=255)
               gACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCttt                                                                                                                                   >  1:1018119/1‑149 (MQ=255)
                aCGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCttt                                                                                                                                   <  4:160659/148‑1 (MQ=255)
                   cTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGa                                                                                                                                >  3:1021372/1‑148 (MQ=255)
                        aaGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATActct                                                                                                                          <  2:965637/149‑1 (MQ=255)
                         aGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATActct                                                                                                                          >  4:201632/1‑148 (MQ=255)
                          gAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctgt                                                                                                                        <  1:376314/149‑1 (MQ=255)
                               gtgtCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGa                                                                                                                    <  4:664222/148‑1 (MQ=255)
                                  tCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATta                                                                                                                 <  3:265551/148‑1 (MQ=255)
                                    ccTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATaa                                                                                                              <  2:1049723/149‑1 (MQ=255)
                                          aGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGa                                                                                                         <  4:49162/148‑1 (MQ=255)
                                           gACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAg                                                                                                       <  2:950204/149‑1 (MQ=255)
                                            aCGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAg                                                                                                       >  3:442319/1‑148 (MQ=255)
                                              gCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTTATTATAATCAGAAGGc                                                                                                     <  4:668419/148‑1 (MQ=255)
                                               cTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCa                                                                                                    <  3:792682/148‑1 (MQ=255)
                                                tGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctgtc                                                                                                                       >  2:864284/1‑128 (MQ=255)
                                                tGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAg                                                                                                   >  3:355052/1‑148 (MQ=255)
                                                    tCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTa                                                                                              >  1:20360/1‑149 (MQ=255)
                                                     cAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTa                                                                                              >  4:514943/1‑148 (MQ=255)
                                                       aaGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATcc                                                                                           <  1:111793/149‑1 (MQ=255)
                                                        aGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATcc                                                                                           <  3:201632/148‑1 (MQ=255)
                                                        aGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCt                                                                                          <  2:1018119/149‑1 (MQ=255)
                                                         ggTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCt                                                                                          >  3:750794/1‑148 (MQ=255)
                                                            cacaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGaa                                                                                      <  2:277008/149‑1 (MQ=255)
                                                            cacaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGaa                                                                                      <  2:187618/149‑1 (MQ=255)
                                                            cacaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGaa                                                                                      <  2:942061/149‑1 (MQ=255)
                                                            cacaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGa                                                                                       <  3:586385/148‑1 (MQ=255)
                                                             acaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAAc                                                                                     >  2:468836/1‑149 (MQ=255)
                                                             acaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAAc                                                                                     >  1:1043379/1‑149 (MQ=255)
                                                               aaTGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAg                                                                                   >  2:651032/1‑149 (MQ=255)
                                                               aaTGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAg                                                                                   >  2:1008793/1‑149 (MQ=255)
                                                                aTGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGa                                                                                  <  1:560635/149‑1 (MQ=255)
                                                                 tGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGa                                                                                  <  3:160497/148‑1 (MQ=255)
                                                                 tGGCTGGGAGGCTGGCGAGAAGCCATAGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGa                                                                                  >  3:882821/1‑148 (MQ=255)
                                                                       ggAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACt                                                                            >  4:265506/1‑148 (MQ=255)
                                                                        gAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACttt                                                                          <  1:761950/149‑1 (MQ=255)
                                                                          ggCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCa                                                                        >  2:855545/1‑149 (MQ=255)
                                                                              ggCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGt                                                                     <  4:1021372/148‑1 (MQ=255)
                                                                                  agaAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctgtctg                                                                                                                     <  3:958394/96‑1 (MQ=255)
                                                                                  agaAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctgtctg                                                                                                                     >  4:958394/1‑96 (MQ=255)
                                                                                      gCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCg                                                             <  4:750794/148‑1 (MQ=255)
                                                                                       ccATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCgg                                                            <  3:179607/148‑1 (MQ=255)
                                                                                             cGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTc                                                                         <  4:72562/129‑1 (MQ=255)
                                                                                             cGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTc                                                                         >  3:72562/1‑129 (MQ=255)
                                                                                               gACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATc                                                    <  3:265506/148‑1 (MQ=255)
                                                                                                  aGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAgac                                                <  2:711575/149‑1 (MQ=255)
                                                                                                    ggTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAgacg                                               <  3:598177/148‑1 (MQ=255)
                                                                                                     gtagtaATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAgacgac                                             <  2:310722/149‑1 (MQ=255)
                                                                                                     gtagtaATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAgacgac                                             <  2:431490/149‑1 (MQ=255)
                                                                                                      tagtaATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATc                                                    <  2:751183/141‑1 (MQ=255)
                                                                                                      tagtaATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATc                                                    >  1:751183/1‑141 (MQ=255)
                                                                                                         taATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATc                                         >  2:25242/1‑149 (MQ=255)
                                                                                                           aTGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCaa                                       >  1:224897/1‑149 (MQ=255)
                                                                                                               ccTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAg                                   <  2:914302/149‑1 (MQ=255)
                                                                                                                cTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCATTTTGATACCCTGTCTGATTATAATCAGAAGCCAGCCGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGc                                  >  1:696054/1‑149 (MQ=255)
                                                                                                                 tGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGc                                  >  3:461535/1‑148 (MQ=255)
                                                                                                                 tGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGc                                  >  4:125636/1‑148 (MQ=255)
                                                                                                                  gATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAga                                                 <  4:587817/132‑1 (MQ=255)
                                                                                                                    ttCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCa                                                   >  2:417040/1‑128 (MQ=255)
                                                                                                                    ttCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCa                                                   <  1:417040/128‑1 (MQ=255)
                                                                                                                    ttCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACt                               >  4:995346/1‑148 (MQ=255)
                                                                                                                         tctcTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCa                         >  1:667104/1‑149 (MQ=255)
                                                                                                                            ctTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATcc                      >  2:854996/1‑149 (MQ=255)
                                                                                                                            ctTTAAAAAGAAACCTCCGCATTGCGGAGGGTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATcc                      <  2:696054/149‑1 (MQ=255)
                                                                                                                             tttAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAgg                                                                      <  1:1014553/100‑1 (MQ=255)
                                                                                                                             tttAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAgg                                                                      >  2:1014553/1‑100 (MQ=255)
                                                                                                                             tttAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCGGATTTTAATCAGAAGGCAGACGTATCCTGAAACAGACCGACTTTCAGGTGGCTGGCGGTATAGATAAGCCGACCATCAACCAGCACTTCGCCATCCg                     >  2:327506/1‑149 (MQ=255)
                                                                                                                              ttAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCg                     >  3:454789/1‑148 (MQ=255)
                                                                                                                               tAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAACCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGc                    >  4:818635/1‑148 (MQ=255)
                                                                                                                                 aaaaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAg                 >  1:47014/1‑149 (MQ=255)
                                                                                                                                 aaaaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAg                 >  2:812562/1‑149 (MQ=255)
                                                                                                                                   aaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCATGc               <  2:140075/149‑1 (MQ=255)
                                                                                                                                   aaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGc               <  1:1008793/149‑1 (MQ=255)
                                                                                                                                   aaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATAATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTGTCATGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAgg                >  4:365698/1‑148 (MQ=255)
                                                                                                                                     gAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGccc             <  2:142829/149‑1 (MQ=255)
                                                                                                                                     gAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGcc              <  4:454789/148‑1 (MQ=255)
                                                                                                                                       aaCCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAgg                <  4:492583/144‑1 (MQ=255)
                                                                                                                                       aaCCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAgg                >  3:492583/1‑144 (MQ=255)
                                                                                                                                        aCCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATa          >  1:815256/1‑149 (MQ=255)
                                                                                                                                           tCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATc       >  2:944127/1‑149 (MQ=255)
                                                                                                                                           tCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATc       >  1:91390/1‑149 (MQ=255)
                                                                                                                                           tCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATc       >  1:232625/1‑149 (MQ=255)
                                                                                                                                            ccGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATc       >  3:189350/1‑148 (MQ=255)
                                                                                                                                            ccGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATc       >  3:862925/1‑148 (MQ=255)
                                                                                                                                            ccGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATc       >  3:826221/1‑148 (MQ=255)
                                                                                                                                                aTTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAgacg  <  1:468836/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACG  >  NC_000913/1015767‑1016059

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