Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,763,609 A→G 5.5% V134A (GTG→GCG)  sufB ← component of SufBCD Fe‑S cluster assembly scaffold

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,763,6090AG 5.5% 255.9 / 3.8 90V134A (GTG→GCG) sufBcomponent of SufBCD Fe‑S cluster assembly scaffold
Reads supporting (aligned to +/- strand):  ref base A (40/45);  new base G (2/3);  total (42/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

TATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGATGGTGCTGAGTAGTAGCTGTAATCCTGATA  >  NC_000913/1763462‑1763757
                                                                                                                                                   |                                                                                                                                                    
tatTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACg                                                                                                                                                    <  2:527977/150‑1 (MQ=255)
tatTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACg                                                                                                                                                    >  1:271226/1‑150 (MQ=255)
 atTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACgg                                                                                                                                                   >  1:486766/1‑150 (MQ=255)
    tACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAATGAACAGAAAATAATTCCCTGCTCCGCCATTTTTTCGCGATAAGTAGTGGCAACCGTAACTTAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTACCGCACGGGaa                                                                                                                                                <  1:50495/150‑1 (MQ=255)
      cGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACg                                                                                                                                              >  1:598384/1‑150 (MQ=255)
      cGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACg                                                                                                                                              <  1:337643/150‑1 (MQ=255)
       gCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGc                                                                                                                                             >  2:1017596/1‑150 (MQ=255)
        cACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCTATAAGTAGTGGCAACCGAAACTGAGTCGAAGATGGCATACACCGCCCCCGATTTGCCTTCCCGCACGGTAAAGcc                                                                                                                                            >  2:403093/1‑150 (MQ=255)
        cACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGTGCGGGAAGGca                                                                                                                                            <  1:1378366/150‑2 (MQ=255)
           cAGTTCCGTGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCTACCGAAACTGTGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGTGCGGGAAGGcaaag                                                                                                                                         <  2:1044852/150‑5 (MQ=255)
           cAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGCAACGCCCaa                                                                                                                                         >  2:1624039/1‑150 (MQ=255)
            aGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGTGCGGGAAGGcaaaga                                                                                                                                        <  2:1211415/150‑6 (MQ=255)
                ccGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCt                                                                                                                                    <  2:1628748/150‑1 (MQ=255)
                ccGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCt                                                                                                                                    <  2:1711420/150‑1 (MQ=255)
                      gATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAAcg                                                                                                                              <  1:1374254/150‑1 (MQ=255)
                       aTCGTGGATCGCCTCACCAAAGGAACAGAAAATACTTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCATCCGAAACTGAGTCGCAACTGGCAGCCACCGCCACCTCTTTGCCTGCCCGCACGGGCTCGCCCAACGGCTCATAcgc                                                                                                                             >  1:269094/1‑150 (MQ=255)
                        tCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAAcgcc                                                                                                                            >  1:1619370/1‑150 (MQ=255)
                          gTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAAcgcctt                                                                                                                          >  1:879300/1‑148 (MQ=255)
                          gTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAAcgccgc                                                                                                                          >  1:1265700/1‑150 (MQ=255)
                            ggATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCt                                                                                                                        >  2:1618240/1‑150 (MQ=255)
                             gATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTc                                                                                                                       >  2:1166887/1‑150 (MQ=255)
                              aTCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTcc                                                                                                                      <  2:822133/150‑1 (MQ=255)
                              aTCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTcc                                                                                                                      <  2:1717217/150‑1 (MQ=255)
                                 gCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCAcc                                                                                                                   >  1:361328/1‑150 (MQ=255)
                                      aCCAAAGGAACAGAAAAAAATTACCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACAGAAACTGAGTCGAAAATGGCATCCACCGCAACCTCTTTGCATTCCAGCACGGGAACTACCAACTGCTCAAACGCCGCATCCAcatatat                                                                                                              >  2:1587166/1‑144 (MQ=255)
                                       ccAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTa                                                                                                             >  2:1550481/1‑150 (MQ=255)
                                        cAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTAc                                                                                                            <  1:95867/150‑1 (MQ=255)
                                        cAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTAc                                                                                                            >  2:932018/1‑150 (MQ=255)
                                        cAAAGGAACAGAAAAAAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTAc                                                                                                            <  2:938402/150‑1 (MQ=255)
                                         aaaGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCGCCTCTTTACt                                                                                                           >  2:1691892/1‑150 (MQ=255)
                                         aaaGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACt                                                                                                           <  1:353090/150‑1 (MQ=255)
                                          aaGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACtt                                                                                                          <  1:1264396/150‑1 (MQ=255)
                                           atgAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCCCCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCAACTCTTTACTTa                                                                                                         >  2:880050/3‑150 (MQ=255)
                                              aaCAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTaaaa                                                                                                      <  2:1225069/150‑1 (MQ=255)
                                              aaCAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTaaaa                                                                                                      <  1:798778/150‑1 (MQ=255)
                                               aCAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTaaaaa                                                                                                     <  1:1447586/150‑1 (MQ=255)
                                               aCAGAAAATAAATCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCATCCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTaaaaa                                                                                                     <  2:784062/150‑1 (MQ=255)
                                               aCAGAAAAAAAATCCCAGCTCCGCCAGTTTTCCGCGAGAAGTAGTGGCAAACGAAACTGAGTCGAAAATGGCATCCACCGCAACCTCTTTGCCTTCCCGCACGGGAACGGCCACCTGCTCAAACGCCGACTCCACCTCTTTACTTaaaaa                                                                                                     <  1:1159558/150‑1 (MQ=255)
                                                 aGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAgg                                                                                                   <  1:882562/150‑1 (MQ=255)
                                                 aGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAgg                                                                                                   <  1:747826/150‑1 (MQ=255)
                                                           tCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGcc                                                                                         >  2:20193/1‑150 (MQ=255)
                                                              cTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGt                                                                                      <  2:656379/150‑1 (MQ=255)
                                                              cTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGt                                                                                      >  1:656379/1‑150 (MQ=255)
                                                               tGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGtt                                                                                     <  2:1265700/150‑1 (MQ=255)
                                                                   ccGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAACATGGCATCCACCGCCACCTCTATGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTtgct                                                                                 <  1:265570/150‑1 (MQ=255)
                                                                   ccGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTtgct                                                                                 <  1:1358543/150‑1 (MQ=255)
                                                                     gCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGACCGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTtgctgc                                                                               >  2:1079381/1‑150 (MQ=255)
                                                                     gCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTtgctgc                                                                               >  1:490092/1‑150 (MQ=255)
                                                                     gCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACAGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCAAGTTtgctgc                                                                               >  2:565734/1‑150 (MQ=255)
                                                                      ccAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCa                                                                              >  2:695743/1‑150 (MQ=255)
                                                                        aGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCAcc                                                                            <  2:361328/150‑1 (MQ=255)
                                                                        aGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCAcc                                                                            >  2:177803/1‑150 (MQ=255)
                                                                              tCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCa                                                                      <  2:271226/150‑1 (MQ=255)
                                                                                   ttatgtagtGGCAATCGATACTGAGTCGATAATGGCATTCACCGCCACCTCTTTGCGTTCCCGCACGGGATCGCCCAACAGCTCAAACGCCGCCTGGAGGTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCAGCGCGGCAGGTTc                                                                 <  2:1096817/146‑1 (MQ=255)
                                                                                     aaGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAg                                                               <  2:1341681/150‑1 (MQ=255)
                                                                                       gtagtGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGAc                                                             <  2:486766/150‑1 (MQ=255)
                                                                                         agtGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGAcgc                                                           >  2:174274/1‑150 (MQ=255)
                                                                                          gtGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCGCGGGAACGCCCAACTGCTCAAACGCCGCCTCCCCCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGGGCAGAcgcg                                                          >  2:329378/1‑150 (MQ=255)
                                                                                          gtGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGAcgcg                                                          <  1:996823/150‑1 (MQ=255)
                                                                                          gtGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGAcgcg                                                          <  1:932018/150‑1 (MQ=255)
                                                                                          gtGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACATAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGCcgcc                                                          >  2:1437799/1‑149 (MQ=255)
                                                                                           tGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACgcgc                                                         >  1:727122/1‑150 (MQ=255)
                                                                                             gCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCCCCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGGTCAGACGCGCaa                                                       >  2:1246714/1‑150 (MQ=255)
                                                                                              cAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAg                                                      >  2:968337/1‑150 (MQ=255)
                                                                                               aaCCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAgt                                                     >  2:1042076/1‑150 (MQ=255)
                                                                                                aCCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCAATTTTGCTGCACCGCGCCAGGTTCCGACGAGCAAgtg                                                    >  1:872368/1‑150 (MQ=255)
                                                                                                 ccGAAACTGAGTCGAAAATTGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAgtgt                                                   <  2:412006/150‑1 (MQ=255)
                                                                                                   gAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTgtcg                                                 >  1:950075/1‑150 (MQ=255)
                                                                                                   gAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTgtcg                                                 <  2:622242/150‑1 (MQ=255)
                                                                                                     aaCTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTgtcgtc                                               <  1:1197959/150‑1 (MQ=255)
                                                                                                             cGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTAc                                       <  1:177803/150‑1 (MQ=255)
                                                                                                              gAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTAcc                                      <  1:547835/150‑1 (MQ=255)
                                                                                                                    ggCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTAGCTGCCCCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAAATACCGCACGa                                >  2:1377633/1‑150 (MQ=255)
                                                                                                                     gCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGAt                               <  1:1618240/150‑1 (MQ=255)
                                                                                                                     gCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGGTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGAt                               >  1:210296/1‑150 (MQ=255)
                                                                                                                      cATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCGCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGATg                              >  2:957104/1‑150 (MQ=255)
                                                                                                                        tCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGATGGt                            >  2:434725/1‑150 (MQ=255)
                                                                                                                          cACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCGCCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGATGGTGc                          <  1:1691892/150‑1 (MQ=255)
                                                                                                                          cACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGATGGTGc                          <  1:1241537/150‑1 (MQ=255)
                                                                                                                          cACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGATGGTGc                          >  2:926636/1‑150 (MQ=255)
                                                                                                                             cGCCACCTCTTTGCCTTCGCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGATGGTGCTGa                       <  2:950075/150‑1 (MQ=255)
                                                                                                                              gCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCACCTCCACCTCTTTACTTAAAAACGCGTTCGCGCGAGTTTGCTGCAGCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGATGGTGCTGag                      <  2:816454/150‑1 (MQ=255)
                                                                                                                               ccACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCACACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGATGGTGCTGagt                     <  1:556310/150‑1 (MQ=255)
                                                                                                                               ccACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGATGGTGCTGagt                     <  1:563113/150‑1 (MQ=255)
                                                                                                                                  ccTCTTTGCCTTCCCTCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCTCCAGTTTGCTGCACCGGGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGATGGTGCTGagtagt                  <  1:1587166/150‑1 (MQ=255)
                                                                                                                                       ttGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGATGGTGCTGAGTAGTAGCTg             <  2:727122/150‑1 (MQ=255)
                                                                                                                                         gcgttCCCGTGCGGGAAGGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGATGGTGCTGAGTAGTAGCTGTa           >  1:615920/4‑150 (MQ=255)
                                                                                                                                          ccTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGATGCACCGCGCCAGGTGCAGCGGCGCAAGTGTCGTCACAATTACCGCACGATGGTGCTGAGTGGTAGCTGTaa          >  1:828742/1‑150 (MQ=255)
                                                                                                                                            ttCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGATGGTGCTGAGTAGTAGCTGTAATc        >  1:715224/1‑150 (MQ=255)
                                                                                                                                              cccGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGATGGTGCTGAGTAGTAGCTGTAATCCt      <  2:1619370/150‑1 (MQ=255)
                                                                                                                                                  aaCGGGAACGCCCAACTGCTCAATCGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGATGGTGCTGAGTAGTAGCTGTAATCCTGATa  <  2:631490/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTTGCCTTCCCGCACGGGAACGCCCAACTGCTCAAACGCCGCCTCCACCTCTTTACTTAAAAAGGCGTTCGCGCCAGTTTGCTGCACCGCGCCAGGTTCAGACGCGCAAGTGTCGTCACAATTACCGCACGATGGTGCTGAGTAGTAGCTGTAATCCTGATA  >  NC_000913/1763462‑1763757

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