Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,742,796 A→C 5.1% E138D (GAA→GAC yfiN → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,742,7960AC 5.1% 381.8 / 7.3 128E138D (GAA→GACyfiNputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (63/58);  new base C (5/2);  total (68/60)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 130 total aligned reads displayed.

TTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAG  >  NC_000913/2742648‑2742944
                                                                                                                                                    |                                                                                                                                                     
tttCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGaaa                                                                                                                                                      <  2:1610410/150‑1 (MQ=255)
 ttCAACTGCAGAAGTAGGTGATTAGGTGGAATATATTGTGGCTTCCTGGCATTACACCCGTTAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCTCTGGCTCTTCGCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAAc                                                                                                                                                     <  1:408373/150‑1 (MQ=255)
 ttCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAAc                                                                                                                                                     >  1:1007310/1‑150 (MQ=255)
 ttCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAAc                                                                                                                                                     >  2:734949/1‑150 (MQ=255)
        gCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGc                                                                                                                                              <  1:1739217/150‑1 (MQ=255)
        gCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGc                                                                                                                                              >  2:3803/1‑150 (MQ=255)
          aGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGa                                                                                                                                            >  1:1797365/1‑150 (MQ=255)
           gAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTTTCACCATTgtgacg                                                                                                                                           <  2:541872/150‑6 (MQ=255)
           gAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGaa                                                                                                                                           <  1:1054168/150‑1 (MQ=255)
                aCGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACg                                                                                                                                      <  1:1428124/150‑1 (MQ=255)
                 cGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCTATTCGTCACAATGGTGAAACCATTGGCGAAGTACGc                                                                                                                                     <  2:1672177/150‑1 (MQ=255)
                  gTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCt                                                                                                                                    >  1:366806/1‑150 (MQ=255)
                  gTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCt                                                                                                                                    >  1:1482528/1‑150 (MQ=255)
                  gTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCt                                                                                                                                    <  1:612766/150‑1 (MQ=255)
                   tGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCtt                                                                                                                                   >  1:1543880/1‑150 (MQ=255)
                     aTAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTaa                                                                                                                                 <  2:1296869/150‑1 (MQ=255)
                       aaGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAAcc                                                                                                                               <  1:1370185/150‑1 (MQ=255)
                           agcaAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTc                                                                                                                           >  1:520095/1‑150 (MQ=255)
                            gcaAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTcg                                                                                                                          <  2:625159/150‑1 (MQ=255)
                            gcaAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCAGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCGTAACCGCTcg                                                                                                                          >  1:743054/1‑150 (MQ=255)
                              aaaaTATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCTGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTcgcg                                                                                                                        >  1:1541874/1‑150 (MQ=255)
                              aaaaTATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTcgcg                                                                                                                        >  1:1700026/1‑150 (MQ=255)
                              aaaaTATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTcgcg                                                                                                                        >  2:213985/1‑150 (MQ=255)
                               aaaTATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATACAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTCCGCTTAACCGCGCGCGa                                                                                                                       >  1:1131458/1‑150 (MQ=255)
                                  tatTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAg                                                                                                                    <  1:828543/150‑1 (MQ=255)
                                   atTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGt                                                                                                                   <  1:1690871/150‑1 (MQ=255)
                                     tCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTc                                                                                                                 <  1:722803/150‑1 (MQ=255)
                                       tGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCaa                                                                                                               <  2:363736/150‑1 (MQ=255)
                                       tGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGAAAGTTCaa                                                                                                               >  1:60992/1‑150 (MQ=255)
                                          cATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTATCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGACGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCa                                                                                                            >  2:1185841/1‑150 (MQ=255)
                                           aTCCTGGCATTACACCCGTAAGGAGCCAGGCGATACTTTCAGTAATTTCATAAGCCACTTGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAg                                                                                                           <  2:285691/150‑1 (MQ=255)
                                            tCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCGGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTAAATCAGc                                                                                                          >  1:150138/1‑150 (MQ=255)
                                             ccTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGcc                                                                                                         >  1:1135670/1‑150 (MQ=255)
                                               tGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCAt                                                                                                       >  2:668555/1‑150 (MQ=255)
                                                   aTTACACCCGTAAGGATCCAGGCGATAATTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCTATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTa                                                                                                   <  2:124882/150‑1 (MQ=255)
                                                      acacCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCTTTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTATCTg                                                                                               >  1:171539/1‑150 (MQ=255)
                                                      acacCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCt                                                                                                >  1:28636/1‑150 (MQ=255)
                                                       cacCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGAGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTg                                                                                               >  1:319887/1‑150 (MQ=255)
                                                        acCCGTAAGGATCCCGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCAAAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTgg                                                                                              <  2:1107376/150‑1 (MQ=255)
                                                        acCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTgg                                                                                              >  2:479673/1‑150 (MQ=255)
                                                         cccGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGt                                                                                             <  2:566952/150‑1 (MQ=255)
                                                         cccGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGt                                                                                             >  1:198436/1‑150 (MQ=255)
                                                         cccGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGt                                                                                             >  1:854929/1‑150 (MQ=255)
                                                         cccGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCCGGt                                                                                             >  2:600451/1‑150 (MQ=255)
                                                             tAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTtc                                                                                         <  1:1652837/150‑1 (MQ=255)
                                                             tAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTtc                                                                                         <  1:69215/150‑1 (MQ=255)
                                                              aaGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTtcg                                                                                        <  2:1051406/150‑1 (MQ=255)
                                                                ggATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTtcgct                                                                                      <  2:226149/150‑1 (MQ=255)
                                                                 gATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCCCGACAGTTCAATCAGCCATTTTATCTGGTTTtcgctc                                                                                     <  1:1011951/150‑1 (MQ=255)
                                                                  aTCCAGTCGAGACTTTAAGTAATTTCATAATCCACTGTCTCTTCACTGCCCCCATCATTCAGCCGATTCCTCTCAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTtcgctcg                                                                                    <  1:1549444/150‑1 (MQ=255)
                                                                  aTCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTtcgctcg                                                                                    <  1:1217079/150‑1 (MQ=255)
                                                                   tCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTtcgctcgc                                                                                   <  2:663030/150‑1 (MQ=255)
                                                                     cAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCg                                                                                 <  1:74299/150‑1 (MQ=255)
                                                                       ggCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCATTTTATCTGGTTTTCGCTCGCCGTAc                                                                              <  2:1373478/150‑1 (MQ=255)
                                                                          gATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCAGATTCGTCACAATGGTGAAACCATTGGCGAAGAACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTg                                                                            >  1:618752/1‑150 (MQ=255)
                                                                            tACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCTTTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTATCTGGTTTTCGCTCGCCGTACTGAcc                                                                         <  2:171539/150‑1 (MQ=255)
                                                                            tACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTgcc                                                                          >  1:33146/1‑148 (MQ=255)
                                                                             aCTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGAcc                                                                         <  1:393888/150‑1 (MQ=255)
                                                                                  cAGTAATTTCATAAGCAACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTTCGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCTCTCGCCGTACTGACCGGTTg                                                                    <  2:537131/150‑1 (MQ=255)
                                                                                    gTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGtt                                                                  >  1:1533049/1‑149 (MQ=255)
                                                                                     tAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTAt                                                                 <  1:216145/150‑1 (MQ=255)
                                                                                      aaTTTCATAAGCCACTGGCTCTTCCCTGCCCCCCTCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCGGGCATTTTCTCTGGTTTTCGCTCGCCGTTCGGACCGGGTGGAtt                                                                >  1:84751/1‑150 (MQ=255)
                                                                                         ttCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTctg                                                             <  2:239922/150‑1 (MQ=255)
                                                                                         ttCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTctg                                                             <  2:60992/150‑1 (MQ=255)
                                                                                            aTAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTctgctg                                                          >  2:1368414/1‑150 (MQ=255)
                                                                                              aaGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGc                                                        >  2:1502326/1‑150 (MQ=255)
                                                                                                  tccgatctctTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCAGTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCa                                                    >  2:1688892/7‑150 (MQ=255)
                                                                                                    cTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGCCAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGCCCGGTTGTATTCTGCTGGCATCAgg                                                  >  1:1207413/1‑150 (MQ=255)
                                                                                                    cTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAgg                                                  >  1:1013153/1‑150 (MQ=255)
                                                                                                      ggCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCa                                                >  1:268803/1‑150 (MQ=255)
                                                                                                              ccTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAAtt                                        >  1:547324/1‑150 (MQ=255)
                                                                                                              ccTGCCCCAATCATACAGCCGATTCGTCAAAATGGGGACACCATTGGCGAAGTACGCTTAACCGCTCGCCCCAGTGCAATCCGCCAGTTTATCTGGGTTTCCCACGCAGAACAGACCGGTTGTAATCTGCTGGAATCAAGCATCGCAAtt                                        >  1:1611434/1‑150 (MQ=255)
                                                                                                               cTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCTTCGCAATTa                                       <  1:973565/150‑1 (MQ=255)
                                                                                                                 gCCCCCATCATTAAGCCGATTAGTCACAAAGGTCAAGCCAATGGCGAAGTACGCTTACCCGGTCGCGACAGTTAAATCCGCCATTATATCTGGTTTTCGCTCGCCGTGCTGCCCGGTTGTATTCTGCTGGCAACAGGCATCGCAATTAcc                                     >  2:105011/1‑150 (MQ=255)
                                                                                                                  cccccATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTAccc                                    <  2:896646/150‑1 (MQ=255)
                                                                                                                   ccccATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCt                                   <  2:1520281/150‑1 (MQ=255)
                                                                                                                   ccccATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCt                                   >  2:849446/1‑150 (MQ=255)
                                                                                                                    cccATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACAGCTCGCGACAGTTCAATCAGCCATGTTATCTGGTTTTCGCTCGCCGTCCTGACCGTATGTATTCTGCTGGCTTCAGGCAACGCAATTACCCTc                                  >  2:1228004/1‑150 (MQ=255)
                                                                                                                     ccATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCCGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCAGCAGGCATCGCAATTACCCTcc                                 >  2:1672363/1‑149 (MQ=255)
                                                                                                                     ccATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCa                                 <  2:1428385/150‑1 (MQ=255)
                                                                                                                       atcatTCAGCCGTTTCGTCACAATGGTGAATCCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCAcc                               <  2:1587587/150‑1 (MQ=255)
                                                                                                                          atTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTATCCTCACCCGc                            <  2:1541874/150‑1 (MQ=255)
                                                                                                                              aGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAATCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCAtt                        >  2:1456728/1‑150 (MQ=255)
                                                                                                                              aGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCAtt                        >  1:1447778/1‑150 (MQ=255)
                                                                                                                                ccGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTg                     <  2:366806/150‑1 (MQ=255)
                                                                                                                                   aTTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGcac                  <  1:1502326/150‑1 (MQ=255)
                                                                                                                                       gTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATgg             >  2:161041/1‑150 (MQ=255)
                                                                                                                                        tCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATgg             >  2:103533/1‑150 (MQ=255)
                                                                                                                                         gccaATGGTTTCACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGc            >  1:1237888/3‑150 (MQ=255)
                                                                                                                                         gccaATGGTTTCACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTAATCTGCTGGCATCAGGCATTGCAATTACCCTCACCCGCCATTTGCACAATGGc            >  2:1664316/3‑150 (MQ=255)
                                                                                                                                         gccaATGGTTTCACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGCTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACTATGGc            >  1:74232/3‑150 (MQ=255)
                                                                                                                                           caATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTGTGCAAATGGcgg          <  1:715963/150‑3 (MQ=255)
                                                                                                                                           caATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCt          <  1:216826/150‑1 (MQ=255)
                                                                                                                                           caATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCt          >  2:1618610/1‑150 (MQ=255)
                                                                                                                                            aaTGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTCTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTGTGCAAATGGcggg         <  1:1103420/150‑4 (MQ=255)
                                                                                                                                               ggTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTa      <  1:503474/150‑1 (MQ=255)
                                                                                                                                                gTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAg     <  1:727576/150‑1 (MQ=255)
                                                                                                                                                 tGAAACCATTGGCGAAGTACACTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTAGCTCTCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGa    >  2:892044/1‑150 (MQ=255)
                                                                                                                                                  gAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGaa   >  2:1779139/1‑150 (MQ=255)
                                                                                                                                                   aaaCCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAg  <  2:1416359/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAG  >  NC_000913/2742648‑2742944

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: