Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,315,439 A→T 5.9% F34I (TTT→ATT) ‡ phnO ← aminoalkylphosphonate N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,315,4390AT 5.9% 717.9 / 23.7 232F34I (TTT→ATT) ‡phnOaminoalkylphosphonate N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (101/116);  new base T (3/11);  total (105/127)
Fisher's exact test for biased strand distribution p-value = 9.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

Only 100 of 236 total aligned reads displayed.

CACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTG  >  NC_000913/4315290‑4315587
                                                                                                                                                     |                                                                                                                                                    
     aCTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCgt                                                                                                                                                 >  2:1649910/1‑150 (MQ=255)
      cTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCgtg                                                                                                                                                <  1:1612000/150‑1 (MQ=255)
        ccTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAATTGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCgtggt                                                                                                                                              >  2:1576890/1‑150 (MQ=255)
        ccTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGTTTCGCgtggg                                                                                                                                              <  1:779069/150‑2 (MQ=255)
          tGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCTAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGtggtgg                                                                                                                                            <  2:1416346/150‑1 (MQ=255)
          tGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGTTTCGCgtgggtt                                                                                                                                            <  1:373606/150‑4 (MQ=255)
          tGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGtggtgg                                                                                                                                            <  2:1201540/150‑1 (MQ=255)
          tGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGtggtgg                                                                                                                                            >  1:1201540/1‑150 (MQ=255)
           gAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGTTTCGCgtgggttt                                                                                                                                           <  2:417658/150‑5 (MQ=255)
           gAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTggtggt                                                                                                                                           <  1:484221/150‑1 (MQ=255)
                tCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAAc                                                                                                                                      <  1:686779/150‑1 (MQ=255)
                tCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAAc                                                                                                                                      <  1:652842/150‑1 (MQ=255)
                 cGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACt                                                                                                                                     >  1:598733/1‑150 (MQ=255)
                  gCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTc                                                                                                                                    >  1:1949364/1‑150 (MQ=255)
                   ccGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTcc                                                                                                                                   >  1:1282403/1‑150 (MQ=255)
                    cGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGTGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCg                                                                                                                                  <  1:1788327/150‑1 (MQ=255)
                      aTCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGcc                                                                                                                                <  2:1793930/150‑1 (MQ=255)
                             tGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTa                                                                                                                         >  2:642469/1‑150 (MQ=255)
                             tGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTa                                                                                                                         <  2:508627/150‑1 (MQ=255)
                             tGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTa                                                                                                                         <  2:1736136/150‑1 (MQ=255)
                              gACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAg                                                                                                                        <  2:1282403/150‑1 (MQ=255)
                                cATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGTCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCt                                                                                                                      <  1:270907/150‑1 (MQ=255)
                                   gatgCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTcac                                                                                                                   <  1:1025389/150‑1 (MQ=255)
                                     tgCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAGCTCCGCCTGTTTTAGCTCACaa                                                                                                                 >  2:1363841/1‑150 (MQ=255)
                                     tgCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACaa                                                                                                                 >  1:327225/1‑150 (MQ=255)
                                      gCAGATGAAACTGCAAATGCAGGTCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACaaa                                                                                                                >  1:148526/1‑150 (MQ=255)
                                      gCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAAATTAGCCATCAAGCAGAGCCAGATGGTAGCGCATGTTTGGGTCGAGCAGATTCGCGTTAAAACCCACGCGAAAAGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACaaa                                                                                                                >  2:849472/1‑150 (MQ=255)
                                         gATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCa                                                                                                             >  2:797105/1‑150 (MQ=255)
                                          aTGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAg                                                                                                            <  1:1846836/150‑1 (MQ=255)
                                            gAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGTTGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAgcg                                                                                                          >  2:559542/1‑150 (MQ=255)
                                                cTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTa                                                                                                      >  1:1486545/1‑150 (MQ=255)
                                                 tGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTaa                                                                                                     >  2:378327/1‑150 (MQ=255)
                                                 tGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGAGGTGGTCAACCTCCGCCTGTTTTAGCTCACAAATCAGCGCGTaa                                                                                                     >  2:535877/1‑150 (MQ=255)
                                                  gCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGTGTCGCGCAGATTGGCTTTAAAAACCACGCGAATCGCGTGGTTGTCTAACTGCGCCTGTTTTAGTTCACAAATCTTCGCGTaat                                                                                                    <  2:1685142/150‑2 (MQ=255)
                                                   cAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAAc                                                                                                   >  2:1548001/1‑150 (MQ=255)
                                                   cAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAAc                                                                                                   <  2:656044/150‑1 (MQ=255)
                                                       tGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACcgcg                                                                                               >  2:1431247/1‑150 (MQ=255)
                                                        gCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGt                                                                                              >  1:287371/1‑150 (MQ=255)
                                                            gCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGgt                                                                                          >  1:1217604/1‑150 (MQ=255)
                                                              cGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGgtgt                                                                                        <  2:1218312/150‑1 (MQ=255)
                                                                     gCCGACATCTTCGCCATCAAGCAGCGCGAGTTGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACGCAGGCGATACGGGTGGTGGTGATACTGCTCCTGTTTTAGCTGACAAATCAGCGCGTAAACCGCGTGTGTGTCGTACTg                                                                                 <  2:903297/150‑1 (MQ=255)
                                                                      ccGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGc                                                                                >  2:680746/1‑150 (MQ=255)
                                                                           aaCTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGc                                                                           <  2:1524094/150‑1 (MQ=255)
                                                                             cTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCg                                                                         >  1:1248311/1‑150 (MQ=255)
                                                                             cTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCg                                                                         >  2:1493203/1‑150 (MQ=255)
                                                                              ttCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCgg                                                                        >  2:1325335/1‑150 (MQ=255)
                                                                               tCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCggg                                                                       <  1:141835/150‑1 (MQ=255)
                                                                               tCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCggg                                                                       <  1:1889789/150‑1 (MQ=255)
                                                                               tCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCggg                                                                       <  1:1967278/150‑1 (MQ=255)
                                                                                  ccATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGa                                                                    >  2:24986/1‑150 (MQ=255)
                                                                                  ccATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGa                                                                    <  2:775750/150‑1 (MQ=255)
                                                                                  ccATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGa                                                                    >  1:1995320/1‑150 (MQ=255)
                                                                                   cATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGaa                                                                   >  2:1641821/1‑150 (MQ=255)
                                                                                    aTCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAg                                                                  <  1:1641821/150‑1 (MQ=255)
                                                                                     tCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGc                                                                 <  2:148526/150‑1 (MQ=255)
                                                                                     tCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGc                                                                 <  1:1193162/150‑1 (MQ=255)
                                                                                       aaGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTc                                                               <  1:1351217/150‑1 (MQ=255)
                                                                                        agcagcGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTca                                                              <  2:751999/150‑1 (MQ=255)
                                                                                         gcagcGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTcac                                                             >  1:1059306/1‑150 (MQ=255)
                                                                                         gcagcGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTcac                                                             >  1:5622/1‑150 (MQ=255)
                                                                                            gcgcCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAg                                                          <  1:859580/150‑1 (MQ=255)
                                                                                             cgcCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGc                                                         <  1:1313449/150‑1 (MQ=255)
                                                                                              gcCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCa                                                        <  1:68260/150‑1 (MQ=255)
                                                                                              gcCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCa                                                        <  1:1441604/150‑1 (MQ=255)
                                                                                                cAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAgg                                                      >  1:332648/1‑150 (MQ=255)
                                                                                                 aGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGc                                                     <  1:411979/150‑1 (MQ=255)
                                                                                                       tAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGGCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGgtg                                               <  2:85854/150‑1 (MQ=255)
                                                                                                               gTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGttt                                            >  1:535390/1‑145 (MQ=255)
                                                                                                                  tGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCt                                    <  1:833565/150‑1 (MQ=255)
                                                                                                                  tGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCt                                    <  1:424080/150‑1 (MQ=255)
                                                                                                                   gggTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGGGTTTCTCCTTCTg                                   <  2:1949364/150‑1 (MQ=255)
                                                                                                                    ggTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGa                                  >  2:144486/1‑150 (MQ=255)
                                                                                                                     gTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGAt                                 <  2:941793/150‑1 (MQ=255)
                                                                                                                       cgcgCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATgg                               <  1:996291/150‑1 (MQ=255)
                                                                                                                       cgcgCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATgg                               >  2:521201/1‑150 (MQ=255)
                                                                                                                        gcgcAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGa                              >  1:1673361/1‑150 (MQ=255)
                                                                                                                          gcAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATc                            <  1:1320992/150‑1 (MQ=255)
                                                                                                                          gcAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATc                            <  1:725264/150‑1 (MQ=255)
                                                                                                                            aGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAg                          >  1:1926859/1‑150 (MQ=255)
                                                                                                                              aTTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCg                        <  2:278010/150‑1 (MQ=255)
                                                                                                                              aTTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCg                        <  1:1081263/150‑1 (MQ=255)
                                                                                                                               ttGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGt                       >  2:552302/1‑150 (MQ=255)
                                                                                                                                tGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTc                      >  1:1898302/1‑150 (MQ=255)
                                                                                                                                 ggCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTca                     <  1:96294/150‑1 (MQ=255)
                                                                                                                                 ggCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTca                     <  1:627232/150‑1 (MQ=255)
                                                                                                                                  gCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTcag                    >  2:836363/1‑150 (MQ=255)
                                                                                                                                   cGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTcac                                                                                                                   >  1:1259857/1‑54 (MQ=255)
                                                                                                                                   cGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGc                                                                      >  2:454018/1‑99 (MQ=255)
                                                                                                                                       aaaaCCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCg               <  1:1694074/150‑1 (MQ=255)
                                                                                                                                       aaaaCCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCg               >  1:93363/1‑150 (MQ=255)
                                                                                                                                         aaCCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCg               >  2:1080076/1‑148 (MQ=255)
                                                                                                                                          caccACGCGTTTCGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCGGCAGCgtgt            >  1:1416347/3‑150 (MQ=255)
                                                                                                                                          aCCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCgtgt            >  1:785918/1‑150 (MQ=255)
                                                                                                                                            ccACGCGTTTCGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTc           <  1:1569636/149‑1 (MQ=255)
                                                                                                                                             cACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTcga         >  2:767691/1‑150 (MQ=255)
                                                                                                                                              aCGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTcgac        <  1:521201/150‑1 (MQ=255)
                                                                                                                                               cgcgAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTcgacc       >  1:218579/1‑150 (MQ=255)
                                                                                                                                                  gAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTcgaccgac    <  2:327225/150‑1 (MQ=255)
                                                                                                                                                   aaaCGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACt   >  2:1426399/1‑150 (MQ=255)
                                                                                                                                                    aaCGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTg  <  2:598733/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTG  >  NC_000913/4315290‑4315587

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: