Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,421,791 G→T 22.4% R4S (CGC→AGC)  ytfE ← iron‑sulfur cluster repair protein RIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,421,7910GT22.4% 854.4 / 194.2 334R4S (CGC→AGC) ytfEiron‑sulfur cluster repair protein RIC
Reads supporting (aligned to +/- strand):  ref base G (109/150);  new base T (32/43);  total (141/193)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 340 total aligned reads displayed.

AGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTCAGTTAAAGGTTT  >  CP009273/4421642‑4421940
                                                                                                                                                     |                                                                                                                                                     
aGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTcgcg                                                                                                                                                       >  2:4240155/1‑150 (MQ=255)
aGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGAGCCAGTTCACCTAAAGGTTGGTcgcg                                                                                                                                                       <  2:4023377/150‑1 (MQ=255)
 gCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGa                                                                                                                                                      >  2:1699620/1‑150 (MQ=255)
  cTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCCTATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGAt                                                                                                                                                     <  2:2297204/150‑1 (MQ=255)
  cTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGAt                                                                                                                                                     <  1:4254434/150‑1 (MQ=255)
  cTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCAGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGACAGTTCACCTAAAGGTTGGTCGCGAt                                                                                                                                                     >  1:4790763/1‑150 (MQ=255)
   ttCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATa                                                                                                                                                    >  1:2923569/1‑150 (MQ=255)
   ttCAATGACCACAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCATACAGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGAGAGGAATAGAGAGCGCAAGTTCACCTAAAGGTTGGTCGCGATa                                                                                                                                                    >  2:2120605/1‑150 (MQ=255)
        tGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCa                                                                                                                                               <  1:960065/150‑1 (MQ=255)
         gACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCAt                                                                                                                                              <  2:1031809/150‑1 (MQ=255)
          aCCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATa                                                                                                                                             >  1:2757688/1‑150 (MQ=255)
          aCCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATa                                                                                                                                             <  1:4970583/150‑1 (MQ=255)
          aCCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATa                                                                                                                                             <  1:603490/150‑1 (MQ=255)
           ccTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAg                                                                                                                                            >  1:1886620/1‑150 (MQ=255)
            cTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGc                                                                                                                                           >  1:1373319/1‑150 (MQ=255)
             tCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCt                                                                                                                                          <  1:1436696/150‑1 (MQ=255)
             tCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCt                                                                                                                                          <  2:2923893/150‑1 (MQ=255)
             tCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCt                                                                                                                                          >  1:1023689/1‑150 (MQ=255)
              cAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTg                                                                                                                                         >  2:4141757/1‑150 (MQ=255)
                aCATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTACTCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGACTAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGAt                                                                                                                                       >  1:529523/1‑150 (MQ=255)
                      aGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCa                                                                                                                                 <  1:2696716/150‑1 (MQ=255)
                      aGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCa                                                                                                                                 >  1:147845/1‑150 (MQ=255)
                        ttCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCAtt                                                                                                                               >  2:4908163/1‑150 (MQ=255)
                         tCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTc                                                                                                                              <  1:4559133/150‑1 (MQ=255)
                         tCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTc                                                                                                                              <  1:369623/150‑1 (MQ=255)
                           tttACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCtt                                                                                                                            >  2:4754758/1‑150 (MQ=255)
                            ttACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTta                                                                                                                           <  2:327819/150‑1 (MQ=255)
                              aCGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTtaat                                                                                                                         <  1:3583750/150‑1 (MQ=255)
                              aCGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTtaat                                                                                                                         <  2:771314/150‑1 (MQ=255)
                                gTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTtaataa                                                                                                                       <  1:1118588/150‑1 (MQ=255)
                                  gccgccGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAataaga                                                                                                                     <  1:3834717/150‑1 (MQ=255)
                                   ccgccgCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAataagat                                                                                                                    <  2:592511/150‑1 (MQ=255)
                                     gccgcGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGACTAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAatcagataa                                                                                                                  >  2:5102060/1‑150 (MQ=255)
                                     gccgcGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAataagataa                                                                                                                  >  2:4668642/1‑150 (MQ=255)
                                       cgcgcgcgCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAAtaagataaga                                                                                                                <  1:4948504/150‑1 (MQ=255)
                                        gcgcgcgcCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATaagataagat                                                                                                               >  1:3916884/1‑150 (MQ=255)
                                         cgcgcgcgAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAATCAGAGCTGATGCGCGAGGAATCGAGAGCGCCAGTTCACCTTAAGGTTGGTCGCGATAAGACATAGCTGATACCTCATTCTTAATAAGATAAGATg                                                                                                              <  1:5081768/150‑1 (MQ=255)
                                         cgcgcgcCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATg                                                                                                              <  2:5053680/150‑1 (MQ=255)
                                           cgcgcCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTa                                                                                                            <  1:3814721/150‑1 (MQ=255)
                                            gcgcCAGCGTCTGCTTACCGCCACAGCAGGAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATCGAGAGCGCCAGTTCACCTACAGGTTGGTCGCGATAAGCCAGAGCTGATACCTCATTCTGAAtaagatcagagggat                                                                                                           >  1:4317682/1‑145 (MQ=255)
                                             cgcCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTAtt                                                                                                          <  2:73083/150‑1 (MQ=255)
                                              gcCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTAttt                                                                                                         >  1:4128192/1‑150 (MQ=255)
                                                 aGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTaa                                                                                                      >  2:1244366/1‑150 (MQ=255)
                                                  gCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTaaa                                                                                                     <  1:4844000/150‑1 (MQ=255)
                                                  gCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTaaa                                                                                                     <  2:2592429/150‑1 (MQ=255)
                                                  gCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTaaa                                                                                                     <  2:19465/150‑1 (MQ=255)
                                                   cGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAAt                                                                                                    <  2:2504730/150‑1 (MQ=255)
                                                     tCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGc                                                                                                  >  2:2989611/1‑150 (MQ=255)
                                                       tGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCAt                                                                                                <  1:2514173/150‑1 (MQ=255)
                                                          ttACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCtt                                                                                             <  1:2447765/150‑1 (MQ=255)
                                                             ccGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTaa                                                                                          <  1:224781/150‑1 (MQ=255)
                                                             ccGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTaa                                                                                          <  2:2226547/150‑1 (MQ=255)
                                                             ccGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTaa                                                                                          <  1:1545194/150‑1 (MQ=255)
                                                               gCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAgg                                                                                        <  2:941128/150‑1 (MQ=255)
                                                                 cacaGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATTCATCTTTAAGGCa                                                                                      >  1:2145556/1‑150 (MQ=255)
                                                                 cacaGCAGTAATCCATATCATATTTACGAAACAGAGCAGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCa                                                                                      <  2:3545506/150‑1 (MQ=255)
                                                                  acaGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCaa                                                                                     >  2:4587481/1‑150 (MQ=255)
                                                                     gcagTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAg                                                                                  >  1:1310779/1‑150 (MQ=255)
                                                                        gTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCta                                                                               <  1:2544787/150‑1 (MQ=255)
                                                                        gTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCta                                                                               <  2:2073179/150‑1 (MQ=255)
                                                                          aaTCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCtata                                                                             <  2:879544/150‑1 (MQ=255)
                                                                           aTCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGACAAAAGCTATaa                                                                            >  1:298049/1‑150 (MQ=255)
                                                                             ccATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAAcc                                                                          >  1:4110372/1‑150 (MQ=255)
                                                                              catatcatatTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAAccc                                                                         >  2:3737203/1‑150 (MQ=255)
                                                                                tatcatatTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCtt                                                                       <  2:1570822/150‑1 (MQ=255)
                                                                                   catatTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACg                                                                    >  1:1080871/1‑150 (MQ=255)
                                                                                    atatTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGt                                                                   <  1:2397273/150‑1 (MQ=255)
                                                                                       tttACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGc                                                                <  1:2456632/150‑1 (MQ=255)
                                                                                       tttACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGc                                                                <  1:2715021/150‑1 (MQ=255)
                                                                                          aCGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCaga                                                             >  1:557209/1‑150 (MQ=255)
                                                                                              aaCAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCt                                                         >  2:624494/1‑150 (MQ=255)
                                                                                                cAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTa                                                       >  2:1396066/1‑150 (MQ=255)
                                                                                                   agCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTAttt                                                    >  2:1975286/1‑150 (MQ=255)
                                                                                                    gCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTa                                                   >  2:27507/1‑150 (MQ=255)
                                                                                                    gCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTa                                                   <  2:3943334/150‑1 (MQ=255)
                                                                                                     cTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTaa                                                  <  2:313962/150‑1 (MQ=255)
                                                                                                      tGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAAc                                                 >  2:3190274/1‑150 (MQ=255)
                                                                                                          gcgcGTGGAATAGAGAGCGCCAGTTCACCTTTAGGTTTGTCGCGATAAGCCTTAGCTGATACATCATTCTTAATAAGATAAGATGTATTTTATATGCATCTTTAAGGCAAAAGCTATTACCCTTACGTAGCAGAGGCTTATTTTTCTCAc                                             <  2:1731659/150‑1 (MQ=255)
                                                                                                           cgcgAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACt                                            <  2:2675398/150‑1 (MQ=255)
                                                                                                              gAGGAATAGAGAGCGCCAGTTCACCTACAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGcc                                         <  2:3065189/150‑1 (MQ=255)
                                                                                                               aGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCa                                        >  2:1139569/1‑150 (MQ=255)
                                                                                                                     agagagCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCgtgt                                  <  2:4344070/150‑1 (MQ=255)
                                                                                                                       agagCGCCAGTTCACCTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCa                                >  2:4519916/1‑150 (MQ=255)
                                                                                                                             ccAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTc                          >  1:3776025/1‑150 (MQ=255)
                                                                                                                                gTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAAATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTga                       <  1:4197634/150‑1 (MQ=255)
                                                                                                                                gTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAAATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTga                       <  1:1699324/150‑1 (MQ=255)
                                                                                                                                  tCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTgaga                     >  1:3285656/1‑150 (MQ=255)
                                                                                                                                      cTAAAGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGa                 >  1:3192344/1‑150 (MQ=255)
                                                                                                                                      cTAAAGGTGGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGa                 >  2:4280400/1‑150 (MQ=255)
                                                                                                                                       tAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGAt                <  1:2086557/150‑1 (MQ=255)
                                                                                                                                         aaGGTTGGTCGCTATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTc              <  2:3421594/150‑1 (MQ=255)
                                                                                                                                         aaGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTc              >  1:3609370/1‑150 (MQ=255)
                                                                                                                                         aaGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTc              >  2:3842198/1‑150 (MQ=255)
                                                                                                                                          aGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTCa             >  1:5037959/1‑150 (MQ=255)
                                                                                                                                           ggTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTCAg            <  1:704617/150‑1 (MQ=255)
                                                                                                                                             ttGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTCAGtt          >  2:5063874/1‑150 (MQ=255)
                                                                                                                                             aaggTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCAGAGACAGATTCAGtt          >  2:696419/3‑150 (MQ=255)
                                                                                                                                              tGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTCAGTTa         <  1:126334/150‑1 (MQ=255)
                                                                                                                                                 tCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTCAGTTAAAg      <  2:2419998/150‑1 (MQ=255)
                                                                                                                                                    cgATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTCAGTTAAAGGtt   >  1:1619825/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGGCAAAAGCTATAACCCTTACGTAGCAGAGGCTTATTTAACTCACTGCCAGCGTGTCAACTTTCTGAGACAGATTCAGTTAAAGGTTT  >  CP009273/4421642‑4421940

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