Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,962,562 T→C 7.6% N68S (AAC→AGC) ‡ ppdC ← putative prepilin peptidase‑dependent protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,962,5620TC 7.6% 292.3 / 12.4 107N68S (AAC→AGC) ‡ppdCputative prepilin peptidase‑dependent protein
Reads supporting (aligned to +/- strand):  ref base T (50/48);  new base C (3/5);  total (53/54)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

Only 100 of 108 total aligned reads displayed.

GGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCGC  >  NC_000913/2962413‑2962710
                                                                                                                                                     |                                                                                                                                                    
ggTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGt                                                                                                                                                      >  2:83292/1‑150 (MQ=255)
  tAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTg                                                                                                                                                    <  2:856168/150‑1 (MQ=255)
    gACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGc                                                                                                                                                  <  1:1383797/150‑1 (MQ=255)
      cccTTAACATAGCGGCTCCTGACTACTGACGAATCGGGCAATGAAGGCAGGTCATATCGCCTTCTCTGCCCCCCGGTGAAACTAGCGTAAAGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGGAg                                                                                                                                                >  1:891938/1‑150 (MQ=255)
       ccTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACCCATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAgg                                                                                                                                               >  1:607686/1‑150 (MQ=255)
       ccTTAACATAGCGGCTCCTGACTACTGACGAATCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGATGACACATCCCGCCTGAGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAgg                                                                                                                                               >  1:1195080/1‑150 (MQ=255)
        cTTAACATAGCGGCTCCTGACTACTGACTATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGt                                                                                                                                              <  1:1111320/150‑1 (MQ=255)
         ttAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtc                                                                                                                                             <  1:662933/150‑2 (MQ=255)
         ttAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtc                                                                                                                                             <  1:1591864/150‑2 (MQ=255)
          tAACTTAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTgg                                                                                                                                            <  1:1328353/150‑1 (MQ=255)
          tAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCAACTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTCCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTgg                                                                                                                                            >  1:360875/1‑150 (MQ=255)
          tAACATACCTGCTCCTGCCTACTGCCGATTCGTGCAATGCAGGCGGGTCATCGCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTgg                                                                                                                                            <  2:106259/150‑1 (MQ=255)
           accATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGGCTGCAATCGGTTGACCTGCCAGTTGGCAGGTGGc                                                                                                                                           <  1:1449364/148‑1 (MQ=255)
            aCATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCTGTTGACCTGCCAACTGGCAGGtcaac                                                                                                                                          <  1:1323113/150‑5 (MQ=255)
            aCATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCg                                                                                                                                          >  2:1211306/1‑150 (MQ=255)
              aTAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGaa                                                                                                                                        >  1:806492/1‑150 (MQ=255)
                aGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAAt                                                                                                                                      <  2:1051343/150‑1 (MQ=255)
                  cGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATcg                                                                                                                                    >  2:86866/1‑150 (MQ=255)
                    gCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATcgcg                                                                                                                                  >  1:626466/1‑150 (MQ=255)
                    gCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCACTTGGCAGGTGGCGAAATcgcg                                                                                                                                  >  1:553132/1‑150 (MQ=255)
                    ccTCCTGCCTACTGACGATTCGGGCACTGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCTTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTGGACCTGCCAGTTGGCAGGTGGCGAAATcgcg                                                                                                                                  <  2:144325/149‑1 (MQ=255)
                      tCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATcgcgcg                                                                                                                                <  1:807646/150‑1 (MQ=255)
                         tGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAg                                                                                                                             >  1:498537/1‑150 (MQ=255)
                            ctactGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGTAATCGCGCGCAGGTg                                                                                                                          >  1:336106/1‑150 (MQ=255)
                             tactGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCGTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACCCATCCCGCCTGCGAGGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGc                                                                                                                         >  1:598873/1‑150 (MQ=255)
                              actGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCg                                                                                                                        >  1:1491948/1‑150 (MQ=255)
                               ctGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGt                                                                                                                       >  1:1097661/1‑150 (MQ=255)
                                 gACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGTGttt                                                                                                                     >  1:1074239/1‑150 (MQ=255)
                                 gACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGttt                                                                                                                     <  1:1416292/150‑1 (MQ=255)
                                   cGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTTACCTGCCAGTTGGAAGGTGGCGAAATCGCGCGCAGTTGCGTTtgc                                                                                                                   >  2:1048795/1‑150 (MQ=255)
                                      ttCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTtgctgc                                                                                                                <  2:1368801/150‑1 (MQ=255)
                                       tCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgcc                                                                                                               <  2:1648770/150‑1 (MQ=255)
                                        cGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgcca                                                                                                              <  1:91890/150‑1 (MQ=255)
                                         gggCAATGCAGGCGGGTCATCTCGGCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgccag                                                                                                             <  2:891938/150‑1 (MQ=255)
                                         gggCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgccag                                                                                                             >  2:1644001/1‑150 (MQ=255)
                                          ggCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgccagc                                                                                                            >  1:887551/1‑150 (MQ=255)
                                            cAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgccagcca                                                                                                          >  1:443470/1‑150 (MQ=255)
                                            cAATGCAGGCGGGTCATCTCGCCTTCTCTGACCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgccagcca                                                                                                          >  2:791914/1‑150 (MQ=255)
                                             aaTGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCAt                                                                                                         <  2:594714/150‑1 (MQ=255)
                                              aTGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATg                                                                                                        >  2:1098574/1‑150 (MQ=255)
                                              aTGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATg                                                                                                        <  2:443470/150‑1 (MQ=255)
                                                 cAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATgccg                                                                                                     <  2:1598134/150‑1 (MQ=255)
                                                           aTCTCGCCTTCTATTCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTTGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGt                                                                                           >  2:748600/1‑150 (MQ=255)
                                                               cgacTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGt                                                                                       <  1:1577583/147‑1 (MQ=255)
                                                                        tGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTc                                                                              <  2:357210/150‑1 (MQ=255)
                                                                        tGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTc                                                                              <  2:626466/150‑1 (MQ=255)
                                                                         gCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTct                                                                             >  2:149022/1‑150 (MQ=255)
                                                                         gCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTct                                                                             <  2:1097661/150‑1 (MQ=255)
                                                                          cccccGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTctg                                                                            <  2:498537/150‑1 (MQ=255)
                                                                            cccGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTctgct                                                                          <  1:524498/150‑1 (MQ=255)
                                                                            cccGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTctgct                                                                          <  2:1257849/150‑1 (MQ=255)
                                                                              cGGGTGACACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGCTGACCTGCCAGATGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCAGCTgg                                                                        <  2:407613/150‑1 (MQ=255)
                                                                              cGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTgg                                                                        >  2:524110/1‑150 (MQ=255)
                                                                               gggTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGc                                                                       >  2:1600165/1‑150 (MQ=255)
                                                                                 gTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCaa                                                                     >  1:1191426/1‑150 (MQ=255)
                                                                                  tGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGACGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAg                                                                    >  1:352283/1‑150 (MQ=255)
                                                                                    aaccTAGCGGACCGCTGATGCTGACCCATCCCGCCTGCGATGTCGGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGAc                                                                  <  1:350778/147‑1 (MQ=255)
                                                                                       cTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGt                                                               >  2:1229268/1‑150 (MQ=255)
                                                                                        tAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGtt                                                              <  1:1572753/150‑1 (MQ=255)
                                                                                        tAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGtt                                                              <  2:761371/150‑1 (MQ=255)
                                                                                          gCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCa                                                            >  2:686861/1‑150 (MQ=255)
                                                                                           ggTAACGCTGATGCTGTCACATTCCGCCTGCGATGTCTGCATTCGGTTTACCTGCCATTTGGCAGGTGGCGAATTCGCGCGCAGGTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCAt                                                           <  2:1081879/149‑1 (MQ=255)
                                                                                           cGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCAt                                                           <  1:149022/150‑1 (MQ=255)
                                                                                                cGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATg                                                      <  1:243025/150‑1 (MQ=255)
                                                                                                cGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATg                                                      >  1:254437/1‑150 (MQ=255)
                                                                                                cGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATg                                                      <  2:755182/150‑1 (MQ=255)
                                                                                                 gCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGt                                                     >  2:1090869/1‑150 (MQ=255)
                                                                                                  cTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGtt                                                    <  1:686861/150‑1 (MQ=255)
                                                                                                      tGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCt                                                <  2:468568/150‑1 (MQ=255)
                                                                                                      tGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCt                                                >  1:107177/1‑150 (MQ=255)
                                                                                                      tGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCt                                                <  2:887551/150‑1 (MQ=255)
                                                                                                      tGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCt                                                >  1:703085/1‑150 (MQ=255)
                                                                                                       gCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTg                                               <  1:579000/150‑1 (MQ=255)
                                                                                                        cTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCTAGACTGTTCATTAATGTTCGCTgg                                              >  2:174711/1‑150 (MQ=255)
                                                                                                           acacATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGaaa                                           <  2:360485/150‑1 (MQ=255)
                                                                                                             acaTCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAAcc                                         <  2:1491948/150‑1 (MQ=255)
                                                                                                                     ccTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACg                                 <  2:39826/150‑1 (MQ=255)
                                                                                                                      cTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTAGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGc                                <  2:202285/150‑1 (MQ=255)
                                                                                                                            tGTCTGCATTCGGTTGACCTGCCAGTTTGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGAc                          >  2:1371032/1‑150 (MQ=255)
                                                                                                                               cTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAat                       >  1:971651/1‑150 (MQ=255)
                                                                                                                                 gCATTCGGTTGAGCTGCCAGTTGGCAGGTGGCGTAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGTGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGTAACCCGATAACGCAGTGACAatca                     <  2:360875/150‑1 (MQ=255)
                                                                                                                                 gCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAatca                     <  1:1503886/150‑1 (MQ=255)
                                                                                                                                 gCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAatca                     <  2:1191426/150‑1 (MQ=255)
                                                                                                                                       ggTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCa               >  2:1645644/1‑150 (MQ=255)
                                                                                                                                        gTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCAt              >  2:443994/1‑150 (MQ=255)
                                                                                                                                         ttGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATc             <  1:1211306/150‑1 (MQ=255)
                                                                                                                                          tGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCa            >  1:1117549/1‑150 (MQ=255)
                                                                                                                                          ccaCCTGCCAACTGGCAGGTGGCGAAATCGCGCGCAGTTNCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCa            >  2:1173532/3‑150 (MQ=255)
                                                                                                                                          ccaCCTGCCAACTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCa            >  1:781043/3‑150 (MQ=255)
                                                                                                                                           gACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCaa           <  2:1367900/150‑1 (MQ=255)
                                                                                                                                             ccTGCCAGTTGGCAGGTGGCGAAATCGCGCTCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCCCCATCAaca         >  2:1232063/1‑150 (MQ=255)
                                                                                                                                             ccTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAaca         >  2:1561539/1‑150 (MQ=255)
                                                                                                                                             ccTGCCAGTTGGCAGGTGGAGAAATCGCGCGAAGTTGCGTTTGCTGCAAGCCATGCCGAAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACATGAATCACAATCAacac        >  1:729950/1‑150 (MQ=255)
                                                                                                                                              cTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAacac        >  2:419940/1‑150 (MQ=255)
                                                                                                                                               tGCCAACTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACAcc       >  2:1674569/1‑150 (MQ=255)
                                                                                                                                                gCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTTCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGATTCTGCTGGCAAGAATGTTCATTAATGTTCGCTGGAAAACCGATTACGCAGTGACAATCATCACCATCAacacaa      >  2:509209/1‑148 (MQ=255)
                                                                                                                                                   aGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGCCTGTTCATTAATGTTCGCTGGAACCCCGATAACGCAGTGCCAATCATCACCTTCAACCCCATCg   >  1:845421/1‑150 (MQ=255)
                                                                                                                                                   aGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATTATCACCATCAACACCATCg   >  1:792916/1‑150 (MQ=255)
                                                                                                                                                   aGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCg   <  2:379616/150‑1 (MQ=255)
                                                                                                                                                    gTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCGc  >  1:1218448/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCATCGC  >  NC_000913/2962413‑2962710

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