Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 580,403 G→T 5.0% N229K (AAC→AAA ompT ← DLP12 prophage; outer membrane protease VII (outer membrane protein 3b)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273580,4030GT 5.0% 1185.6 / 24.6 298N229K (AAC→AAAompTDLP12 prophage; outer membrane protease VII (outer membrane protein 3b)
Reads supporting (aligned to +/- strand):  ref base G (144/139);  new base T (5/10);  total (149/149)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 302 total aligned reads displayed.

TATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGAT  >  CP009273/580255‑580552
                                                                                                                                                    |                                                                                                                                                     
  ttCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCTTTa                                                                                                                                                    <  2:4547175/150‑1 (MQ=255)
   tCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTAt                                                                                                                                                   >  2:2145201/1‑150 (MQ=255)
    cGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATc                                                                                                                                                  <  2:4995439/150‑1 (MQ=255)
    cGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATc                                                                                                                                                  <  2:1104626/150‑1 (MQ=255)
     gTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCa                                                                                                                                                 >  1:4841681/1‑150 (MQ=255)
     gTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCAGGTCATAGTGTTCATCGTTATCa                                                                                                                                                 >  2:3917337/1‑150 (MQ=255)
        aCCCGATTCCATGCGCCTTCAACCTAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGGCATAGTGTTCATCGTTATCAgat                                                                                                                                              <  1:447418/150‑1 (MQ=255)
         cccGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAgatg                                                                                                                                             >  1:596019/1‑150 (MQ=255)
         cccGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAgatg                                                                                                                                             >  1:2721210/1‑150 (MQ=255)
         cccGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAgatg                                                                                                                                             >  1:3274018/1‑150 (MQ=255)
         cccGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAgatg                                                                                                                                             >  2:2504996/1‑150 (MQ=255)
            gATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGAtt                                                                                                                                          >  2:3849707/1‑150 (MQ=255)
            gATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGAtt                                                                                                                                          >  1:3932864/1‑150 (MQ=255)
             aTTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTc                                                                                                                                         <  1:2437481/150‑1 (MQ=255)
             aTTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTc                                                                                                                                         <  1:3769949/150‑1 (MQ=255)
              ttCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTcc                                                                                                                                        >  1:1774078/1‑150 (MQ=255)
              ttCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTcc                                                                                                                                        >  1:1520259/1‑150 (MQ=255)
                ccATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCTTTATCAGATGATTCCAc                                                                                                                                      >  2:3258681/1‑150 (MQ=255)
                 cATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCAcc                                                                                                                                     >  1:3722380/1‑150 (MQ=255)
                     cgacTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGc                                                                                                                                 >  1:4694954/4‑150 (MQ=255)
                      gcCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGcc                                                                                                                                <  2:1329543/150‑1 (MQ=255)
                       ccTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCg                                                                                                                               >  2:4130010/1‑150 (MQ=255)
                       ccTTCAACATAAACTTTTGATTTATGTGTTACGTAATAACCAGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTAATCGTTATCAGATGATTCCACCCAGCCg                                                                                                                               <  1:2979830/150‑1 (MQ=255)
                        cTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGc                                                                                                                              >  1:2850174/1‑150 (MQ=255)
                            aaCATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTa                                                                                                                          <  1:3799491/150‑1 (MQ=255)
                            aaCATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTa                                                                                                                          <  1:3800002/150‑1 (MQ=255)
                             aCATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTAt                                                                                                                         <  1:3034182/150‑1 (MQ=255)
                             aCATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTAt                                                                                                                         >  2:3399155/1‑150 (MQ=255)
                              cATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTAtt                                                                                                                        >  2:642970/1‑150 (MQ=255)
                              cATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTAtt                                                                                                                        <  1:3601189/150‑1 (MQ=255)
                                 aaaCTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTaa                                                                                                                     <  1:4396861/150‑1 (MQ=255)
                                  aaCTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTaaa                                                                                                                    <  1:5093431/150‑1 (MQ=255)
                                  aaCTTTTGCGTTAGGTGTGACGTAATAACCTGCATAGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCAACGTTATCAGATGATTCCACCCAGCCGCTGTATTTaaa                                                                                                                    >  2:1190085/1‑150 (MQ=255)
                                      tttGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAAtgtg                                                                                                                <  2:1449929/150‑1 (MQ=255)
                                      tttGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAAtgtg                                                                                                                <  2:3588960/150‑1 (MQ=255)
                                       ttGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGc                                                                                                               >  2:2553234/1‑150 (MQ=255)
                                           gTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCAcc                                                                                                           >  1:4048295/1‑150 (MQ=255)
                                            ttAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCTTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCg                                                                                                          <  1:2626793/150‑1 (MQ=255)
                                            ttAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCg                                                                                                          >  2:4830388/1‑150 (MQ=255)
                                                gtgtGACGTAATAACCTGCATTGACTTCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTCCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGtt                                                                                                      <  1:3032816/150‑1 (MQ=255)
                                                     aCGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCaaaa                                                                                                 >  1:4549036/1‑150 (MQ=255)
                                                       gTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATc                                                                                               <  1:2497626/150‑1 (MQ=255)
                                                         aataaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCtt                                                                                             <  1:3845163/150‑1 (MQ=255)
                                                           taaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCa                                                                                           >  2:2997373/1‑150 (MQ=255)
                                                           taaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCa                                                                                           <  1:3123851/150‑1 (MQ=255)
                                                            aaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCAt                                                                                          <  2:181998/150‑1 (MQ=255)
                                                             aCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATa                                                                                         >  1:1800695/1‑150 (MQ=255)
                                                               cTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAAc                                                                                       <  2:2627056/150‑1 (MQ=255)
                                                                 gCATTGACTGCAACAGAATAGTATTTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAAGGa                                                                                     <  1:2165547/150‑1 (MQ=255)
                                                                        cTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTc                                                                              >  2:4535765/1‑150 (MQ=255)
                                                                         tGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTcc                                                                             >  1:2868637/1‑150 (MQ=255)
                                                                          gCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCa                                                                            <  1:1394176/150‑1 (MQ=255)
                                                                          gCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCa                                                                            <  2:713871/150‑1 (MQ=255)
                                                                           cAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAg                                                                           <  2:1966983/150‑1 (MQ=255)
                                                                           cAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAg                                                                           >  1:1138267/1‑150 (MQ=255)
                                                                              cAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCa                                                                        >  1:1181676/1‑150 (MQ=255)
                                                                              cAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCa                                                                        <  2:3782628/150‑1 (MQ=255)
                                                                               aGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCaa                                                                       >  2:1501340/1‑150 (MQ=255)
                                                                                 aaTAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGc                                                                     <  1:1318172/150‑1 (MQ=255)
                                                                                  aTAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCCCCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGcc                                                                    >  2:917179/1‑150 (MQ=255)
                                                                                    aGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCaa                                                                  <  1:5021974/150‑1 (MQ=255)
                                                                                     gTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAAt                                                                 <  2:3681512/150‑1 (MQ=255)
                                                                                      tAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATg                                                                <  1:917119/150‑1 (MQ=255)
                                                                                         ttttGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAg                                                             >  1:4225885/1‑150 (MQ=255)
                                                                                          tttGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAgg                                                            <  1:2084219/150‑1 (MQ=255)
                                                                                          tttGGTCTTTGACCTTACTGCGATAAGTGATTCTTTGTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAgg                                                            <  2:4162050/150‑1 (MQ=255)
                                                                                           ttGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAggg                                                           >  2:1025683/1‑150 (MQ=255)
                                                                                            tGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCTTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGc                                                          >  1:4091105/1‑150 (MQ=255)
                                                                                            tGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGc                                                          >  2:1225685/1‑150 (MQ=255)
                                                                                              gTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCAt                                                        >  1:538945/1‑150 (MQ=255)
                                                                                                cTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCAttt                                                      <  2:3118361/150‑1 (MQ=255)
                                                                                                cTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCAttt                                                      >  2:516811/1‑150 (MQ=255)
                                                                                                cTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCAttt                                                      >  2:448589/1‑150 (MQ=255)
                                                                                                cTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCAttt                                                      >  2:4475645/1‑150 (MQ=255)
                                                                                                 tttGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCAtttt                                                     >  2:921389/1‑150 (MQ=255)
                                                                                                      ccTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAAc                                                <  1:756164/150‑1 (MQ=255)
                                                                                                            tGCGATAAGTGATTCTTTTTCCCGGGTCAGAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACgttgtt                                          <  1:1831498/150‑1 (MQ=255)
                                                                                                                aTAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTa                                      <  1:1066545/150‑1 (MQ=255)
                                                                                                                aTAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTa                                      <  1:2315902/150‑1 (MQ=255)
                                                                                                                 tAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAg                                     <  1:5093369/150‑1 (MQ=255)
                                                                                                                 tAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAg                                     >  2:4369711/1‑150 (MQ=255)
                                                                                                                   aGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGcc                                   <  1:1663613/150‑1 (MQ=255)
                                                                                                                   aGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGcc                                   >  2:1696662/1‑150 (MQ=255)
                                                                                                                   aGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGcc                                   >  2:1257704/1‑150 (MQ=255)
                                                                                                                      gATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGAt                                <  1:1054776/150‑1 (MQ=255)
                                                                                                                       aTTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGAtt                               <  1:2989782/150‑1 (MQ=255)
                                                                                                                        ttCTTTTTCCCGGGTCATAGTGTTCATCTTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTg                              <  2:2867017/150‑1 (MQ=255)
                                                                                                                             tttCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCtt                         >  2:574968/1‑150 (MQ=255)
                                                                                                                             tttCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCtt                         >  2:1552414/1‑150 (MQ=255)
                                                                                                                              ttCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCttt                        >  1:4137369/1‑150 (MQ=255)
                                                                                                                               tCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTtc                       <  1:4624906/150‑1 (MQ=255)
                                                                                                                                 ccGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTtctc                     >  2:2366916/1‑150 (MQ=255)
                                                                                                                                   gggTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCa                   <  2:3972774/150‑1 (MQ=255)
                                                                                                                                    ggTCATTGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCAt                  <  1:3424858/150‑1 (MQ=255)
                                                                                                                                      tCATAGTGTTCATCTTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTc                <  1:1932796/150‑1 (MQ=255)
                                                                                                                                      tCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTc                >  2:2444931/1‑150 (MQ=255)
                                                                                                                                          aGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGa            <  1:3137059/150‑1 (MQ=255)
                                                                                                                                           gtgtTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGaa           <  2:4472458/150‑1 (MQ=255)
                                                                                                                                                 aTCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGcc     <  1:4010840/150‑1 (MQ=255)
                                                                                                                                                   cGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGa   <  2:2977240/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTTTGTAGCCGATTGCTCTTTCTCCATTCGGGAAGGAGCCGAT  >  CP009273/580255‑580552

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: