Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,326,970 C→A 5.6% H556N (CAC→AAC)  topA → DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,326,9700CA 5.6% 940.9 / 27.3 248H556N (CAC→AAC) topADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base C (125/109);  new base A (10/4);  total (135/113)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 253 total aligned reads displayed.

ATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAA  >  CP009273/1326821‑1327119
                                                                                                                                                     |                                                                                                                                                     
aTCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATc                                                                                                                                                       >  1:2450546/1‑150 (MQ=255)
    tGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGa                                                                                                                                                   >  1:781358/1‑150 (MQ=255)
      gCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAg                                                                                                                                                 >  2:5164455/1‑150 (MQ=255)
       cTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTACGCGAGTTAATGAACTACGACTTTACCGCGCCGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGc                                                                                                                                                >  1:4875731/1‑150 (MQ=255)
        taagTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACAGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCa                                                                                                                                               <  2:5041774/147‑1 (MQ=255)
           gTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCagag                                                                                                                                            <  1:2268733/150‑1 (MQ=255)
           gTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCagag                                                                                                                                            <  2:208685/150‑1 (MQ=255)
            tGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGt                                                                                                                                           <  2:565016/150‑1 (MQ=255)
              cGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGAgg                                                                                                                                         >  1:4868153/1‑150 (MQ=255)
               gtgTAGAAAATCGTCGTTTATATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAATACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGa                                                                                                                                        <  1:4049755/148‑1 (MQ=255)
                   aGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGc                                                                                                                                    >  2:3024075/1‑150 (MQ=255)
                   aGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGc                                                                                                                                    >  1:4067642/1‑150 (MQ=255)
                     aaaaTCGTCGTTTCTCTGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTg                                                                                                                                  <  2:5185860/150‑1 (MQ=255)
                     aaaaTCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTg                                                                                                                                  <  1:4736957/150‑1 (MQ=255)
                        aTCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTAc                                                                                                                               <  1:1537611/150‑1 (MQ=255)
                        aTCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTAc                                                                                                                               <  2:2764655/150‑1 (MQ=255)
                         tcgtcgTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACt                                                                                                                              <  2:3331465/150‑1 (MQ=255)
                         tcgtcgTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACt                                                                                                                              <  1:371187/150‑1 (MQ=255)
                         tcgtcgTTTATATGGGGAATAATTGGGCGAACTCTTCACGGATCGCCTTGAAGAAAATTTCCGCGAGTGAATGATATTCGACTTTACCGCGCTGATGGAATACAGCCTCGACCAGGTGGCAAATAACGAAGCAGATTGGTAAGGTGTACt                                                                                                                              <  2:2403897/150‑1 (MQ=255)
                               tttCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAGAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCa                                                                                                                        >  1:4966467/1‑150 (MQ=255)
                                 tCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCAct                                                                                                                      <  1:3472930/150‑1 (MQ=255)
                                   tATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCActtc                                                                                                                    >  1:1561226/1‑150 (MQ=255)
                                    aTGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCActtct                                                                                                                   >  2:2150497/1‑150 (MQ=255)
                                        ggAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTtctc                                                                                                               <  1:4930399/150‑1 (MQ=255)
                                        ggAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTtctc                                                                                                               <  2:3723759/150‑1 (MQ=255)
                                         gAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCg                                                                                                              <  2:1082504/150‑1 (MQ=255)
                                           aaaaaTGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGa                                                                                                            >  1:4672856/1‑150 (MQ=255)
                                            aaaaTGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATAACTTCTTATCGGAt                                                                                                           >  1:4322311/1‑150 (MQ=255)
                                              aaTGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGAttt                                                                                                         >  2:1433131/1‑150 (MQ=255)
                                                    cGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCca                                                                                                   >  1:3225022/1‑150 (MQ=255)
                                                     gAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCcag                                                                                                  >  1:3187493/1‑150 (MQ=255)
                                                       aaTCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCcagca                                                                                                >  2:872855/1‑150 (MQ=255)
                                                       aaTCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCcagca                                                                                                >  2:2967712/1‑150 (MQ=255)
                                                        aTCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCcagcag                                                                                               <  1:3906035/150‑1 (MQ=255)
                                                          cGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGtt                                                                                             >  1:3024145/1‑150 (MQ=255)
                                                           gTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTa                                                                                            <  1:3982938/150‑1 (MQ=255)
                                                           gTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTa                                                                                            >  1:1337554/1‑150 (MQ=255)
                                                            tCACCGATCGCCTAGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAg                                                                                           >  2:718522/1‑150 (MQ=255)
                                                             cACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGa                                                                                          >  1:4580392/1‑150 (MQ=255)
                                                              aCCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGAt                                                                                         >  1:1243879/1‑150 (MQ=255)
                                                               ccGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATa                                                                                        <  1:4970762/150‑1 (MQ=255)
                                                                cGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATaa                                                                                       >  1:1041061/1‑150 (MQ=255)
                                                                 gATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATaaa                                                                                      >  2:4395711/1‑150 (MQ=255)
                                                                 gATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATAACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATaaa                                                                                      >  1:5191483/1‑150 (MQ=255)
                                                                     gCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATAACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTg                                                                                  <  2:4615861/150‑1 (MQ=255)
                                                                     gCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATAACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTg                                                                                  <  2:2011052/150‑1 (MQ=255)
                                                                      ccTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGa                                                                                 <  2:2368886/150‑1 (MQ=255)
                                                                       cTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGan                                                                                <  2:1579512/150‑2 (MQ=255)
                                                                          gaagaaAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGaaaaa                                                                             <  2:2541071/150‑1 (MQ=255)
                                                                           aagaaAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAg                                                                            <  1:607461/150‑1 (MQ=255)
                                                                           aagaaAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAg                                                                            <  1:2909800/150‑1 (MQ=255)
                                                                           aagaaAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAg                                                                            <  2:1611764/150‑1 (MQ=255)
                                                                           aagaaAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAg                                                                            >  1:1067257/1‑150 (MQ=255)
                                                                             gaaAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGAt                                                                          <  1:3004899/150‑1 (MQ=255)
                                                                             gaaAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGAt                                                                          >  2:4762984/1‑150 (MQ=255)
                                                                              aaaaTTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATAACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATc                                                                         >  1:3436174/1‑150 (MQ=255)
                                                                               aaaTTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATcc                                                                        <  2:4520446/150‑1 (MQ=255)
                                                                               aaaTTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATcc                                                                        <  2:1914569/150‑1 (MQ=255)
                                                                                aaTTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCg                                                                       >  1:2939960/1‑150 (MQ=255)
                                                                                   ttCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGaa                                                                    >  2:3431771/1‑150 (MQ=255)
                                                                                      cgcgAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAagag                                                                 >  2:413716/1‑150 (MQ=255)
                                                                                       gcgAGTTTATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAgg                                                                <  2:4850609/150‑1 (MQ=255)
                                                                                         gAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGggt                                                              >  2:3327325/1‑150 (MQ=255)
                                                                                          aGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGggtg                                                             <  2:1736598/150‑1 (MQ=255)
                                                                                           gTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGggtgg                                                            <  2:1201307/150‑1 (MQ=255)
                                                                                           gTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGggtgg                                                            <  1:1927304/150‑1 (MQ=255)
                                                                                             tAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGACGAGGGTGGTa                                                          >  2:2686770/1‑150 (MQ=255)
                                                                                             tAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATAACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTa                                                          >  1:1498838/1‑150 (MQ=255)
                                                                                                tGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATgc                                                       <  2:1881420/150‑1 (MQ=255)
                                                                                                tGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATgc                                                       >  2:3020915/1‑150 (MQ=255)
                                                                                                         aCTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAAcc                                              >  2:57942/1‑150 (MQ=255)
                                                                                                         aCTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAAcc                                              <  2:301941/150‑1 (MQ=255)
                                                                                                         aCTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATAACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCAGAAcc                                              >  1:279907/1‑150 (MQ=255)
                                                                                                          cTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCa                                             <  1:4818877/150‑1 (MQ=255)
                                                                                                            ttACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGa                                           <  1:5126705/150‑1 (MQ=255)
                                                                                                                 gcgcAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGtt                                      <  2:769646/150‑1 (MQ=255)
                                                                                                                  cgcAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTc                                     >  2:3445434/1‑150 (MQ=255)
                                                                                                                  cgcAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTc                                     >  1:3616305/1‑150 (MQ=255)
                                                                                                                     aGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGa                                  <  2:2897103/150‑1 (MQ=255)
                                                                                                                           aaaaCAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCa                            <  2:4723080/150‑1 (MQ=255)
                                                                                                                           aaaaCAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCa                            >  1:4985337/1‑150 (MQ=255)
                                                                                                                            aaaCAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACAAGCAt                           >  2:3625291/1‑150 (MQ=255)
                                                                                                                             aaCAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCAtt                          >  1:1949066/1‑150 (MQ=255)
                                                                                                                             aaCAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCAtt                          >  1:4838110/1‑150 (MQ=255)
                                                                                                                             aaCAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCAtt                          <  2:2513410/150‑1 (MQ=255)
                                                                                                                             aaCAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCAtt                          >  2:1769221/1‑150 (MQ=255)
                                                                                                                             aaCAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTACTGACCAGCAtt                          >  2:1993249/1‑150 (MQ=255)
                                                                                                                              aCAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTg                         >  1:4187073/1‑150 (MQ=255)
                                                                                                                                aGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGAc                       >  2:3258186/1‑150 (MQ=255)
                                                                                                                                         cAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACt              <  2:2773348/150‑1 (MQ=255)
                                                                                                                                         cAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACt              >  2:1117763/1‑150 (MQ=255)
                                                                                                                                           ggTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTtg            >  1:5144238/1‑150 (MQ=255)
                                                                                                                                           ggTGCCACAACACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTtg            <  1:2752095/150‑1 (MQ=255)
                                                                                                                                            gtgtCAAATCACGAAGCATAGTGGAAAGCTGTACTGGATCACTTGTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGACCGACTtgt           <  1:1709584/150‑1 (MQ=255)
                                                                                                                                                cAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTc       <  2:4649873/150‑1 (MQ=255)
                                                                                                                                                cAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTc       <  1:1433168/150‑1 (MQ=255)
                                                                                                                                                cAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTc       >  1:2039925/1‑150 (MQ=255)
                                                                                                                                                    tCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCaa   >  1:2221843/1‑150 (MQ=255)
                                                                                                                                                     cACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCaaa  >  1:2099859/1‑150 (MQ=255)
                                                                                                                                                     cACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCaaa  >  1:2099887/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTCACCGATCGCCTTGAAGAAAATTTCCGCGAGTTAATGAACTACGACTTTACCGCGCAGATGGAAAACAGCCTCGACCAGGTGGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAA  >  CP009273/1326821‑1327119

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