Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,386,129 A→G 5.2% 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.2% 373.9 / 6.4 114M1081T (ATG→ACG) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base A (50/58);  new base G (4/2);  total (54/60)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 115 total aligned reads displayed.

GGCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATC  >  NC_000913/4385984‑4386277
                                                                                                                                                 |                                                                                                                                                    
ggCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCt                                                                                                                                                   <  4:247172/149‑1 (MQ=255)
ggCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCt                                                                                                                                                   <  3:915763/149‑1 (MQ=255)
 gCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTg                                                                                                                                                  >  2:974211/1‑149 (MQ=255)
 gCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTg                                                                                                                                                  >  2:60761/1‑149 (MQ=255)
 gCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCAGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTg                                                                                                                                                  >  4:704853/1‑149 (MQ=255)
     ttttttCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGATGCGTCTGGaaa                                                                                                                                              <  3:485344/149‑2 (MQ=255)
     ttttttCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGATGCGTCTGGaaa                                                                                                                                              <  3:547460/149‑2 (MQ=255)
        tttCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGc                                                                                                                                           <  3:182734/149‑1 (MQ=255)
              cgcTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAAt                                                                                                                                     <  1:449737/149‑1 (MQ=255)
                cTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATgg                                                                                                                                   <  2:355828/149‑1 (MQ=255)
                  aTCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCa                                                                                                                                 <  2:736228/149‑1 (MQ=255)
                      gTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATc                                                                                                                             <  1:602330/149‑1 (MQ=255)
                      cttttGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATc                                                                                                                             <  1:974211/148‑1 (MQ=255)
                       ttttGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCg                                                                                                                            >  1:1063292/1‑149 (MQ=255)
                       ttttGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTGTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCg                                                                                                                            >  2:908752/1‑149 (MQ=255)
                        tttgtttgTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGa                                                                                                                           <  1:302619/149‑1 (MQ=255)
                           gtttgtGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATAc                                                                                                                        >  4:1009105/1‑149 (MQ=255)
                           gtttgtGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATAc                                                                                                                        >  4:273157/1‑149 (MQ=255)
                               gtgAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTg                                                                                                                    <  3:974717/149‑1 (MQ=255)
                                tgAGCCGGATTGGTTCATCCGGTACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTgg                                                                                                                   <  3:663367/149‑1 (MQ=255)
                                tgAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTgg                                                                                                                   <  4:1072798/149‑1 (MQ=255)
                                    ccGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCAt                                                                                                               <  4:393434/149‑1 (MQ=255)
                                        aTTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGTGGGAATGGCATATCGATACCGTGGGCATGGaa                                                                                                           <  4:1142/149‑1 (MQ=255)
                                        aTTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGaa                                                                                                           <  1:60761/149‑1 (MQ=255)
                                        aTTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGaa                                                                                                           <  3:273157/149‑1 (MQ=255)
                                        aTTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGaa                                                                                                           <  2:230223/149‑1 (MQ=255)
                                            gTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagcca                                                                                                       >  2:205293/1‑149 (MQ=255)
                                            gTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagcca                                                                                                       >  1:225067/1‑149 (MQ=255)
                                             ttCATCCGGCCCACCAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATAc                                                                                                                        <  2:216634/131‑1 (MQ=255)
                                             ttCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagccag                                                                                                      >  3:744014/1‑149 (MQ=255)
                                             ttCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagccag                                                                                                      >  2:949164/1‑149 (MQ=255)
                                             ttCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGGTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATAc                                                                                                                        >  1:216634/1‑131 (MQ=255)
                                               cATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagccagcc                                                                                                    <  3:959077/149‑1 (MQ=255)
                                                 tCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGccagccag                                                                                                  >  3:29266/1‑149 (MQ=255)
                                                 tCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGccagccag                                                                                                  >  2:675718/1‑149 (MQ=255)
                                                    ggCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAAt                                                                                               <  2:817655/149‑1 (MQ=255)
                                                       tcacaAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAg                                                                                            >  2:39824/2‑149 (MQ=255)
                                                       acacaAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAg                                                                                            <  3:169780/149‑1 (MQ=255)
                                                         acaaacttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCt                                                                                          >  4:249463/1‑149 (MQ=255)
                                                            aacttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtggt                                                                                       <  2:829193/149‑1 (MQ=255)
                                                            aacttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtggt                                                                                       <  1:326314/149‑1 (MQ=255)
                                                               ttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGaa                                                                                                           <  1:46417/126‑1 (MQ=255)
                                                               ttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGaa                                                                                                           >  2:46417/1‑126 (MQ=255)
                                                                tacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGAt                                                                                   >  1:586294/1‑149 (MQ=255)
                                                                  caaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGAtct                                                                                 >  1:1095915/1‑149 (MQ=255)
                                                                      cttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTCCCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATg                                                                             <  4:1037675/149‑1 (MQ=255)
                                                                      cttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATg                                                                             <  3:327700/149‑1 (MQ=255)
                                                                        tCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGc                                                                           <  1:949102/149‑1 (MQ=255)
                                                                        tCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGc                                                                           <  1:1094527/149‑1 (MQ=255)
                                                                         cccGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCg                                                                          >  4:547211/1‑149 (MQ=255)
                                                                              cTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCg                                                                     >  3:795112/1‑149 (MQ=255)
                                                                               tGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCgg                                                                    <  2:80604/149‑1 (MQ=255)
                                                                                     ccTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGTATACg                                                              <  1:908752/149‑1 (MQ=255)
                                                                                      cTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACgg                                                             <  4:1048411/149‑1 (MQ=255)
                                                                                      cTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACgg                                                             >  3:511713/1‑149 (MQ=255)
                                                                                        ttGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTg                                                           <  1:160460/149‑1 (MQ=255)
                                                                                        ttGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTg                                                           <  4:292871/149‑1 (MQ=255)
                                                                                             cGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagcca                                                                                                       >  1:129360/1‑100 (MQ=255)
                                                                                             cGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagcca                                                                                                       <  2:129360/100‑1 (MQ=255)
                                                                                              gCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCAt                                                     <  3:249463/149‑1 (MQ=255)
                                                                                                     tCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAgcg                                              <  3:618962/149‑1 (MQ=255)
                                                                                                      cAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtggt                                                                                       <  4:8048/107‑1 (MQ=255)
                                                                                                      cAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtggt                                                                                       >  3:8048/1‑107 (MQ=255)
                                                                                                          gTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCTGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTa                                         >  2:18824/1‑149 (MQ=255)
                                                                                                           ttCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTaa                                        >  2:509212/1‑149 (MQ=255)
                                                                                                            tctcCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTaaa                                       <  2:586294/149‑1 (MQ=255)
                                                                                                              tcCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATa                                     <  4:719424/149‑1 (MQ=255)
                                                                                                               ccccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATAc                                    <  3:717467/149‑1 (MQ=255)
                                                                                                                cccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAgccagcca                                                                                                   <  2:154985/85‑1 (MQ=255)
                                                                                                                cccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAgccagcca                                                                                                   >  1:154985/1‑85 (MQ=255)
                                                                                                                cccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACg                                   >  1:767815/1‑149 (MQ=255)
                                                                                                                 ccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGTGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGc                                  <  3:704853/149‑1 (MQ=255)
                                                                                                                  cGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCa                                 <  1:675718/149‑1 (MQ=255)
                                                                                                                  cGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCa                                 >  4:826894/1‑149 (MQ=255)
                                                                                                                   gtttgtttACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAGATACGCAg                                >  4:620926/1‑149 (MQ=255)
                                                                                                                   gtttgtttACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAg                                <  4:795112/149‑1 (MQ=255)
                                                                                                                   gtttgtttACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAg                                <  4:724219/149‑1 (MQ=255)
                                                                                                                        tttACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGa                           >  3:466495/1‑149 (MQ=255)
                                                                                                                          tACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGaaa                         >  4:619314/1‑149 (MQ=255)
                                                                                                                              gTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCt                     <  1:1017013/149‑1 (MQ=255)
                                                                                                                              gTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCt                     >  3:309177/1‑149 (MQ=255)
                                                                                                                              gTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGAGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCt                     <  1:193119/149‑1 (MQ=255)
                                                                                                                                 gAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGaa                  >  4:811798/1‑149 (MQ=255)
                                                                                                                                  aGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAAt                 <  1:950944/149‑1 (MQ=255)
                                                                                                                                   ggCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATc                >  4:116265/1‑149 (MQ=255)
                                                                                                                                   gcccttCCAGATGCGTCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATc                >  1:852349/5‑149 (MQ=255)
                                                                                                                                   gcccttCCAGATGCGTCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCAGGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATc                >  4:57541/5‑149 (MQ=255)
                                                                                                                                    gCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAgcggc                                            <  3:781786/120‑1 (MQ=255)
                                                                                                                                    gCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAgcggc                                            >  4:781786/1‑120 (MQ=255)
                                                                                                                                     ccttCCAGATGCGTCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATccc              >  4:559756/3‑149 (MQ=255)
                                                                                                                                     cTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCa                                 >  4:326944/1‑130 (MQ=255)
                                                                                                                                     cTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCa                                 <  3:326944/130‑1 (MQ=255)
                                                                                                                                      cttCCAGATGCGTCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGTAGCTCGAAAATCTGAATcccc             >  4:782850/2‑149 (MQ=255)
                                                                                                                                       ttCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCt            <  2:928828/149‑1 (MQ=255)
                                                                                                                                       ttCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCt            >  3:661178/1‑149 (MQ=255)
                                                                                                                                        tCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTg           <  2:277886/149‑1 (MQ=255)
                                                                                                                                           aGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTg        >  3:378111/1‑149 (MQ=255)
                                                                                                                                           aGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTg        >  4:112643/1‑149 (MQ=255)
                                                                                                                                           aGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTg        >  3:243622/1‑149 (MQ=255)
                                                                                                                                            gACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGc       >  2:666296/1‑149 (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 < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: