Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,079,951 G→T 6.9% D374Y (GAT→TAT)  efeB → deferrrochelatase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,079,9510GT 6.9% 1106.8 / 34.2 293D374Y (GAT→TAT) efeBdeferrrochelatase, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (148/124);  new base T (15/6);  total (163/130)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 298 total aligned reads displayed.

CGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCGGGT  >  CP009273/1079805‑1080099
                                                                                                                                                  |                                                                                                                                                    
   tGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTg                                                                                                                                                 >  1:2730852/1‑149 (MQ=255)
   tGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTg                                                                                                                                                 <  1:4102077/149‑1 (MQ=255)
   tGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTg                                                                                                                                                 <  1:1773809/149‑1 (MQ=255)
   tGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGAGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTg                                                                                                                                                 <  1:1711817/149‑1 (MQ=255)
   tGGAAAGCCAAATCCGGCTGGCGAATCCCCGCACGGCGGAGAGAGAGTCCATCCTGATGCTGCGTCGTGGCTACAGTTAATCACTGGGCGTCACCAACTCCGGGCAACTCGCTATGGGGATGCTGTTTGTCAGCTACCAACACGATCtc                                                                                                                                                 >  2:1584790/1‑148 (MQ=255)
     gACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGa                                                                                                                                               >  1:4291436/1‑149 (MQ=255)
       cAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGaaa                                                                                                                                             >  1:623952/1‑149 (MQ=255)
         gCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGaaaaa                                                                                                                                           <  1:79075/149‑1 (MQ=255)
            atatCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGc                                                                                                                                        >  2:668135/1‑149 (MQ=255)
              atCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCtt                                                                                                                                      >  1:116235/1‑149 (MQ=255)
                ccGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTcc                                                                                                                                    <  2:3960363/149‑1 (MQ=255)
                   gCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGa                                                                                                                                 >  2:2140603/1‑149 (MQ=255)
                     tGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACa                                                                                                                               >  2:953971/1‑149 (MQ=255)
                     tGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACa                                                                                                                               >  1:3130969/1‑149 (MQ=255)
                           aTCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACaa                                                                                                                         >  2:2031061/1‑149 (MQ=255)
                            tCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACaaa                                                                                                                        >  1:4586367/1‑149 (MQ=255)
                             ccccGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACaaaa                                                                                                                       >  2:306655/1‑149 (MQ=255)
                                   aCGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTAGTCACTGGGCGTCACCAACTCCTGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCt                                                                                                                 <  2:3576171/149‑1 (MQ=255)
                                    cggcggAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTc                                                                                                                <  2:2863851/149‑1 (MQ=255)
                                       cggAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAAt                                                                                                             <  2:1434681/149‑1 (MQ=255)
                                           gagtgagtCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCg                                                                                                         <  2:1682041/149‑1 (MQ=255)
                                             gtgagtCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGaa                                                                                                       <  1:4409800/149‑1 (MQ=255)
                                             gtgagtCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGaa                                                                                                       >  1:1021817/1‑149 (MQ=255)
                                                 gtCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAATGCTCAATGGCGAAgcgc                                                                                                   >  2:2032686/1‑149 (MQ=255)
                                                 gtCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAgcgc                                                                                                   >  1:4560136/1‑149 (MQ=255)
                                                     aGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTgga                                                                                               >  1:3748900/1‑149 (MQ=255)
                                                     aGCCTGATGCTGCGTCGTGGCTACAGTTAGTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTgga                                                                                               <  2:247912/149‑1 (MQ=255)
                                                          gATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAAt                                                                                          <  2:2435636/149‑1 (MQ=255)
                                                             gctgcGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACg                                                                                       >  2:2805558/1‑149 (MQ=255)
                                                              ctgcGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGt                                                                                      <  2:4677753/149‑1 (MQ=255)
                                                                 cgtcgtGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTaa                                                                                   <  2:1111790/149‑1 (MQ=255)
                                                                   tcgtGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAAc                                                                                 >  2:1228474/1‑149 (MQ=255)
                                                                   tcgtGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAAc                                                                                 >  2:692952/1‑149 (MQ=255)
                                                                     gtGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCt                                                                               >  1:988620/1‑149 (MQ=255)
                                                                     gtGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCt                                                                               <  2:2314241/149‑1 (MQ=255)
                                                                       ggCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTAt                                                                             >  2:4290183/1‑149 (MQ=255)
                                                                        gCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCACCTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATc                                                                            >  1:1385737/1‑149 (MQ=255)
                                                                        gCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATc                                                                            >  1:1386196/1‑149 (MQ=255)
                                                                            cAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATcggcg                                                                        >  1:2425618/1‑149 (MQ=255)
                                                                             aGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATcggcgg                                                                       >  2:1955080/1‑149 (MQ=255)
                                                                                tATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATcggcggcgg                                                                    >  2:489197/1‑149 (MQ=255)
                                                                                  ttCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGtt                                                                  >  1:1061088/1‑149 (MQ=255)
                                                                                     aCTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAtt                                                               <  1:3413714/149‑1 (MQ=255)
                                                                                      cTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAttt                                                              <  2:860634/149‑1 (MQ=255)
                                                                                         ggCGTCACCATCTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAtttttt                                                           <  2:1016736/149‑1 (MQ=255)
                                                                                         ggCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAtttttt                                                           >  2:3350524/1‑149 (MQ=255)
                                                                                         ggCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAtttttt                                                           >  1:1717426/1‑149 (MQ=255)
                                                                                            gTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTgc                                                        <  1:2547339/149‑1 (MQ=255)
                                                                                               aCCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAGTACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCt                                                     >  1:350756/1‑149 (MQ=255)
                                                                                                ccAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTg                                                    >  1:858656/1‑149 (MQ=255)
                                                                                                 cAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGTTGTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGc                                                   <  2:397680/149‑1 (MQ=255)
                                                                                                 cAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGc                                                   <  2:2730852/149‑1 (MQ=255)
                                                                                                 cAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGc                                                   <  2:2928834/149‑1 (MQ=255)
                                                                                                  aaCTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTTCGCTGcc                                                  >  1:3197719/1‑149 (MQ=255)
                                                                                                   aCTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTCTTTTTTTGCGCTGCCg                                                 >  2:4694018/1‑149 (MQ=255)
                                                                                                     tCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCggg                                               >  1:329058/1‑149 (MQ=255)
                                                                                                      ccGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCgggg                                              >  2:4089805/1‑149 (MQ=255)
                                                                                                      ccGGGCAACTGGAGATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCgggg                                              >  2:547751/1‑149 (MQ=255)
                                                                                                        gggCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGt                                            <  1:2023559/149‑1 (MQ=255)
                                                                                                        gggCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGt                                            <  1:295600/149‑1 (MQ=255)
                                                                                                        gggCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGt                                            >  2:3766244/1‑149 (MQ=255)
                                                                                                         ggCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTg                                           <  2:3136260/149‑1 (MQ=255)
                                                                                                             aCTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAgg                                       >  2:1921338/1‑149 (MQ=255)
                                                                                                              cTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGa                                      >  2:739648/1‑149 (MQ=255)
                                                                                                              cTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGa                                      >  1:156102/1‑149 (MQ=255)
                                                                                                               tGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGAc                                     <  1:1235258/149‑1 (MQ=255)
                                                                                                                ggATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGAcg                                    >  2:2067235/1‑149 (MQ=255)
                                                                                                                ggATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGAcg                                    >  2:1613802/1‑149 (MQ=255)
                                                                                                                 gATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGAcgc                                   >  1:4264497/1‑149 (MQ=255)
                                                                                                                 gATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGAcgc                                   >  2:1658136/1‑149 (MQ=255)
                                                                                                                     tGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAAc                               >  1:2573044/1‑149 (MQ=255)
                                                                                                                     tGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAAc                               <  2:2908337/149‑1 (MQ=255)
                                                                                                                     tGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAAc                               <  2:1826166/149‑1 (MQ=255)
                                                                                                                      ggggTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTTCGCTGCCCGGGGTGAAGGACGCCAACg                              >  1:3597071/1‑149 (MQ=255)
                                                                                                                      ggggTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACg                              <  1:2513966/149‑1 (MQ=255)
                                                                                                                          ttGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGatta                          >  2:1856886/1‑149 (MQ=255)
                                                                                                                          ttGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGatta                          <  1:1610876/149‑1 (MQ=255)
                                                                                                                            gCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGattatt                        >  2:4201265/1‑149 (MQ=255)
                                                                                                                            gCTGTTTGTCTGCTACCAACACGATCTGGAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTAttt                       <  1:4712677/149‑1 (MQ=255)
                                                                                                                              tGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTc                      >  2:742677/1‑149 (MQ=255)
                                                                                                                              tGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGACCGATTATTTc                      >  1:49333/1‑149 (MQ=255)
                                                                                                                               gTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCg                     >  2:4449988/1‑149 (MQ=255)
                                                                                                                               gTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCg                     >  1:2934139/1‑149 (MQ=255)
                                                                                                                               gTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCg                     >  1:1283559/1‑149 (MQ=255)
                                                                                                                                tttGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCgg                    >  2:824399/1‑149 (MQ=255)
                                                                                                                                 ttGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGa                   <  1:1845011/149‑1 (MQ=255)
                                                                                                                                 ttGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGa                   <  2:2246263/149‑1 (MQ=255)
                                                                                                                                  tGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGaa                  <  2:1959685/149‑1 (MQ=255)
                                                                                                                                      tGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAgcgc              >  2:3859436/1‑149 (MQ=255)
                                                                                                                                       gCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGcgcg             >  1:1953309/1‑149 (MQ=255)
                                                                                                                                        cTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGt            <  1:692952/149‑1 (MQ=255)
                                                                                                                                         tACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGtt           >  2:3696513/1‑149 (MQ=255)
                                                                                                                                          aCCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTa          >  2:174031/1‑149 (MQ=255)
                                                                                                                                          aCCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTa          >  1:3961066/1‑149 (MQ=255)
                                                                                                                                             aaCACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTg       >  1:2268483/1‑149 (MQ=255)
                                                                                                                                              acacGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGc      >  2:3230979/1‑149 (MQ=255)
                                                                                                                                               cacGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCg     >  2:4591153/1‑149 (MQ=255)
                                                                                                                                               cacGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCg     >  2:1312140/1‑149 (MQ=255)
                                                                                                                                                  gATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGTGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCGGGt  >  2:1265200/1‑149 (MQ=255)
                                                                                                                                                  gATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCGGGt  <  1:3812964/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCGGGT  >  CP009273/1079805‑1080099

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: