Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,962,561 G→A 8.2% N68N (AAC→AAT) ‡ 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,5610GA 8.2% 366.4 / 14.6 122N68N (AAC→AAT) ‡ppdCputative prepilin peptidase‑dependent protein
Reads supporting (aligned to +/- strand):  ref base G (51/61);  new base A (3/7);  total (54/68)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

Only 100 of 122 total aligned reads displayed.

TAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCA  >  NC_000913/2962415‑2962706
                                                                                                                                                  |                                                                                                                                                 
tAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTg                                                                                                                                                <  1:695742/150‑1 (MQ=255)
tAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGGCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTg                                                                                                                                                <  2:1350513/150‑1 (MQ=255)
  gACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGc                                                                                                                                              >  1:1456689/1‑150 (MQ=255)
  gACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGc                                                                                                                                              >  1:1390555/1‑150 (MQ=255)
     ccTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAgg                                                                                                                                           <  2:924327/150‑1 (MQ=255)
     ccTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAgg                                                                                                                                           <  2:1334739/150‑1 (MQ=255)
        tAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtca                                                                                                                                        <  1:1629903/150‑3 (MQ=255)
         aaCATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGc                                                                                                                                       <  2:1650453/150‑1 (MQ=255)
          aCATATCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCg                                                                                                                                      <  2:1342486/150‑1 (MQ=255)
          aCATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtcaac                                                                                                                                      <  1:223639/150‑5 (MQ=255)
           cctAGCGGCACCTGACTCCTGACGATTCCGGCAAGGCAGGCGGGGCATCTCGCCTTCGCTGCCCCCGGGAGAACCTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtcaacc                                                                                                                                     <  2:1009129/148‑6 (MQ=255)
                cGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATcg                                                                                                                                >  1:357353/1‑150 (MQ=255)
                  gCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGCCATCTCGCCTTCTCTTCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATcgcg                                                                                                                              <  1:286688/150‑1 (MQ=255)
                   cTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGGCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCTTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATcgcgc                                                                                                                             <  2:1036538/150‑1 (MQ=255)
                          ctactGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAATTCGCGCGCAGTTg                                                                                                                      >  2:1287221/1‑150 (MQ=255)
                          ctactGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTg                                                                                                                      >  2:492652/1‑150 (MQ=255)
                            actGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCg                                                                                                                    >  2:1245554/1‑150 (MQ=255)
                            actGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCg                                                                                                                    >  2:776023/1‑150 (MQ=255)
                            actGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCg                                                                                                                    <  2:878258/150‑1 (MQ=255)
                              tGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCACTGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGtt                                                                                                                  <  1:1575737/150‑1 (MQ=255)
                               gACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGttt                                                                                                                 >  2:1393814/1‑150 (MQ=255)
                                 cGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTtgc                                                                                                               <  1:131063/150‑1 (MQ=255)
                                  gATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTTGCAGGTGGCGAAATCGCGCGCAGTTGCGTTtgct                                                                                                              >  2:361999/1‑150 (MQ=255)
                                  gATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTtgct                                                                                                              <  2:37011/150‑1 (MQ=255)
                                    tgcGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTtgctgc                                                                                                            <  2:909337/148‑1 (MQ=255)
                                      cGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTCAAACTAGCGTAACGCTGCTGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgcca                                                                                                          >  2:77760/1‑150 (MQ=255)
                                       gggCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgccag                                                                                                         <  1:1359703/150‑1 (MQ=255)
                                        ggCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGCCACATCCCGCCTGCGATGTCTGCATTCGGTTGCCCTGCCAGTTGGCTGGTGGCGAAGTCGCGCGCTGTTGCGTTTGCTgccagc                                                                                                        >  2:82441/1‑150 (MQ=255)
                                        ggCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgccagc                                                                                                        <  2:1663783/150‑1 (MQ=255)
                                         gCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgccagcc                                                                                                       >  1:1098136/1‑150 (MQ=255)
                                         gCAATGCAGGCGGGTCATCACGCCTTCTGTGCCCCCGGGTGAAACTAGAGTAACGCTGATGCTGACACATCACGCCTGCGAAGTATGCATTCGGTTGACCTGCCAGTAGGCAGGTGGCGAAATACCGAGCAGTTGAGTTTGATgccagcc                                                                                                       >  1:121697/1‑150 (MQ=255)
                                            aTGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTCGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATg                                                                                                    <  1:1074515/150‑1 (MQ=255)
                                               cAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATgccg                                                                                                 >  1:749906/1‑150 (MQ=255)
                                                      gTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGc                                                                                          <  2:238675/150‑1 (MQ=255)
                                                      gTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGc                                                                                          <  1:1023000/150‑1 (MQ=255)
                                                      gTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCCGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGc                                                                                          >  2:141154/1‑150 (MQ=255)
                                                          tctcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGtt                                                                                      <  2:1648039/150‑1 (MQ=255)
                                                          tctcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGtt                                                                                      <  1:1685021/150‑1 (MQ=255)
                                                           ctcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTg                                                                                     >  2:964181/1‑150 (MQ=255)
                                                            tcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGCCACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTgg                                                                                    <  1:387092/150‑1 (MQ=255)
                                                            tcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTgg                                                                                    >  2:1496794/1‑150 (MQ=255)
                                                             cGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTggtggt                                                                                   >  1:1572436/1‑150 (MQ=255)
                                                             cGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGt                                                                                   <  2:1522249/150‑1 (MQ=255)
                                                             cGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGt                                                                                   <  2:1098136/150‑1 (MQ=255)
                                                              gCCTTCTCTGCCCCCGGGGGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTa                                                                                  <  2:505839/150‑1 (MQ=255)
                                                                   ctctGCCCCCGGGGGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGt                                                                             <  2:1252998/150‑1 (MQ=255)
                                                                     ctGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGttt                                                                           >  1:44336/1‑150 (MQ=255)
                                                                       gCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTct                                                                         >  1:197728/1‑150 (MQ=255)
                                                                          cccGGGTGAAACTAGCGTAACGCTGATGCTGCCACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGCGCTGTTGGTACTGGTTTctgct                                                                      >  2:251010/1‑150 (MQ=255)
                                                                          cccGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCGGTTGGTACTGGTTTctgct                                                                      <  2:1653032/150‑1 (MQ=255)
                                                                           ccGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTctgctg                                                                     <  1:483069/150‑1 (MQ=255)
                                                                           ccGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTctgctg                                                                     >  1:1168973/1‑150 (MQ=255)
                                                                            cGGTTGAACCTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTgg                                                                    <  2:223711/150‑1 (MQ=255)
                                                                              ggTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCGGCTGGCa                                                                  >  1:848584/1‑150 (MQ=255)
                                                                               gTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCaa                                                                 >  1:539303/1‑150 (MQ=255)
                                                                                    aCTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTg                                                            <  2:521549/150‑1 (MQ=255)
                                                                                    aCTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTg                                                            <  2:610837/150‑1 (MQ=255)
                                                                                       aGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTTCCAGTTGGCATGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCTCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTc                                                         >  1:325198/1‑150 (MQ=255)
                                                                                          gTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCAtt                                                      >  1:487041/1‑150 (MQ=255)
                                                                                          gTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCAtt                                                      <  2:1439597/150‑1 (MQ=255)
                                                                                          gTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCAtt                                                      <  2:1432951/150‑1 (MQ=255)
                                                                                              cGCTGATGCTGACACATCCCGCCTGCGATTGCTGAATTCGTTTGACATGCCAGTTGCAAGGTGCGCAATTCCCGCGAGGTTGCTTTTGGTGCCCGCCATCCCGCCACACCGGTTGGACCTGGTTACGGCGGGAAAGCATGTACATTAATg                                                  >  2:26332/1‑150 (MQ=255)
                                                                                              cGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATg                                                  <  2:1624277/150‑1 (MQ=255)
                                                                                              cGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATg                                                  <  2:1526283/150‑1 (MQ=255)
                                                                                                    tGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCt                                            <  1:672807/150‑1 (MQ=255)
                                                                                                      cTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTgg                                          <  2:357353/150‑1 (MQ=255)
                                                                                                      cTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTgg                                          >  1:438949/1‑150 (MQ=255)
                                                                                                      cTGAAACATCCCGCGTGAGATGTCTGCATTCGGTTGACCTGCCAGTTGTCAGGTGGCGAAATGGCGCGGAGTTGCGTTTGCTGCCAGCCATGCCGCCAGATCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCttg                                          <  2:121697/150‑3 (MQ=255)
                                                                                                        gACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGaa                                        >  1:329982/1‑150 (MQ=255)
                                                                                                        gACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGaa                                        >  1:1707053/1‑150 (MQ=255)
                                                                                                        gACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGaa                                        <  1:1509080/150‑1 (MQ=255)
                                                                                                        gACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGaa                                        >  1:1266935/1‑150 (MQ=255)
                                                                                                        gACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGaa                                        >  1:538604/1‑150 (MQ=255)
                                                                                                          cacaTCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAAc                                      >  1:176661/1‑150 (MQ=255)
                                                                                                           acaTCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCCCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAAcc                                     <  1:1518308/150‑1 (MQ=255)
                                                                                                             aTCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCg                                   <  2:1059640/150‑1 (MQ=255)
                                                                                                             aTCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCg                                   >  1:762253/1‑150 (MQ=255)
                                                                                                             aTCCCGCCTGCGATGTCGGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCg                                   <  2:242422/150‑1 (MQ=255)
                                                                                                               cccGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGAt                                 <  1:420134/150‑1 (MQ=255)
                                                                                                                ccGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATa                                <  2:213067/150‑1 (MQ=255)
                                                                                                                 cGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATaa                               <  1:492652/150‑1 (MQ=255)
                                                                                                                       cGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGt                         <  2:913113/150‑1 (MQ=255)
                                                                                                                           gTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACa                     >  1:548131/1‑150 (MQ=255)
                                                                                                                             cTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAat                   <  1:397840/150‑1 (MQ=255)
                                                                                                                                    cGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCAcc            <  2:1707053/150‑1 (MQ=255)
                                                                                                                                    cGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCAcc            >  1:574992/1‑150 (MQ=255)
                                                                                                                                    cGGTGGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCAcc            <  2:1266935/150‑1 (MQ=255)
                                                                                                                                     ggTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCa           >  1:1261244/1‑150 (MQ=255)
                                                                                                                                     ggTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCa           >  2:1045361/1‑150 (MQ=255)
                                                                                                                                     ggTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGACAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTATGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCCCCa           >  2:81639/1‑150 (MQ=255)
                                                                                                                                      gTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCAt          <  2:1390555/150‑1 (MQ=255)
                                                                                                                                       gccaCCTGCCAACTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATc         >  1:98490/4‑150 (MQ=255)
                                                                                                                                       gccaCCTGCCAACTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATc         >  1:415214/4‑150 (MQ=255)
                                                                                                                                        ccaCCTGCCAACTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCa        >  2:1323651/3‑150 (MQ=255)
                                                                                                                                         gACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCaa       >  1:712758/1‑150 (MQ=255)
                                                                                                                                         gACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCaa       >  1:1038399/1‑150 (MQ=255)
                                                                                                                                          aCCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAac      >  1:1389292/1‑150 (MQ=255)
                                                                                                                                           ccTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAACCCCGATAACGCAGTGACAATCCTCACCATCAaca     >  1:899980/1‑150 (MQ=255)
                                                                                                                                           ccTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAaca     <  1:998121/150‑1 (MQ=255)
                                                                                                                                              gCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCa  <  1:1017154/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCACCATCAACACCA  >  NC_000913/2962415‑2962706

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