Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,113,450 A→T 100% I161I (ATT→ATA yceI ← secreted protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,113,4500AT99.1% 731.5 / ‑2.7 221I161I (ATT→ATAyceIsecreted protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base T (133/86);  total (135/86)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 221 total aligned reads displayed.

CAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTAATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGTCTTTC  >  CP009273/1113305‑1113594
                                                                                                                                                 |                                                                                                                                                
cAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTAt                                                                                                                                                >  1:3273160/1‑148 (MQ=255)
cAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTAt                                                                                                                                                <  1:842716/148‑1 (MQ=255)
 aaCAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTACTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTAtt                                                                                                                                               >  2:2431138/1‑148 (MQ=255)
  aCAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTAttt                                                                                                                                              >  2:1284557/1‑148 (MQ=255)
   cAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCATGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTAtttt                                                                                                                                             >  2:2678317/1‑148 (MQ=255)
   cAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTAtttt                                                                                                                                             >  2:2977838/1‑148 (MQ=255)
    aGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTg                                                                                                                                            <  1:2181186/148‑1 (MQ=255)
     gAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGc                                                                                                                                           >  2:1671905/1‑148 (MQ=255)
      aataaataTAATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGcc                                                                                                                                          <  1:173405/148‑1 (MQ=255)
        taaataTCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCtt                                                                                                                                        >  1:2164096/1‑148 (MQ=255)
          aataTCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCg                                                                                                                                      >  1:3089948/1‑148 (MQ=255)
           atatCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTAATTTTGCCTTCgg                                                                                                                                     >  2:864816/1‑148 (MQ=255)
             atCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGcc                                                                                                                                   <  2:335186/148‑1 (MQ=255)
              tCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCt                                                                                                                                  <  2:326396/148‑1 (MQ=255)
               cATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTc                                                                                                                                 >  2:2935068/1‑148 (MQ=255)
                 ttGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGa                                                                                                                               >  2:509392/1‑148 (MQ=255)
                 ttGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGa                                                                                                                               >  2:1511388/1‑148 (MQ=255)
                   gCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAg                                                                                                                             <  1:871954/148‑1 (MQ=255)
                    ccGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGc                                                                                                                            <  1:963499/148‑1 (MQ=255)
                    ccGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGc                                                                                                                            >  2:1590769/1‑148 (MQ=255)
                       gATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTg                                                                                                                         >  1:1180919/1‑148 (MQ=255)
                       gATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTg                                                                                                                         >  1:2931304/1‑148 (MQ=255)
                       gATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCGTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTg                                                                                                                         >  1:1855189/1‑148 (MQ=255)
                        aTGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGc                                                                                                                        >  2:349677/1‑148 (MQ=255)
                        aTGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGc                                                                                                                        >  2:2224959/1‑148 (MQ=255)
                         tGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCa                                                                                                                       >  2:3030261/1‑148 (MQ=255)
                            gTGCCATCCGGCAACATTTGACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGt                                                                                                                    <  2:884181/148‑1 (MQ=255)
                                      gCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCAcccc                                                                                                          >  1:2994509/1‑148 (MQ=255)
                                       cAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCa                                                                                                         >  2:657945/1‑148 (MQ=255)
                                       cAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCa                                                                                                         >  2:2270594/1‑148 (MQ=255)
                                        aaCATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCAt                                                                                                        >  1:2051364/1‑148 (MQ=255)
                                          cATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATgg                                                                                                      <  1:500193/148‑1 (MQ=255)
                                          cATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATgg                                                                                                      >  1:695650/1‑148 (MQ=255)
                                           aTTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATggg                                                                                                     >  1:3179082/1‑148 (MQ=255)
                                            tttCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGgt                                                                                                    <  1:1835761/148‑1 (MQ=255)
                                            tttCACGCTTACTTCAGCTGTACTCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTACGGCATCGAAGCCTGCACGTTTACCACCCCATGGgt                                                                                                    >  2:681020/1‑148 (MQ=255)
                                             ttCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGgtc                                                                                                   >  2:3185194/1‑148 (MQ=255)
                                             ttCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGgtc                                                                                                   >  2:1456445/1‑148 (MQ=255)
                                             ttCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGgtc                                                                                                   >  1:1762303/1‑148 (MQ=255)
                                             ttCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGgtc                                                                                                   >  2:3438623/1‑148 (MQ=255)
                                                    ttACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTATTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGccc                                                                                            <  2:3091408/148‑1 (MQ=255)
                                                       cTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGa                                                                                         >  2:2688107/1‑148 (MQ=255)
                                                         tCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGAcc                                                                                       >  2:2311336/1‑148 (MQ=255)
                                                             ctgTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCaatt                                                                                   >  1:600585/1‑148 (MQ=255)
                                                                tACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCaattaat                                                                                <  1:3243363/148‑1 (MQ=255)
                                                                  cGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAAttt                                                                              >  1:1376880/1‑148 (MQ=255)
                                                                  cGCCTTCCACTGAAATAATCAGATCCACTTACTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAAttt                                                                              <  1:2431138/148‑1 (MQ=255)
                                                                    ccTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAGGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCg                                                                            >  2:1443124/1‑148 (MQ=255)
                                                                    ccTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCg                                                                            >  2:1448583/1‑148 (MQ=255)
                                                                      ttCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCt                                                                          <  2:780805/148‑1 (MQ=255)
                                                                      ttCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCt                                                                          >  1:2456048/1‑148 (MQ=255)
                                                                        ccACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTc                                                                        <  2:174855/148‑1 (MQ=255)
                                                                           cTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAg                                                                     >  1:1293744/1‑148 (MQ=255)
                                                                               aataatCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTg                                                                 <  1:2318947/148‑1 (MQ=255)
                                                                               aataatCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTg                                                                 <  1:349677/148‑1 (MQ=255)
                                                                               aataatCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTg                                                                 <  2:1297159/148‑1 (MQ=255)
                                                                                  aatCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAAGTAATTTCGCTTCCAGCGTGACg                                                              >  2:504720/1‑148 (MQ=255)
                                                                                   atCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACgg                                                             >  2:1723353/1‑148 (MQ=255)
                                                                                   atCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGATTCCAGCGTGACgg                                                             >  2:2608476/1‑148 (MQ=255)
                                                                                    tCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACggg                                                            <  2:1253296/148‑1 (MQ=255)
                                                                                     cAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGt                                                           <  1:2935068/148‑1 (MQ=255)
                                                                                       gATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGttt                                                         >  2:253571/1‑148 (MQ=255)
                                                                                             cTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCACTTAATTTCGCTTCCAGCGTGACGGGTTTGGTTAc                                                   >  2:1425309/1‑148 (MQ=255)
                                                                                             cTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTAc                                                   >  2:1425124/1‑148 (MQ=255)
                                                                                              ttCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACg                                                  <  2:938002/148‑1 (MQ=255)
                                                                                                ccTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGcc                                                <  1:2130047/148‑1 (MQ=255)
                                                                                                 cTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCa                                               <  2:1069458/148‑1 (MQ=255)
                                                                                                  tGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCAt                                              >  1:1064718/1‑148 (MQ=255)
                                                                                                         cTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGt                                       <  1:2482432/148‑1 (MQ=255)
                                                                                                         cTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGt                                       >  2:1474815/1‑148 (MQ=255)
                                                                                                         cTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGt                                       <  1:2862218/148‑1 (MQ=255)
                                                                                                          tGGACCTAAATCTGTCTTGATATTGAAGTCCTTTAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTc                                      <  2:154332/148‑1 (MQ=255)
                                                                                                          tGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTc                                      >  1:223268/1‑148 (MQ=255)
                                                                                                          tGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTc                                      <  1:2678317/148‑1 (MQ=255)
                                                                                                             aCCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGa                                   <  2:2006759/148‑1 (MQ=255)
                                                                                                              ccTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGAt                                  >  2:2574153/1‑148 (MQ=255)
                                                                                                                tAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCg                                >  2:2371070/1‑148 (MQ=255)
                                                                                                                   aTCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCg                             >  1:641839/1‑148 (MQ=255)
                                                                                                                   aTCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCg                             >  2:2741179/1‑148 (MQ=255)
                                                                                                                   aTCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCCCCg                             >  1:3428993/1‑148 (MQ=255)
                                                                                                                      tGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTa                          >  2:1844390/1‑148 (MQ=255)
                                                                                                                       gTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTaa                         >  1:25036/1‑148 (MQ=255)
                                                                                                                          ttGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAatat                      <  2:1046454/148‑1 (MQ=255)
                                                                                                                            gATATTTAAGTCCTTGAGCTTAATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATcc                    >  2:3101747/1‑148 (MQ=255)
                                                                                                                            gATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTGTGGTTACGCCCGTCAGCGGCAGATCGCCGCTAATATcg                    >  2:1837471/1‑147 (MQ=255)
                                                                                                                             atatTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCa                   >  1:1984434/1‑148 (MQ=255)
                                                                                                                              tatTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAg                  >  2:2882336/1‑148 (MQ=255)
                                                                                                                                ttGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGtt                >  2:2778348/1‑148 (MQ=255)
                                                                                                                                  gAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCg              >  1:1008676/1‑148 (MQ=255)
                                                                                                                                  gAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCg              >  1:1009084/1‑148 (MQ=255)
                                                                                                                                   aaGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGt             >  1:1412659/1‑148 (MQ=255)
                                                                                                                                    aGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTc            <  1:1274350/148‑1 (MQ=255)
                                                                                                                                     gTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCa           >  2:1919725/1‑148 (MQ=255)
                                                                                                                                       ccTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCAcc         >  1:1118596/1‑148 (MQ=255)
                                                                                                                                         ttGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGt       <  2:3374546/148‑1 (MQ=255)
                                                                                                                                          tGAGCTTTAGTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACGCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGtc      <  2:1894331/148‑1 (MQ=255)
                                                                                                                                           gAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGtct     <  2:1631152/148‑1 (MQ=255)
                                                                                                                                           gAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGtct     <  2:2759910/148‑1 (MQ=255)
                                                                                                                                            aGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGtctt    <  1:1677125/148‑1 (MQ=255)
                                                                                                                                              cTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGtctttc  <  2:1678965/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTAATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGTCTTTC  >  CP009273/1113305‑1113594

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: