Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,561,205 A→T 5.3% intergenic (+307/‑163) yffL → / → yffM CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,561,2050AT 5.3% 332.9 / 7.9 94intergenic (+307/‑163)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (48/41);  new base T (2/3);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

TATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTA  >  NC_000913/2561057‑2561351
                                                                                                                                                    |                                                                                                                                                  
taTTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCa                                                                                                                                                    >  1:445279/1‑149 (MQ=255)
   tAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGc                                                                                                                                                 <  1:460895/149‑1 (MQ=255)
    aaGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCt                                                                                                                                                >  2:563160/1‑149 (MQ=255)
     aGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAg                                                                                                                                                  <  1:74797/146‑1 (MQ=255)
     aGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAg                                                                                                                                                  >  4:119861/1‑146 (MQ=255)
     aGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAg                                                                                                                                                  <  3:119861/146‑1 (MQ=255)
     aGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAg                                                                                                                                                  >  2:74797/1‑146 (MQ=255)
     aGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTg                                                                                                                                               >  1:760294/1‑149 (MQ=255)
          tCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAAc                                                                                                                                          <  1:337312/149‑1 (MQ=255)
          tCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAAc                                                                                                                                          >  3:618910/1‑149 (MQ=255)
              ggCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTTGCTGGaaatgatt                                                                                                                                      <  2:300639/149‑6 (MQ=255)
              ggCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTTGCTGGaaatgatt                                                                                                                                      <  2:243174/149‑6 (MQ=255)
                   cGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtca                                                                                                                                 <  4:618910/149‑1 (MQ=255)
                    gggAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcac                                                                                                                                <  4:770017/149‑1 (MQ=255)
                    gggAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcac                                                                                                                                <  2:708867/149‑1 (MQ=255)
                      gAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcactc                                                                                                                              <  1:753284/149‑1 (MQ=255)
                      gAGCTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGCTAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcactc                                                                                                                              >  1:768021/1‑149 (MQ=255)
                          tttttAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGc                                                                                                                          >  1:149455/1‑149 (MQ=255)
                           ttttAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGca                                                                                                                         <  3:164028/149‑1 (MQ=255)
                                  cTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGCATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTg                                                                                                                  >  4:586555/1‑149 (MQ=255)
                                    cccGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGaa                                                                                                                <  1:654473/149‑1 (MQ=255)
                                     ccGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAc                                                                                                               <  3:703370/149‑1 (MQ=255)
                                     ccGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAc                                                                                                               <  1:10871/149‑1 (MQ=255)
                                        gtcggtcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagc                                                                                                            >  3:756640/1‑149 (MQ=255)
                                            gtcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGAc                                                                                                        >  2:472150/1‑149 (MQ=255)
                                             tcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGAcc                                                                                                       >  2:581940/1‑149 (MQ=255)
                                                tAGAGAGTCGCGTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACcaac                                                                                                    >  3:691741/1‑149 (MQ=255)
                                                 agagagTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAAc                                                                                                                     <  3:301701/131‑1 (MQ=255)
                                                 agagagTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAAc                                                                                                                     >  4:301701/1‑131 (MQ=255)
                                                 agagagTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACACCCATTCACTCACGCACAACCTGAACAGCAGACcaaca                                                                                                   >  4:61336/1‑149 (MQ=255)
                                                  gagagTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACcaacaa                                                                                                  >  2:561043/1‑149 (MQ=255)
                                                    gagTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAAc                                                                                                >  2:559149/1‑149 (MQ=255)
                                                     agTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAAcc                                                                                               >  4:559400/1‑149 (MQ=255)
                                                      gTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACca                                                                                              <  3:78383/149‑1 (MQ=255)
                                                            tAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATgg                                                                                        >  4:328058/1‑149 (MQ=255)
                                                              gagTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCa                                                                                      >  3:662824/1‑149 (MQ=255)
                                                                 ttAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACcaac                                                                                                    >  1:794306/1‑132 (MQ=255)
                                                                 ttAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACcaac                                                                                                    <  2:794306/132‑1 (MQ=255)
                                                                 ttAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTc                                                                                   >  4:465206/1‑149 (MQ=255)
                                                                 ttAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTc                                                                                   >  2:215994/1‑149 (MQ=255)
                                                                     aCAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTc                                                                               <  4:46169/149‑1 (MQ=255)
                                                                       agaagaATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAAc                                                                                                >  4:441977/1‑130 (MQ=255)
                                                                       agaagaATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAt                                                                             >  1:793767/1‑149 (MQ=255)
                                                                       agaagaATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCACATTCACTCACGCACAACCTGAACAGCAGACCAACAAAc                                                                                                <  3:441977/130‑1 (MQ=255)
                                                                           gaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATgggg                                                                         >  1:421554/1‑149 (MQ=255)
                                                                            aaTTATTTCGCGGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGTGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGTAGAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGa                                                                        >  4:72259/1‑149 (MQ=255)
                                                                            aaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGa                                                                        <  4:756640/149‑1 (MQ=255)
                                                                              ttattTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAgc                                                                      >  3:615971/1‑149 (MQ=255)
                                                                                  ttCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCgg                                                                  <  2:445279/149‑1 (MQ=255)
                                                                                    cgccgcCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGta                                                                <  1:559149/149‑1 (MQ=255)
                                                                                       cgccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTatat                                                             <  3:715581/149‑1 (MQ=255)
                                                                                        gccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATAtc                                                            >  4:12965/1‑149 (MQ=255)
                                                                                         ccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATAtct                                                           <  1:581940/149‑1 (MQ=255)
                                                                                             aTAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACTAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCaaa                                                       <  4:662824/149‑1 (MQ=255)
                                                                                             aTAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCaaa                                                       <  3:12965/149‑1 (MQ=255)
                                                                                              tAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAg                                                      <  1:561043/149‑1 (MQ=255)
                                                                                               aGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAgg                                                     >  1:262341/1‑149 (MQ=255)
                                                                                               aGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAgg                                                     <  1:563160/149‑1 (MQ=255)
                                                                                                   cGGCTGGTTGCGAGCATTAACCATTCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGTCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCTTCATGGCATTCGTTCATGGGTAGCGCGGTATATCTCAAAGGGTGa                                                 <  2:768021/149‑1 (MQ=255)
                                                                                                    ggCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGAt                                                >  4:460006/1‑149 (MQ=255)
                                                                                                      cTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGAttt                                              <  3:465206/149‑1 (MQ=255)
                                                                                                       tGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTc                                             <  1:215994/149‑1 (MQ=255)
                                                                                                              gAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtgg                                      >  2:60328/1‑149 (MQ=255)
                                                                                                                 catcatccATGCTGATTAGCGCAATCACTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtggtgg                                   <  4:691741/143‑1 (MQ=255)
                                                                                                                 cATTAACCATGCTCATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAACACATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtgttgg                                   <  4:507141/149‑1 (MQ=255)
                                                                                                                  aTTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTggg                                  >  2:814502/1‑149 (MQ=255)
                                                                                                                     aaCCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCt                               >  3:118661/1‑149 (MQ=255)
                                                                                                                     aaCCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTCGGCCt                               >  3:417852/1‑149 (MQ=255)
                                                                                                                      aCCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCtt                              >  2:765727/1‑149 (MQ=255)
                                                                                                                       ccATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTc                             <  3:559400/149‑1 (MQ=255)
                                                                                                                       ccATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTc                             <  1:413519/149‑1 (MQ=255)
                                                                                                                        cATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCg                            >  2:218213/1‑149 (MQ=255)
                                                                                                                         aTGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCgg                           >  4:79870/1‑149 (MQ=255)
                                                                                                                          tGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGa                          >  1:440607/1‑149 (MQ=255)
                                                                                                                          tGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGa                          >  2:70586/1‑149 (MQ=255)
                                                                                                                          tGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGa                          >  4:381269/1‑149 (MQ=255)
                                                                                                                           gCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGaa                         <  3:79870/149‑1 (MQ=255)
                                                                                                                           gCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGaa                         <  4:615971/149‑1 (MQ=255)
                                                                                                                             tGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAgg                       >  4:224860/1‑149 (MQ=255)
                                                                                                                               aTTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCt                     >  4:389367/1‑149 (MQ=255)
                                                                                                                                 tAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTc                   <  2:262341/149‑1 (MQ=255)
                                                                                                                                   gcgcAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCaaa                                                       >  4:381246/1‑111 (MQ=255)
                                                                                                                                   gcgcAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCaaa                                                       <  3:381246/111‑1 (MQ=255)
                                                                                                                                   gcgcAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGa                 <  3:61336/149‑1 (MQ=255)
                                                                                                                                   gcgcAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGa                 <  3:72259/149‑1 (MQ=255)
                                                                                                                                        aTCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAAt            <  3:328058/149‑1 (MQ=255)
                                                                                                                                        aTCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAAt            >  2:240487/1‑149 (MQ=255)
                                                                                                                                         tCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAAtt           >  2:564471/1‑149 (MQ=255)
                                                                                                                                           gtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGc         >  2:374883/2‑149 (MQ=255)
                                                                                                                                            tttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtggtg                                    >  4:564131/1‑121 (MQ=255)
                                                                                                                                            tttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtggtg                                    <  3:564131/121‑1 (MQ=255)
                                                                                                                                               ccAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTc                             >  2:496942/1‑125 (MQ=255)
                                                                                                                                               ccAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTc                             <  1:496942/125‑1 (MQ=255)
                                                                                                                                               ccAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCCGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGc     >  2:105745/1‑149 (MQ=255)
                                                                                                                                                  gCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTa  >  3:204353/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTA  >  NC_000913/2561057‑2561351

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