Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,403,534 C→T 100% intergenic (+224/‑105) prmA → / → dusB methyltransferase for 50S ribosomal subunit protein L11/tRNA‑dihydrouridine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,403,5340CT96.6% 759.4 / 12.5 233intergenic (+224/‑105)prmA/dusBmethyltransferase for 50S ribosomal subunit protein L11/tRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base C (6/2);  new base T (125/100);  total (131/102)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 233 total aligned reads displayed.

CAAATTTGCGTGCCAAAATTTTCATTCATAAAGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCG  >  CP009273/3403391‑3403678
                                                                                                                                               |                                                                                                                                                
cAAATTTGCGTGCCAAAATTTTCATTCATAAAGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTg                                                                                                                                              >  2:2445601/1‑148 (MQ=255)
cAAATTTGCGTGCCAAAATTTTCATTCATAAAGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTg                                                                                                                                              <  1:2732759/148‑1 (MQ=255)
  aaTTTGCGTGCCAAAATTTTCATTCATAAAGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGAt                                                                                                                                            <  2:3068627/148‑1 (MQ=255)
  aaTTTGCGTGCCAAAATTTTCATTCATAAAGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGAt                                                                                                                                            <  2:125398/148‑1 (MQ=255)
  aaTTTGCGTGCCAAAATTTTCATTCATAAAGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGAt                                                                                                                                            <  2:1440679/148‑1 (MQ=255)
  aaTGTGCGTGCCAAAATTTTCATTCATAAAGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGAt                                                                                                                                            <  2:3352730/148‑1 (MQ=255)
      tGCGTGCCAAAATTTTCATTCATAAAGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCAt                                                                                                                                        <  2:218283/148‑1 (MQ=255)
                 aTTTTCATTCATAAAGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGAtt                                                                                                                             >  2:1747603/1‑148 (MQ=255)
                  ttttCATTCATAAAGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTg                                                                                                                            >  1:447204/1‑148 (MQ=255)
                      cattcatAAAGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCa                                                                                                                        >  1:3510325/1‑148 (MQ=255)
                         tcatAAAGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAg                                                                                                                     >  2:1272135/1‑148 (MQ=255)
                          catAAAGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGt                                                                                                                    <  1:2369714/148‑1 (MQ=255)
                             aaaGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTg                                                                                                                 >  2:3074552/1‑148 (MQ=255)
                             aaaGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTg                                                                                                                 >  1:1182107/1‑148 (MQ=255)
                             aaaGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTg                                                                                                                 >  2:3526875/1‑148 (MQ=255)
                             aaaGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCAGTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCACAGTTTg                                                                                                                 >  2:616114/1‑148 (MQ=255)
                              aaGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTgg                                                                                                                <  2:517097/148‑1 (MQ=255)
                               aGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGc                                                                                                               >  1:1541284/1‑148 (MQ=255)
                               aGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGc                                                                                                               <  1:1291340/148‑1 (MQ=255)
                                gAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGcc                                                                                                              <  2:3473691/148‑1 (MQ=255)
                                gAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGcc                                                                                                              <  2:1714900/148‑1 (MQ=255)
                                   aaaTTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCttt                                                                                                           <  2:2002608/148‑1 (MQ=255)
                                   aaaTTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCttt                                                                                                           >  2:1902152/1‑148 (MQ=255)
                                   aaaTTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCttt                                                                                                           <  2:3140518/148‑1 (MQ=255)
                                     aTTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCa                                                                                                         <  1:2899798/148‑1 (MQ=255)
                                       tGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCAtc                                                                                                       >  1:3239225/1‑148 (MQ=255)
                                          gaACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCg                                                                                                    <  1:463925/148‑1 (MQ=255)
                                          gaACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCg                                                                                                    >  1:1682026/1‑148 (MQ=255)
                                          gaACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCg                                                                                                    <  1:2003266/148‑1 (MQ=255)
                                                    aaaTTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAAt                                                                                          <  2:504778/148‑1 (MQ=255)
                                                      aTTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGc                                                                                        >  1:2349823/1‑148 (MQ=255)
                                                             tGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAatat                                                                                 <  2:250153/148‑1 (MQ=255)
                                                             tGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAatat                                                                                 <  2:617710/148‑1 (MQ=255)
                                                                gtgtAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATAcg                                                                              <  2:3179096/148‑1 (MQ=255)
                                                                  gtAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATAcgcc                                                                            <  2:2311978/148‑1 (MQ=255)
                                                                    taaTTTTAGTCACTATTTTCTACTATGATGATTTTTATGATTAATTATCGCACGACGCTCATTTTAAATGCTATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATAcgccgc                                                                          <  2:1550846/147‑1 (MQ=255)
                                                                    aaaTTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATAcgccgc                                                                          <  2:3303167/148‑1 (MQ=255)
                                                                    aaaTTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATAcgccgc                                                                          >  1:762040/1‑148 (MQ=255)
                                                                    aaaTTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATAcgccgc                                                                          <  1:3526875/148‑1 (MQ=255)
                                                                    aaaTTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATAcgccgc                                                                          <  1:1568177/148‑1 (MQ=255)
                                                                    aaaTTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATAcgccgc                                                                          <  1:3074552/148‑1 (MQ=255)
                                                                    aaaTTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATAcgccgc                                                                          <  1:3073902/148‑1 (MQ=255)
                                                                      aTTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCt                                                                        <  2:1395650/148‑1 (MQ=255)
                                                                       ttttAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCtt                                                                       >  2:812720/1‑148 (MQ=255)
                                                                        tttAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTg                                                                      >  1:2937567/1‑148 (MQ=255)
                                                                         ttAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGc                                                                     <  1:1345262/148‑1 (MQ=255)
                                                                          tAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCa                                                                    <  2:236057/148‑1 (MQ=255)
                                                                             tCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTc                                                                 <  2:1649406/148‑1 (MQ=255)
                                                                              cACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTca                                                                >  2:49264/1‑148 (MQ=255)
                                                                              cACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTca                                                                >  1:1868703/1‑148 (MQ=255)
                                                                               aCTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCTGCCTTGCAGTcac                                                               <  1:1747603/148‑1 (MQ=255)
                                                                               aCTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTcac                                                               <  1:2614777/148‑1 (MQ=255)
                                                                               aCTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTcac                                                               <  1:3167952/148‑1 (MQ=255)
                                                                                cTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACTCCGCCTTGCAGTcaca                                                              >  1:2978880/1‑148 (MQ=255)
                                                                                cTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTcaca                                                              >  1:45303/1‑148 (MQ=255)
                                                                                cTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTcaca                                                              >  1:1202943/1‑148 (MQ=255)
                                                                                cTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTcaca                                                              <  1:1447219/148‑1 (MQ=255)
                                                                                  aTTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGt                                                            <  2:1488603/148‑1 (MQ=255)
                                                                                   ttttCTCATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTa                                                           <  2:2073744/148‑1 (MQ=255)
                                                                                     ttCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATg                                                         <  2:2355510/148‑1 (MQ=255)
                                                                                         aaTATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCa                                                     >  2:901120/1‑148 (MQ=255)
                                                                                          atatGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCAt                                                    >  1:1033964/1‑148 (MQ=255)
                                                                                          atatGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCAt                                                    >  2:2261662/1‑148 (MQ=255)
                                                                                           tatGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCAtt                                                   <  1:621126/148‑1 (MQ=255)
                                                                                           tatGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCAtt                                                   >  1:2523110/1‑148 (MQ=255)
                                                                                           tatGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCAtt                                                   >  2:1110851/1‑148 (MQ=255)
                                                                                            atgatgATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCAttt                                                  <  1:1151934/148‑1 (MQ=255)
                                                                                              gatgatTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCt                                                >  1:182732/1‑148 (MQ=255)
                                                                                              gatgatTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCt                                                >  1:879707/1‑148 (MQ=255)
                                                                                               atgatTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCtt                                               >  2:48431/1‑148 (MQ=255)
                                                                                                tgatTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTa                                              <  1:2489827/148‑1 (MQ=255)
                                                                                                   tttttATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACt                                           <  1:912938/148‑1 (MQ=255)
                                                                                                    ttttATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTc                                          <  2:381578/148‑1 (MQ=255)
                                                                                                    ttttATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTc                                          >  2:121651/1‑148 (MQ=255)
                                                                                                         tGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATgcg                                     <  2:762040/148‑1 (MQ=255)
                                                                                                          gAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATgcgc                                    <  2:3239225/148‑1 (MQ=255)
                                                                                                          gAGTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATgcgc                                    >  2:2917806/1‑148 (MQ=255)
                                                                                                          gAGTAATTATCGCACCACGCTCATTTTAAATGCAATGTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCNTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATgcgc                                    >  1:1802319/1‑148 (MQ=255)
                                                                                                            gTAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCAt                                  <  1:34821/148‑1 (MQ=255)
                                                                                                             tAATTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATc                                 >  1:2398345/1‑148 (MQ=255)
                                                                                                              aaTTATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCg                                >  1:275438/1‑148 (MQ=255)
                                                                                                              aaTTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCg                                >  2:2296715/1‑148 (MQ=255)
                                                                                                                ttATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGa                              >  2:2855632/1‑148 (MQ=255)
                                                                                                                 tATCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGAc                             >  2:2366230/1‑148 (MQ=255)
                                                                                                                  aTCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACa                            >  1:3250959/1‑148 (MQ=255)
                                                                                                                  aTCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACa                            >  1:1288485/1‑148 (MQ=255)
                                                                                                                   tCGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACaa                           >  2:1741556/1‑148 (MQ=255)
                                                                                                                    cGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAat                          >  1:279719/1‑148 (MQ=255)
                                                                                                                    cGCACCACGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAat                          <  2:249292/148‑1 (MQ=255)
                                                                                                                        ccacGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCa                      >  2:1896972/1‑148 (MQ=255)
                                                                                                                         cacGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAg                     >  2:3032467/1‑148 (MQ=255)
                                                                                                                         cacGCTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAg                     >  2:2649874/1‑148 (MQ=255)
                                                                                                                             cTCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCa                 <  2:1182107/148‑1 (MQ=255)
                                                                                                                              tCATTTTAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAg                >  2:2147595/1‑148 (MQ=255)
                                                                                                                                 ttttAAATGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGaaa             >  1:3083286/1‑148 (MQ=255)
                                                                                                                                     aaaTGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGc         <  1:553222/148‑1 (MQ=255)
                                                                                                                                      aaTGCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGcc        <  2:1141992/148‑1 (MQ=255)
                                                                                                                                         gCAATTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGa     >  2:1301201/1‑148 (MQ=255)
                                                                                                                                           aaTTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATc   >  2:144776/1‑148 (MQ=255)
                                                                                                                                            aTTTTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCg  >  2:3007400/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CAAATTTGCGTGCCAAAATTTTCATTCATAAAGAAAAATTGAGAACTTACTCAAATTTCTTTGAGTGTAAATTTTAGTCACTATTTTCTAATATGATGATTTTTATGAGTAATTATCGCACCACGCTCATTTTAAATGCAATTCTTTGATCCATCTCAGAGGATTGGTCAAAGTTTGGCCTTTCATCTCGTGCAAAAAATGCGTAATATACGCCGCCTTGCAGTCACAGTATGGTCATTTCTTAACTCATGCGCATCGGACAATATCAGCTCAGAAATCGCCTGATCG  >  CP009273/3403391‑3403678

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: