Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,386,129 A→G 5.7% M1081T (ATG→ACG)  mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,386,1290AG 5.7% 386.5 / 10.4 122M1081T (ATG→ACG) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base A (51/64);  new base G (4/3);  total (55/67)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-02

Only 100 of 125 total aligned reads displayed.

GGCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATC  >  NC_000913/4385984‑4386277
                                                                                                                                                 |                                                                                                                                                    
ggCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCt                                                                                                                                                   <  3:928533/149‑1 (MQ=255)
 gCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTg                                                                                                                                                  >  3:325386/1‑149 (MQ=255)
    gTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGaa                                                                                                                                               <  4:453071/149‑1 (MQ=255)
    gTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGaa                                                                                                                                               <  4:71497/149‑1 (MQ=255)
     ttttttCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAg                                                                                                                                              >  4:884226/1‑149 (MQ=255)
        tttCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGATGCGTCTGGaaagcc                                                                                                                                           <  1:720744/149‑5 (MQ=255)
        tttCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGc                                                                                                                                           <  3:111136/149‑1 (MQ=255)
          tCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCgg                                                                                                                                         <  2:860698/149‑1 (MQ=255)
            tGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGATGCGTCTGGaaa                                                                                                                                              >  2:835436/1‑141 (MQ=255)
            tGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGATGCGTCTGGaaa                                                                                                                                              <  1:835436/142‑2 (MQ=255)
            tGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGa                                                                                                                                       >  4:2462/1‑149 (MQ=255)
             gcgcTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGaa                                                                                                                                      <  4:603250/149‑1 (MQ=255)
                cTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATgg                                                                                                                                   <  4:62374/149‑1 (MQ=255)
                cTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATgg                                                                                                                                   >  3:578299/1‑149 (MQ=255)
                  aTCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATgg                                                                                                                                   <  1:204665/147‑1 (MQ=255)
                   tCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCat                                                                                                                                >  2:41405/1‑149 (MQ=255)
                      gTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATc                                                                                                                             >  3:437243/1‑149 (MQ=255)
                       ttttGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCg                                                                                                                            >  1:60065/1‑149 (MQ=255)
                       ttttGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCg                                                                                                                            >  1:274518/1‑149 (MQ=255)
                       ttttGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCg                                                                                                                            >  3:231067/1‑149 (MQ=255)
                         ttgtttgtGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGAt                                                                                                                          <  2:939298/149‑1 (MQ=255)
                               gtgAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTg                                                                                                                    <  3:108179/149‑1 (MQ=255)
                                     cGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATg                                                                                                              <  4:744020/149‑1 (MQ=255)
                                      ggATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATgg                                                                                                             <  2:288380/149‑1 (MQ=255)
                                        aTTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGaa                                                                                                           <  3:738173/149‑1 (MQ=255)
                                        aTTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGaa                                                                                                           >  2:174204/1‑149 (MQ=255)
                                         ttGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAag                                                                                                          >  1:761012/1‑149 (MQ=255)
                                            gTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagcca                                                                                                       >  2:602230/1‑149 (MQ=255)
                                            gTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGCCATATCGATACCGTGGGCATGGAagcca                                                                                                       <  3:71002/149‑1 (MQ=255)
                                             ttCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagccag                                                                                                      >  1:856233/1‑149 (MQ=255)
                                              tCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagccagc                                                                                                     <  2:60065/149‑1 (MQ=255)
                                                aTCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAgg                                                                                                                                             >  4:939808/1‑107 (MQ=255)
                                                aTCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAgg                                                                                                                                             <  3:939808/107‑1 (MQ=255)
                                                aTCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGaa                                                                                                           <  1:685483/141‑1 (MQ=255)
                                                aTCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGaa                                                                                                           >  2:685483/1‑141 (MQ=255)
                                                aTCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAgccagcca                                                                                                   >  1:94637/1‑149 (MQ=255)
                                                aTCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAgccagcca                                                                                                   >  1:11157/1‑149 (MQ=255)
                                                 tCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGATTGTTTACCGTTGAGGCTTTCCAGACGTATCTGGAAGGGCGTGAATGGCATATCGATACCGTGGGCATGGAAGccagccag                                                                                                  >  2:1027683/1‑149 (MQ=255)
                                                   cGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGaa                                                                                                <  1:502823/149‑1 (MQ=255)
                                                     gCACACAGACTTACAAACGTCCCGCCGGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATc                                                                                              <  1:196643/149‑1 (MQ=255)
                                                     gCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGaa                                                                                                           >  4:288665/1‑136 (MQ=255)
                                                     gCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGaa                                                                                                           <  3:288665/136‑1 (MQ=255)
                                                      cacacaAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCa                                                                                             <  2:825611/149‑1 (MQ=255)
                                                       acacacacTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAg                                                                                            >  2:186936/1‑149 (MQ=255)
                                                          caaacttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtg                                                                                         <  2:94637/149‑1 (MQ=255)
                                                           aaacttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtgg                                                                                        >  3:534824/1‑149 (MQ=255)
                                                             acttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtggtg                                                                                      <  4:519764/149‑1 (MQ=255)
                                                                   aaacttCCCGCCTGGCGACCTTTGTCGGCAGACGTCAGCGTTCTCCCCGTTTGTTGACCGTTGCGGCTTTCCAGACGCATCTAGAAGGTCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGAtctc                                                                                <  1:363307/149‑1 (MQ=255)
                                                                   aaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCata                                                                                                                               >  2:548933/1‑102 (MQ=255)
                                                                   aaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCata                                                                                                                               <  1:548933/102‑1 (MQ=255)
                                                                   aaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGAtctc                                                                                <  4:37120/149‑1 (MQ=255)
                                                                   aaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGAtctc                                                                                >  1:1043810/1‑149 (MQ=255)
                                                                   aaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGAtctc                                                                                <  2:523708/149‑1 (MQ=255)
                                                                   aaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCA‑‑TGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCAt                                                                              >  3:582442/1‑149 (MQ=255)
                                                                    aacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCa                                                                               <  1:115514/149‑1 (MQ=255)
                                                                    aacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCa                                                                               <  3:456688/149‑1 (MQ=255)
                                                                    aacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCa                                                                               <  3:614787/149‑1 (MQ=255)
                                                                    aacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCa                                                                               <  2:199468/149‑1 (MQ=255)
                                                                      cttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATg                                                                             <  1:461180/149‑1 (MQ=255)
                                                                       ttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATgg                                                                            >  3:928500/1‑149 (MQ=255)
                                                                             ccTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGc                                                                      >  1:1036557/1‑149 (MQ=255)
                                                                             ccTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGc                                                                      <  2:334497/149‑1 (MQ=255)
                                                                                   gacCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATa                                                                >  2:44464/1‑149 (MQ=255)
                                                                                   gacCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATa                                                                >  3:154630/1‑149 (MQ=255)
                                                                                      cTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATg                                                                             <  1:723873/133‑1 (MQ=255)
                                                                                      cTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATg                                                                             >  2:723873/1‑133 (MQ=255)
                                                                                      cTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACgg                                                             <  1:710679/149‑1 (MQ=255)
                                                                                         tGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGa                                                          <  3:537147/149‑1 (MQ=255)
                                                                                          gCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGAc                                                         <  3:916595/149‑1 (MQ=255)
                                                                                               cAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGt                                                            <  1:896800/141‑1 (MQ=255)
                                                                                               cAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGt                                                            >  2:896800/1‑141 (MQ=255)
                                                                                                   cGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAg                                                <  2:838886/149‑1 (MQ=255)
                                                                                                       aGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGaa                                                                                                <  3:588417/97‑1 (MQ=255)
                                                                                                       aGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGaa                                                                                                >  4:588417/1‑97 (MQ=255)
                                                                                                              tcCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGTAATGGCATATCGATATCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATCCGGTGACCCATCTCAGCCGCGTAAATa                                     >  4:377828/1‑149 (MQ=255)
                                                                                                              tcCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATa                                     >  4:858098/1‑149 (MQ=255)
                                                                                                               ccccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATAc                                    <  4:44149/149‑1 (MQ=255)
                                                                                                                cccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACg                                   >  2:449715/1‑149 (MQ=255)
                                                                                                                  cGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCa                                 <  2:856233/149‑1 (MQ=255)
                                                                                                                  cGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCa                                 <  4:437243/149‑1 (MQ=255)
                                                                                                                   gtttgtttACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATAc                                    <  3:89791/145‑1 (MQ=255)
                                                                                                                   gtttgtttACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATAc                                    >  4:89791/1‑145 (MQ=255)
                                                                                                                       gtttACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCg                            >  2:382694/1‑149 (MQ=255)
                                                                                                                           aCCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGaaaa                        <  3:84764/149‑1 (MQ=255)
                                                                                                                              gTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCt                     <  4:689518/149‑1 (MQ=255)
                                                                                                                              gTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCt                     >  2:274684/1‑149 (MQ=255)
                                                                                                                               ttGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATc                                                                                                                             <  1:630236/44‑1 (MQ=255)
                                                                                                                               ttGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATc                                                                                                                             >  2:630236/1‑44 (MQ=255)
                                                                                                                                  aGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAAt                 <  1:862003/149‑1 (MQ=255)
                                                                                                                                     cTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATccc              <  2:11157/149‑1 (MQ=255)
                                                                                                                                      tttCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCATAATCAGCTGGTGGATCTCATGGCGTAGCCGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATcccc             <  3:310992/149‑1 (MQ=255)
                                                                                                                                      cttCCAGATGCGTCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATcccc             >  2:1044599/2‑149 (MQ=255)
                                                                                                                                       ttCCAGATGCGTCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAgccagcca                                                                                                   >  4:281056/1‑62 (MQ=255)
                                                                                                                                         ccAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGt          <  1:663813/149‑1 (MQ=255)
                                                                                                                                          cAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGtt         >  1:701019/1‑149 (MQ=255)
                                                                                                                                               gCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGa    <  4:231067/149‑1 (MQ=255)
                                                                                                                                                cATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGAt   >  2:117612/1‑149 (MQ=255)
                                                                                                                                                cATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGAt   <  2:252349/149‑1 (MQ=255)
                                                                                                                                                 aTCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCAt                                                                              <  3:322474/73‑1 (MQ=255)
                                                                                                                                                 aTCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATc  <  2:902452/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATC  >  NC_000913/4385984‑4386277

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