Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,412,385 G→T 100% N34K (AAC→AAA kilR ← Rac prophage; inhibitor of ftsZ, killing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,412,3850GT98.1% 834.3 / 5.3 260N34K (AAC→AAAkilRRac prophage; inhibitor of ftsZ, killing protein
Reads supporting (aligned to +/- strand):  ref base G (5/0);  new base T (128/127);  total (133/127)
Fisher's exact test for biased strand distribution p-value = 6.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 261 total aligned reads displayed.

AAGCACATTTTATTTTCTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTGTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACTTTGTTGATGCGGTGTC  >  CP009273/1412241‑1412529
                                                                                                                                                |                                                                                                                                                
 aGCACATTTTATTTTCTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTg                                                                                                                                              >  1:1020647/1‑148 (MQ=255)
  gCACATTTTATTTTCTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGc                                                                                                                                             <  1:2282482/148‑1 (MQ=255)
  gCACATTTTATTTTCTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGc                                                                                                                                             <  2:2111163/148‑1 (MQ=255)
   cacaTTTTATTTTCTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGcc                                                                                                                                            >  2:812522/1‑148 (MQ=255)
       ttttattttCTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAgg                                                                                                                                        >  2:3444978/1‑148 (MQ=255)
         ttattttCTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCa                                                                                                                                      <  1:1288697/148‑1 (MQ=255)
         ttattttCTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCa                                                                                                                                      <  1:1845249/148‑1 (MQ=255)
         ttattttCTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCa                                                                                                                                      <  1:698766/148‑1 (MQ=255)
           attttCTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAat                                                                                                                                    <  2:2659107/148‑1 (MQ=255)
           attttCTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAat                                                                                                                                    >  2:1820082/1‑148 (MQ=255)
            ttttCTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAata                                                                                                                                   <  2:2501843/148‑1 (MQ=255)
             tttCTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGACTGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAatat                                                                                                                                  >  1:3016315/1‑148 (MQ=255)
              ttCTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAatata                                                                                                                                 >  1:2564338/1‑148 (MQ=255)
               tCTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAatatat                                                                                                                                <  2:1561192/148‑1 (MQ=255)
                cTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATg                                                                                                                               >  2:608106/1‑148 (MQ=255)
                cTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATg                                                                                                                               >  1:3319496/1‑148 (MQ=255)
                    gTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCg                                                                                                                           >  2:2064633/1‑148 (MQ=255)
                        tcaCCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCt                                                                                                                       <  1:1942434/148‑1 (MQ=255)
                           ccATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCa                                                                                                                    <  1:873248/148‑1 (MQ=255)
                            cATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCAt                                                                                                                   >  2:55385/1‑148 (MQ=255)
                              tGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATga                                                                                                                 >  2:2011841/1‑148 (MQ=255)
                                aCTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATgaag                                                                                                               >  1:774707/1‑148 (MQ=255)
                                    cGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACa                                                                                                           >  1:2634442/1‑148 (MQ=255)
                                     gCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACAt                                                                                                          <  2:40429/148‑1 (MQ=255)
                                      ccTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATa                                                                                                         <  2:694657/148‑1 (MQ=255)
                                       cTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATaa                                                                                                        <  1:2478062/148‑1 (MQ=255)
                                       cTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATaa                                                                                                        <  1:2854941/148‑1 (MQ=255)
                                       cTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATaa                                                                                                        <  1:941147/148‑1 (MQ=255)
                                       cTTTACAGGTAAACCATCACGACCGAGGAAGACTATAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATaa                                                                                                        <  2:1352912/148‑1 (MQ=255)
                                                tAAACCATCACGACCGAGTAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGa                                                                                               <  2:1904456/148‑1 (MQ=255)
                                                   aCCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTc                                                                                            >  1:137304/1‑148 (MQ=255)
                                                      aTCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTcaca                                                                                         >  2:1510425/1‑148 (MQ=255)
                                                      aTCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTcaca                                                                                         <  1:832174/148‑1 (MQ=255)
                                                       tCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTcacac                                                                                        <  1:2313512/148‑1 (MQ=255)
                                                       tCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTcacac                                                                                        >  1:69228/1‑148 (MQ=255)
                                                        cACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTcacaca                                                                                       >  1:1195918/1‑148 (MQ=255)
                                                          cGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACt                                                                                     <  1:2018774/148‑1 (MQ=255)
                                                            aCCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGa                                                                                   >  2:648317/1‑148 (MQ=255)
                                                                aGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTg                                                                               >  2:3162302/1‑148 (MQ=255)
                                                                aGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTg                                                                               >  1:1226675/1‑148 (MQ=255)
                                                                  gAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGATGTGGt                                                                             >  1:2148388/1‑148 (MQ=255)
                                                                   aaGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGtc                                                                            >  2:3233962/1‑147 (MQ=255)
                                                                    aGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTAt                                                                           >  1:99904/1‑148 (MQ=255)
                                                                      aCTTTAATCATGCGGTCAGTAATTAATGTTTTTGTGGTCAGGTTACTAATATATAGTTTTCGATTTTTAATATTTTTTGCCGAGGAAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTAttt                                                                         <  1:1510023/148‑1 (MQ=255)
                                                                       cTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTTGTATTtc                                                                        <  2:1665818/148‑1 (MQ=255)
                                                                         ttAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTtcat                                                                      >  2:2584633/1‑148 (MQ=255)
                                                                          tAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATACTCGCCAGGAGTCACACACTGACGTGGTATTtcatc                                                                     <  1:3301445/148‑1 (MQ=255)
                                                                            aTCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAg                                                                   >  2:127942/1‑148 (MQ=255)
                                                                               aTGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTc                                                                <  2:2573997/148‑1 (MQ=255)
                                                                                    gTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAg                                                           <  2:274185/148‑1 (MQ=255)
                                                                                      cAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTg                                                         <  1:1110390/148‑1 (MQ=255)
                                                                                      cAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTg                                                         <  1:1909392/148‑1 (MQ=255)
                                                                                      cAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTg                                                         <  1:369459/148‑1 (MQ=255)
                                                                                        gTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGAt                                                       <  1:3466335/148‑1 (MQ=255)
                                                                                         tAATGAATTTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGAtg                                                      <  1:2340294/148‑1 (MQ=255)
                                                                                           atgaatgTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGAtgtg                                                    >  1:1001987/1‑148 (MQ=255)
                                                                                             gaatgTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCa                                                  >  2:608874/1‑148 (MQ=255)
                                                                                             gaatgTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCa                                                  <  1:681184/148‑1 (MQ=255)
                                                                                             gaatgTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCa                                                  >  2:453107/1‑148 (MQ=255)
                                                                                             gaatgTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCa                                                  <  1:575547/148‑1 (MQ=255)
                                                                                                tgTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATc                                               <  1:2343544/148‑1 (MQ=255)
                                                                                                tgTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATc                                               <  1:2205935/148‑1 (MQ=255)
                                                                                                 gTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCa                                              >  2:3053136/1‑148 (MQ=255)
                                                                                                  tttttGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCAt                                             >  1:211022/1‑148 (MQ=255)
                                                                                                   ttttGTGGTCAGGTTACGAATATATATTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATa                                            <  2:561465/148‑1 (MQ=255)
                                                                                                    tttGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATaa                                           <  1:2064633/148‑1 (MQ=255)
                                                                                                    tttGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATaa                                           <  1:3076263/148‑1 (MQ=255)
                                                                                                        tgGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTAt                                       >  2:3395502/1‑148 (MQ=255)
                                                                                                             aGGTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCa                                  <  2:2888473/148‑1 (MQ=255)
                                                                                                              ggTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCAt                                 >  2:217402/1‑148 (MQ=255)
                                                                                                              ggTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCAt                                 <  1:2904171/148‑1 (MQ=255)
                                                                                                              ggTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCAt                                 >  2:2978108/1‑148 (MQ=255)
                                                                                                               gTTACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCAtt                                <  2:1737195/148‑1 (MQ=255)
                                                                                                                ttACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCAttt                               >  1:2782453/1‑148 (MQ=255)
                                                                                                                ttACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCAttt                               >  2:1842015/1‑148 (MQ=255)
                                                                                                                 tACGAATATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCAtttt                              >  2:2640690/1‑148 (MQ=255)
                                                                                                                    gAATATATAGTTTTCGCTTTTTAAAATTTTTTGCAGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGAGCTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTAc                           >  2:1109320/1‑148 (MQ=255)
                                                                                                                     aaTATATAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTAACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACa                          >  2:429825/1‑148 (MQ=255)
                                                                                                                      atatatAGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACaa                         >  2:2492924/1‑148 (MQ=255)
                                                                                                                       tatataGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACaaa                        <  2:3318930/148‑1 (MQ=255)
                                                                                                                          ataGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGa                     <  2:1627766/148‑1 (MQ=255)
                                                                                                                           taGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGaa                    <  1:3391519/148‑1 (MQ=255)
                                                                                                                            aGTTTTCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAAc                   <  1:1786829/148‑1 (MQ=255)
                                                                                                                              ttttCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACtt                 <  1:3244450/148‑1 (MQ=255)
                                                                                                                              ttttCGCTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACtt                 >  2:2244235/1‑148 (MQ=255)
                                                                                                                                   gCTTTTTAATATTTTGTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACTttgttg            <  2:3432576/148‑1 (MQ=255)
                                                                                                                                    cTTTTTAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACTTTGTTGa           >  2:1392463/1‑148 (MQ=255)
                                                                                                                                    cTTTTTAATATTTTTTGCCGAGGCAATATATGGCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACTTTGTTGa           <  1:53043/148‑1 (MQ=255)
                                                                                                                                      ttttAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACTTTGTTGATg         <  2:1226675/148‑1 (MQ=255)
                                                                                                                                       tttAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACTTTGTTGATga        >  2:3181083/1‑147 (MQ=255)
                                                                                                                                       tttAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACTTTGTTGATGc        >  2:2445231/1‑148 (MQ=255)
                                                                                                                                       tttAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACTTTGTTGATGc        >  1:1477481/1‑148 (MQ=255)
                                                                                                                                       tttAATATTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACTTTGTTGATGc        >  1:2535159/1‑148 (MQ=255)
                                                                                                                                           atatTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACTTTGTTGATGCGgtg    <  2:1028445/148‑1 (MQ=255)
                                                                                                                                           atatTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACTTTGTTGATGCGgtg    <  2:334530/148‑1 (MQ=255)
                                                                                                                                           atatTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACTTTGTTGATGCGgtg    <  1:63241/148‑1 (MQ=255)
                                                                                                                                           atatTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACTTTGTTGATGCGgtg    <  1:62278/148‑1 (MQ=255)
                                                                                                                                           atatTTTTTGCCGAGGCAATATATGTCCGGCCTTAATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACTTTGTTGATGCGgtg    <  2:2782453/148‑1 (MQ=255)
                                                                                                                                           atatTGTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACTTTGTTGATGCGgtg    >  1:142066/1‑148 (MQ=255)
                                                                                                                                             atTTTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACTTTGTTGATGCGGTGTc  >  2:931941/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AAGCACATTTTATTTTCTTAGTCATCACCATGACTCCGCCTTTACAGGTAAACCATCACGACCGAGGAAGACTTTAATCATGCGGTCAGTAATGAATGTTTTTGTGGTCAGGTTACGAATATATAGTTTTCGCTTTTTAATATTGTTTGCCGAGGCAATATATGTCCGGCCTTCATGAAGAACATAATCGCCAGGAGTCACACACTGACGTGGTATTTCATCAGTTCCGAAGTGATGTGCAATCATAATTATCTCCATTTTTACAAATGAACTTTGTTGATGCGGTGTC  >  CP009273/1412241‑1412529

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: