Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,531,118 Δ1 bp 7.3% coding (534/1122 nt) sgcX ← KpLE2 phage‑like element; putative endoglucanase with Zn‑dependent exopeptidase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,531,1180T. 7.3% 266.6 / 13.6 82coding (534/1122 nt)sgcXKpLE2 phage‑like element; putative endoglucanase with Zn‑dependent exopeptidase domain
Reads supporting (aligned to +/- strand):  ref base T (35/41);  new base . (3/3);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACTGCG  >  NC_000913/4530974‑4531266
                                                                                                                                                 |                                                                                                                                                    
tACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCttt                                                                                                                                                 >  1:120777/1‑150 (MQ=255)
    cGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTg                                                                                                                                             >  1:495768/1‑150 (MQ=255)
    cGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTg                                                                                                                                             >  1:905311/1‑150 (MQ=255)
     ggCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGATCATGGCATTACAGCCTAGACGATCGTCTAGTGCTTTACTGc                                                                                                                                            <  2:227713/150‑1 (MQ=255)
      gCGTAATACGGTAACAATGCCGCGGATATTAAATTCTTCCTGTACCTAAGCTACCAGATAAACGGCGATATCAAGCTACATAGTACTGATAGCATCAGCTACACCGAGCAGGGCAGTACATCCTAGACGATCGTCTAGTGCTTTACTGca                                                                                                                                           <  2:67296/150‑1 (MQ=255)
      gCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGca                                                                                                                                           <  1:19360/150‑1 (MQ=255)
        gTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGccct                                                                                                                                         <  1:858216/150‑4 (MQ=38)
         tAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGccctg                                                                                                                                        <  1:1102511/150‑5 (MQ=38)
          aaTACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                                                                       >  2:850837/1‑150 (MQ=255)
                ggAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCg                                                                                                                                 <  1:473915/150‑1 (MQ=255)
                ggAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCg                                                                                                                                 >  1:1403402/1‑150 (MQ=255)
                   aCAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGttt                                                                                                                              <  2:161839/150‑1 (MQ=255)
                        gCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAg                                                                                                                         <  2:446000/150‑1 (MQ=255)
                                tatTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcg                                                                                                                 <  2:750899/150‑1 (MQ=255)
                                   tAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcggcg                                                                                                              >  2:968755/1‑150 (MQ=255)
                                    aaaTTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGcggcgg                                                                                                             <  2:1000950/150‑1 (MQ=255)
                                     aaTTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCggg                                                                                                            <  2:92497/150‑1 (MQ=255)
                                     aaTTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCggg                                                                                                            >  2:206049/1‑150 (MQ=255)
                                     aaTTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCggg                                                                                                            <  1:531764/150‑1 (MQ=255)
                                      aTTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGt                                                                                                           <  2:391790/150‑1 (MQ=255)
                                          ttcCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTtata                                                                                                       <  2:58676/150‑1 (MQ=255)
                                          ttcCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTtata                                                                                                       >  2:778920/1‑150 (MQ=255)
                                           tcCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATAc                                                                                                      >  1:681354/1‑150 (MQ=255)
                                            ccTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACa                                                                                                     >  1:1024313/1‑150 (MQ=255)
                                             cTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACGCCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACa                                                                                                     >  1:267921/1‑149 (MQ=255)
                                             cTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACag                                                                                                    >  2:911580/1‑150 (MQ=255)
                                             cTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACa                                                                                                     <  2:267921/149‑1 (MQ=255)
                                              tGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAg                                                                                                                         >  1:537450/1‑128 (MQ=255)
                                              tGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAg                                                                                                                         <  2:537450/128‑1 (MQ=255)
                                                tACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGt                                                                                                 <  1:1081484/150‑1 (MQ=255)
                                                 aCCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTg                                                                                                >  2:1282464/1‑150 (MQ=255)
                                                  ccGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGa                                                                                               <  2:120777/150‑1 (MQ=255)
                                                   cGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGAc                                                                                              <  2:215727/150‑1 (MQ=255)
                                                      aGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAgg                                                                                           <  1:870143/150‑1 (MQ=255)
                                                      aGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAgg                                                                                           >  1:596639/1‑150 (MQ=255)
                                                      aGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAgg                                                                                           <  1:1077789/150‑1 (MQ=255)
                                                        ccaccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCg                                                                                         <  1:333591/150‑1 (MQ=255)
                                                         caccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGt                                                                                        >  1:847805/1‑150 (MQ=255)
                                                                tAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAcc                                                                                 >  1:849603/1‑150 (MQ=255)
                                                                  aaCGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTCCACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                                                                               >  1:112196/1‑150 (MQ=255)
                                                                  aaCGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                                                                               <  2:495768/150‑1 (MQ=255)
                                                                        gATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGcc                                                                         <  2:681354/150‑1 (MQ=255)
                                                                          tatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACGCCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCa                                                                       >  1:1145027/1‑150 (MQ=255)
                                                                          tatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCa                                                                       >  1:415818/1‑150 (MQ=255)
                                                                            tCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATc                                                                     <  2:525514/150‑1 (MQ=255)
                                                                            tCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATc                                                                     <  2:1403402/150‑1 (MQ=255)
                                                                                 cTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGtat                                                                <  2:415818/150‑1 (MQ=255)
                                                                                 cTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACATGAATGCCCATGCGtat                                                                <  2:192051/150‑1 (MQ=255)
                                                                                   ccATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAg                                                              >  1:994397/1‑150 (MQ=255)
                                                                                    cATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGc                                                             >  1:923/1‑150 (MQ=255)
                                                                                      tAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGcgt                                                           <  1:206049/150‑1 (MQ=255)
                                                                                       aGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGcgtc                                                          >  2:1002052/1‑150 (MQ=255)
                                                                                       aGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGcgtc                                                          <  2:596639/150‑1 (MQ=255)
                                                                                            tGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctt                                                     >  1:354233/1‑150 (MQ=255)
                                                                                                 gCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTAAAGCGTCgtctttgtct                                                >  1:660431/1‑150 (MQ=255)
                                                                                                    tCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGtctttgtctttg                                             <  1:968755/150‑1 (MQ=255)
                                                                                                       gCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCa                                          >  1:1429505/1‑150 (MQ=255)
                                                                                                        cTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCAc                                         <  2:994397/150‑1 (MQ=255)
                                                                                                             ccGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAata                                    >  2:1079816/1‑150 (MQ=255)
                                                                                                               gAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTTAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATc                                  <  2:1236928/150‑1 (MQ=255)
                                                                                                               gAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATc                                  <  1:911580/150‑1 (MQ=255)
                                                                                                                   tgtgCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATc                              >  1:659870/4‑150 (MQ=38)
                                                                                                                      gCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCcac                           >  1:77602/1‑150 (MQ=255)
                                                                                                                      gCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCcac                           >  2:1406206/1‑150 (MQ=39)
                                                                                                                           aCAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGttt                      >  1:798709/1‑150 (MQ=255)
                                                                                                                           aCAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGttt                      >  1:1285667/1‑150 (MQ=255)
                                                                                                                             aGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGt                    <  2:847805/150‑1 (MQ=255)
                                                                                                                              gcacTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTc                   <  1:1406206/147‑1 (MQ=255)
                                                                                                                              gCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTc                   <  1:1002052/150‑1 (MQ=255)
                                                                                                                              gCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTc                   <  2:1024313/150‑1 (MQ=255)
                                                                                                                              gCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTc                   >  2:633519/1‑150 (MQ=255)
                                                                                                                                 cTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCGTTGGCACCAATATCAATCCACAGTTTGTCGa                 >  1:531177/1‑150 (MQ=255)
                                                                                                                                 cTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGAc                >  1:1132705/1‑150 (MQ=255)
                                                                                                                                  tAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCt                                                                                <  2:1294798/86‑1 (MQ=255)
                                                                                                                                  tAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCt                                                                                >  1:1294798/1‑86 (MQ=255)
                                                                                                                                  tAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGAc                <  2:849603/150‑1 (MQ=255)
                                                                                                                                        aTCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAgg          <  2:1429505/150‑1 (MQ=255)
                                                                                                                                        aTCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAgg          >  2:480983/1‑150 (MQ=255)
                                                                                                                                         tCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTAAAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGt         <  2:660431/150‑1 (MQ=255)
                                                                                                                                            tCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGAc      <  1:850837/150‑1 (MQ=255)
                                                                                                                                                gTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACtgcg  >  1:1372398/1‑150 (MQ=255)
                                                                                                                                                gTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACtgcg  <  1:1282464/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACTGCG  >  NC_000913/4530974‑4531266

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